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This  investigation  was  the  last  part  of  a  program  conducted 
for  the  purpose  of  developing  a  predictive  model  of  heat  flow  and 
solidification  of  aluminum  alloys  produced  under  conditions  of  high 
cooling  rates.  It  addresses  the  characterization  of  the  thermal 
history  of  rapid  solidification  in  metal  droplets,  and  its  effects 
on  powder  microst ructure .  The  main  effort  was  focused  on  modelling 
the  heat  flow  during  solidification,  and  relationships  were 
established  between  the  atomization  parameters,  the  growth 
kinetics,  the  interface  velocity  and  undercooling,  and  other 
important  variables.  Numerical  solutions  based  on  the  enthalpy 
model  were  developed,  and  their  results  compared  to  the  trends 
predicted  from  the  Newtonian  model.  The  analysis  covered 
situations  of  isothermal  solidification  at  the  melting  temperature, 
as  well  as  those  where  significant  undercoolings  are  necessary  for 
nuc leaf  ion  and  growth.  The  implications  of  single  vs.  multiple 
nucleation  were  also  discussed. 

In  general,  the  results  indicate  that  when  substantial 
undercoolings  are  achieved  in  the  droplet  prior  to  nucleation,  the 
thermal  history  consists  of  two  distinct  solidification  regimes, 
in  the  first  one  the  interface  velocities  are  high,  the  droplet 
absorbs  most  of  the  latent  heat  released,  and  the  external  cooling 
usually  pLays  a  minor  role.  The  second  regime  is  one  of  slower 
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growth,  and  strongly  depends  on  the  heat  extraction  at  the  droplet 
surface.  The  extent  of  "rapid  solidification",  as  determined  from 
the  fraction  of  material  solidified  above  a  certain  critical 
undercooling,  is  a  function  of  the  nucleation  temperature,  the 
particle  size,  a  kinetic  parameter  and  the  heat  transfer 

coefficient.  Substantial  departures  from  the  Newtonian  model  were 

“4 

calculated  for  Biot  numbers  as  low  as  10  . 

The  concepts  developed  from  the  heat  flow  analysis  were 
coupled  to  microstructural  observat io  is  in  aluminum  alloy  powders, 
mostly  in  the  submicron  size  range.  It  was  shown  that'  reducing  the 
particle  size  decreases  the  extent  of  segregation,  promotes 
multiple  nucleation  and  the  formation  of  twins  during 
solidification.  Homogeniously  solidified  powders  were  found  in 
Al-370Si  and  670Si  below  1  pm,  in  diameter;  their  incidence  is 
enhanced  by  reducing  the  solute  content  and  increasing  the  heat 
transfer  coefficient. 
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CHAPTER  1 

INTRODUCTION 

Atomization  is  perhaps  the  process  with  the  highest 
potential  for  large  scale  production  of  rapidly  solidified 
metals.  However,  early  methods  for  generation  of  metal  powders 
were  limited  in  their  achievable  cooling  rates  to  about 
104  K  s-i,  and  their  benefits  in  terms  of  property  enhancement 
were  mainly  due  to  microstructural  refinement  and  homogeneous 
distribution  of  second  phases  and  inclusions.  Extensive  research 
on  splat  cooling  and  related  techniques,  where  cooling  rates  are 
estimated  to  range  from  10°  K  s  to  10  K  s  ,  has  led  to  a 
variety  of  metastable  modifications  in  alloy  microstructure. 

These  results  broaden  the  scenario  for  alloy  development,  for 
example,  through  extension  in  solid  solubility  of  elements  in  the 
primary  phase,  as  in  the  cases  of  highly  ductile  aluminum-silicon 
alloys  and  precipitation  hardening  aluminum-iron  alloys.  An 
important  goal  in  metal  powder  production  is  then  to  increase  the 
achievable  cooling  rates,  and  newer  processes  like  centrifugal 
atomization  are  claimed  to  surpass  10^  K  s-"1"  by  improvements  in 
the  droplet-gas  heat  transfer  coefficient,  resulting  from  high 
droplet  velocities  and  controlled  particle  size  in  the  range  of 
100  um  or  less.  However,  in  order  to  optimize  the  capabilities 
of  atomization  it  is  necessary  to  understand  the  solidification 
mechanism  in  small  metal  droplets  and  its  relationship  to 
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progressively  higher  rates  of  heat  extraction. 

There  is  also  a  theoretical  interest  in  solidification  of 
metal  powders.  Atomization  is  unique  among  the  rapid  cooling 
processes  because  the  interactions  with  molds,  substrates  or 
other  solidifying  particles  are  practically  absent.  Conceivably, 
then,  most  of  the  liquid  droplets  could  achieve  high  undercooling 
levels  prior  to  nucleation--provided  the  purity  level  is 
adequate-~because  the  active  nucleation  catalysts  will  be 
isolated  in  a  small  fraction  of  the  total  volume.  In  addition  to 
this  undercooling,  the  very  efficient  heat  transfer  is  expected 
to  partially  relieve  the  recalescence  effects  and  increase  the 
driving  force  for  microstr uctural  modifications.  Furthermore, 
the  development  of  electrohydrodynamic  atomization  has  provided  a 
tool  for  the  controlled  generation  of  powders  below  1  pm  in 
diameter.  These  powders,  while  not  of  commercial  importance,  are 
especially  suited  for  fundamental  studies  because  they  are 
electron-transparent ,  and  therefore  permit  the  three-dimensional 
observation  of  a  complete  solidification  structure. 

It  is  evident  from  the  current  literature  that  the 
microstructures  resulting  from  rapid  solidification  of  alloy 
powders  can  not  always  be  explained  with  available  theories  of 
nonequilibrium  dendritic  growth.  The  purpose  of  this 
investigation  was  to  enhance  our  understanding  of  the 
relationship  between  the  process  variables  in  atomization,  the 
undercooling  and  the  growth  kinetics,  and  the  powder 
microstructure  and  microchemistry. 
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There  were  two  main  areas  of  interest  in  the  scope  of  this 
investigation.  The  first  one  is  basically  theoretical  and  was 
aimed  to  modelling  the  heat  transfer  during  cooling  and 
solidification  of  atomized  metal  droplets.  The  second  is 
experimental  and  involved  the  microstr uctural  characterization  of 
aluminum  alloy  powders  produced  by  two  different  atomization 
techniques.  The  formation  of  these  structures  is  discussed  in 
terms  of  the  heat  flow  results  and  the  current  understanding  of 
the  solidification  mechanisms. 
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CHAPTER  2 

LITERATURE  REVIEW 

In  the  search  for  improved  properties  in  aluminum  alloys  a 
variety  of  techniques  which  involve  the  direct  quenching  from  the 
liquid  state  have  been  investigated.  The  general  approach  in 
most  of  them  has  been  to  reduce  the  volume  or  thickness  of 
material  being  solidified  at  one  time.  A  second  goal  in  their 
development  has  been  to  increase  the  heat  extraction  efficiency 
between  the  "casting"  and  the  "heat  sink".  The  overall  effect  is 
the  achievement  of  cooling  rates  estimated  in  the  range  of 
10-*  to  109  K  s-1. 

Solidification  under  the  increasingly  higher  cooling  rates 
— and  conceivably  higher  undercoolings — typical  of  atomization, 
splat  cooling  or  any  of  their  related  techniques,  commonly 
results  in  progressive  departure  from  the  microstructures 
produced  by  conventional  casting.  The  first  part  of  this  review 
summarizes  the  structural  modifications  reported  for  aluminum 
alloys  processed  by  rapid  solidification  (RSP) ,  whereas  the 
second  part  is  devoted  to  the  thermodynamic  and  kinetic 
considerations  involved  in  the  current  interpretation  of  these 
phenomena . 
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2.1  Metastable  Effects  in  Rapidly  Solidified  Aluminum  Alloys 

Rapid  solidification  of  aluminum  alloys  commonly  results  in 
one  or  more  of  the  following  modifications: 

1.  microstructural  refinement,  manifested  as  smaller  grain 
size  and  dendrite  arm  spacings; 

2.  extension  in  terminal  solid  solubility  of  the  primary 
u-Al  phase; 

3.  morphological  changes  of  the  eutectic  or  the  primary 
phase ; 

4.  formation  of  non-equilibrium  phases; 

5.  coupled  eutectic  growth  at  off-eutectic  compositions; 
and 

6.  vacancy  super saturation. 

Although  all  these  effects  have  been  observed  in  aluminum 
alloys,  some  even  at  the  moderate  cooling  rates  of  10  to 

O  _  I 

10  K  s  ,  only  the  first  four  are  reviewed  herein  because  they 
are  the  most  important  in  terms  of  improving  the  properties  of 
the  final  product. 

2.1.1  Microstructural  refinement 

In  general,  the  fineness  of  dendritic  microstructures  can  be 
correlated  to  the  average  cooling  rate  during  solidif ication,  r, 
or  the  local  solidification  time  (i.e.  time  available  for 
coarsening),  tf,  by: 


where  (TL  -  Tg)  is  the  solidification  temperature  range,  DAS  is 
the  dendrite  arm  spacing,  a  and  b  are  constants.  Relationships 
of  this  type  have  been  established  for  a  series  of  aluminum 
alloys  [1-4],  although  experimental  difficulties  have  usually 
limited  the  range  of  accurately  measured  cooling  rates  to  less 

-3  _1  + 

than  10  K  s  .  For  example,  studies  in  Al-4.5%Cu  at  cooling 
-5  3  -1 

rates  between  10  and  10  K  s  yield  values  for  a  =  40  and 
b  =  0.39  [2,3].  Values  for  a  =  47  and  b  =  0.33  have  also  been 
determined  for  thin  castings  of  Al-10.5%Si  at  cooling  rates  up  to 
105  K  s*1  [4]  . 

Dendritic  spacings  of  0.01  to  0.5  um  [5-8]  and  eutectic 
spacings  in  the  same  order  of  magnitude  [9]  have  been  reported 
for  electron  transparent  areas  of  A1  alloy  gun  splats.  When 
these  data  were  plotted  versus  the  corresponding  average  cooling 
rates  estimated  from  heat  flow  calculations ,  they  were  found  in 
reasonable  agreement  with  the  extrapolation  of  Equation  (2.1)  for 
a  and  b  established  at  much  lower  cooling  rates  [3,5,6].  These 
results  gave  then  origin  to  the  common  practice  of  estimating 
cooling  rates  based  on  dendrite  arm  spacing  measurements  in 
rapidly  solidified  materials.  It  should  be  noted,  however,  that 
caution  must  be  exercised  in  following  this  practice,  since  the 
dendritic  morphologies  change  sometimes  markedly  as  cooling  rate 
increases.  Furthermore,  the  constants  in  Equation  (2.1)  vary 
with  alloy  system  and  composition,  although  there  is  some 

^weight  %  unless  otherwise  specified 
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disagreement  in  the  literature  as  to  the  influence  of  this  latter 
variable.  Studies  on  the  effect  of  alloy  composition  on  structure 
at  fixed  average  cooling  rate  for  A1  binary  alloys  with  Cu,  Mg, 

Si  or  Zn  have  shown  that  increasing  the  solute  content  refines 
the  dendrite  arm  spacing,  especially  at  low  solute  concentrations 
[1,10]  . 

Refinement  of  secondary  phases  including  eutectic 
constituents  also  occurs  with  increasing  cooling  rate  during 
solidification.  Dispersion  hardening  associated  with  this  effect 
has  been  reported  in  Al-Si  [11]  and  Al-Fe  [12]  base  alloys. 

2.1.2  Extension  of  solid  solubility 

Extension  of  the  equilibrium  terminal  solid  solubility  has 
been  obtained  at  sufficiently  high  cooling  rates  in  almost  all 
the  aluminum  alloy  systems  investigated,  Al-Zn  being  the  only 
notable  exception  [13,14].  Table  1  shows  the  maximum 
super saturations  reported  for  several  elements  in  aluminum. 

Increased  solid  solubility  via  rapid  solidification  permits 
higher  alloying  levels  resulting  in  superior  strength-ductility 
combinations  in  a  number  of  systems.  Enhanced  solid  solution 
hardening  and  dispersion  strengthening  after  aging  have  been 
reported  for  Al-Cu  compositions  beyond  the  equilibrium  solubility 
limit  [15],  The  possibility  of  precipitation  hardening  in  Al-Fe 
base  alloys  has  been  demonstrated  [16] ,  and  considerable  interest 
exists  in  the  substantial  ductilities  attainable  in  high-silicon 
aluminum  alloys  [17,18]. 
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Table  1 

Extension  of  Solid  Solubility  in  Binary  Aluminum 


Alloys 

Quenched 

from  the  Melt 

[14] 

Solute 

Maximum  at 

Equilibrium 

RSP  Data 

te<k> 

X  g  ( a  t  % ) 

(at%) 

Cr 

934 

0.44 

5-6 

Cu 

821 

2.48 

17-18 

Fe 

928 

0.025 

4-6 

Mg 

724 

16.3 

36.8-40 

Mn 

932 

0.9 

6-9 

Ni 

913 

0.023 

1.2-7. 7 

Si 

850 

1.59 

10-16 

Zn 

655 

66.5 

38 

There  are  indications  that  the  extent  of  solute 
super saturation  increases  with  the  severity  of  the  quench  and 
increased  nominal  alloy  composition  [14,15].  Furthermore,  it  has 
been  shown  that  additions  especially  prone  to  solubility 
extension  in  binary  alloys  (e.g.  Mn)  can  be  used  in  ternaries  to 
increase  the  solubility  limit  of  less  susceptible  elements  (e.g. 
Fe,  Co,  Ni)  [19]. 

2.1.3  Morphological  modifications 

A  variety  of  grain  shapes  and  dendrite  morphologies  have 
been  reported  in  alloys  quenched  from  the  liquid  state  [20,21]. 
The  primary  phase  dendrites  in  some  alloys  tend  to  assume  a 
cellular  or  rod-like  configuration  with  progressively  faster 
cooling  rates  [6,21].  Similar  observations  have  been  reported 
for  bulk  specimens  of  nickel  and  iron  base  alloys  with  increasing 
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supercooling  prior  to  nucleation  of  the  solid  phase  [22],  In 
addition,  rapid  solidification  can  sometimes  suppress  formation 
of  an  equilibrium  primary  phase  from  the  melt.  For  example, 
formation  of  a-Al  dendrites  instead  of  equilibrium  Al^Fe  phase 
has  been  reported  in  hypereutectic  Al-Fe  alloys  both  in 
chill-cast  [23]  and  splat-cooled  specimens  [24] . 

Rapid  solidification  may  also  alter  the  composition  range  of 
cooperative  eutectic  growth,  the  eutectic  constituents  and  their 
morphologies.  For  example,  quenching  liquid  Al-Fe  alloys  in  the 
range  of  10  to  10^  K  s--*-  has  resulted  in  eutectic  morphologies 
changing  from  irregular  Al-Al3pe  to  regular  Al-Al6Fe,  as  well  as 
coupled  eutectic  growth  at  noneutectic  compositions  [23,25], 
Changes  in  the  interdendr i tic  eutectic  from  lamellar  to  rod-like 
have  been  shown  to  occur  at  high  cooling  rates  [24] .  In 
addition,  degenerate  and  radial  eutectic  structures  in 
Al-17.3  at%Cu  and  structural  changes  as  a  function  of  Cu  content 
have  been  studied  via  splat  cooling  [26] . 

2.1.4  Nonequilibrium  phases 

Formation  of  nonequilibrium  phases  during  rapid 
solidification  of  aluminum  alloys  is  not  as  common  as  the  effects 
noted  above.  Nevertheless,  a  certain  number  of  cases  are 
reported  in  the  literature  and  have  been  classified  by  Jones  [14] 
in  three  types: 
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I  -  Phases  existing  at  equilibrium  but  not  stable  at  the 

temperature  and  alloy  compositions  in  which  they  were 
observed.  Examples  include  the  formation  of  Al^Mn  and 
Al^Cr  in  alloys  where  the  phases  in  equilibrium  with 
«-Al  are  AlgMn  and  Al7Cr  [27,28]. 

II  -  Phases  not  formed  from  the  melt,  which  either  appear 

during  quenching  through  the  solid  state,  or  are 
precipitated  upon  further  heat  treatment.  Aging  of 
rapidly  solidified  Al-Fe  alloys  shows  a  sequence  of 
precipitation  similar  to  the  Al-Cu  system  [29,30]. 
Supersaturated  Al-Si  solid  solutions  also  exhibit  a 
gradual  decomposition,  starting  with  GP  zones  which 
are  nucleated  by  clustering  of  Si-vacancy  pairs  [31]. 
Although  some  vacancies  are  annihilated  to  accomodate 
the  misfit  between  the  matrix  and  the  cluster, 
coherent  GP  zones  do  not  grow  appreciably  larger  than 
1.5  to  2  nm  [32].  Subsequent  precipitation  has  been 
reported  to  produce  spherical,  acicular  and  lamellar 
morphologies,  as  well  as  rods  and  distorted  cubes, 
depending  on  the  alloy  composition  and  the  thermal 
history  [33-36] .  Si  needles  and  plates  commonly  form 
on  the  {ill}  and  { 100 }  planes  of  the  A1  parent  crystal 
[32]  . 

III  -  Phases  resulting  from  rapid  solidification,  but  not 

known  to  exist  under  equilibrium  conditions  at  any 
other  composition  within  the  alloy  system.  Typical  of 
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this  type  are  the  formation  of  Al^Fe  in  contact  with 
n-Al,  replacing  the  equilibrium  phase  Al^Fe  [37]. 
AI3CU2  has  been  produced  by  splat  quenching 
Al-45at%Cu,  whose  constituents  at  equilibrium  are  0 
(CuAl2)  and  n2  (~CuAl)  [38].  Furthermore, 
non-crystalline  structures  have  been  reported  in 
Al-17 . 3at%Cu  [9],  Al-30at%Ge  [39],  Al-Si  and  Al-Pd 
alloys  [40]. 

It  is  commonly  observed  that  particular  nonequilibrium 
effects  do  not  happen  above  a  critical  thickness  in  splat  cooling 
[41],  i.e.,  below  a  given  cooling  or  growth  rate  during 
solidification.  For  example,  Al6pe  can  displace  the  equilibrium 
phase  Al-jFe  at  cooling  rates  as  low  as  3  K  s  however,  cooling 
rates  of  the  order  of  104  K  s-1  are  required  to  avoid  formation 
of  FCC  Al^Mg2  in  Al-Mg  alloys  [14] . 

2.2  Thermodynamic  and  Kinetic  Considerations 

In  normal  solidification  processes  it  is  usually  assumed 
that  the  thermodynamic  conditions  at  the  interface,  i.e. 
temperature  and  composition,  correspond  fairly  closely  to  those 
dictated  by  the  equilibrium  phase  diagram,  and  that  the  growth 
rate  is  mainly  controlled  by  the  dissipation  of  the  latent  heat 
away  from  the  moving  front.  The  results  of  theories  assuming 
local  equilibrium  at  the  interface,  which  can  predict  the 
nonequilibrium  solidification  path  and  phase  composition  in  the 
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resulting  microstructures,  have  been  considered  as  supportive  of 
the  cited  assumption.  However,  when  dealing  with  processes  where 
the  heat  extraction  mechanisms  are  very  efficient  and  the  volumes 
of  material  involved  are  small  (splat  cooling  or  atomization), 
the  theories  based  on  interface  equilibrium  do  not  always 
satisfactorily  explain  the  metastable  structures  obtained,  and 
kinetic  considerations  must  be  brought  in  to  understand  the 
solidification  process.  The  following  paragraphs  summarize  the 
basic  thermodynamic  and  kinetic  concepts  involved  in  the  current 
interpretation  of  rapid  solidification  phenomena.  Since  the 
extension  of  solid  solubility  is  of  special  interest  to  the 
field,  part  of  this  review  is  devoted  to  that  particular  topic. 

2.2.1  Nucleation 

Solidification  under  rapid  cooling  rates  is  normally  a 
nucleation  and  growth  process,  except  in  the  formation  of 
amorphous  metals.  This  latter  phenomenon  is  usually  ascribed  to 
the  existence  of  a  large  energy  barrier  for  the  atoms  to  join  an 
embryo  (i.e.  the  viscosity  is  high  or  the  diffusivity  is  low), 
therby  suppressing  nucleation  until  the  metal  reaches  its  glass 
transition  temperature  [42] .  Recent  investigations  [43]  have 
shown  that  amorphous  solidification  can  be  obtained  in  some 
eutectic  alloys  even  in  the  presence  of  the  crystalline  phase  (no 
need  for  nucleation) .  If  the  alloy  is  thermodynamically  unable 
to  form  a  single  phase,  and  the  solidification  velocity  is 
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increased  enough  to  overcome  the  diffusion-controlled 
partitioning  process  at  the  interface,  the  growing  eutectic 
exhibits  a  transition  to  glass.  Nevertheless,  amorphous 
solidification  is  not  of  common  occurrence  in  the  atomization 
processes  and  alloys  of  interest  to  this  work,  and  therefore  will 
not  be  further  discussed. 

Since  the  growth  rate  is  considerable  at  relatively  small 
undercoolings  for  metals  and  other  liquids  of  high  fluidity, 
nucleation  is  usually  the  controlling  step  in  the  crystallization 
process.  Although  a  number  of  theories  for  homogeneous  and 
heterogeneous  nucleation  based  on  the  classical  droplet  model 
have  been  developed  and  discussed  elsewhere  [44] ,  their 
application  is  limited  by  the  lack  of  data,  especially  regarding 
the  surface  energies.  Classical  models  usually  assume  isothermal 
conditions  and  a  steady  state  distribution  of  embryos  and  their 
applicability  to  condensed  systems  has  been  questioned  in  the 
past.  Turnbull  [45]  proposed  the  concept  of  an  "incubation 
period"  (whose  length  depends  on  the  diffusivity)  for  the  steady 
state  distribution  of  embryos  to  develop  when  the  sample  is 
rapidly  cooled  from  above  to  below  the  melting  temperature. 

Reiss  and  Katz  [46]  have  analyzed  the  nucleation  problem  during 
continuous  cooling,  and  assuming  that  the  liquid  self-diffusion 
coefficient,  D^,  obeys  an  Ar rhenius- type  law,  they  concluded  that 
cooling  rates  in  excess  to  10  K  s  are  necessary  to  avoid  the 
steady  state  distribution  in  liquid  copper.  However,  they  also 
claim  that  the  actual  value  of  DL  may  be  actually  more  than  six 
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orders  of  magnitude  smaller  than  the  value  predicted  by  the 
Arrhenius  extrapolation  and,  therefore,  cooling  rates  in  excess 
of  10^  K  s-1  would  be  large  enough  to  produce  significant 
deviations  of  the  nucleation  rate  from  that  predicted  by  the 
steady  state  model.  In  general,  for  nucleation  to  be  suppressed 
until  a  desired  undercooling  is  reached  in  an  assembly  of  volume 
V,  the  total  number  of  nuclei  N,  given  by 

t 

N  =  V  /  J(t)  dt  (2.2) 

0 

should  be  less  than  unity.  In  this  expression  the  nucleation 
frequency  J  is  a  function  of  time  t  due  to  the  variations  in 
temperature  and  cooling  rate  during  quenching. 

Another  factor  which  may  be  questionable  regarding  the 
applicability  of  classical  expressions  to  nucleation  in  small 
droplets  is  the  number  of  atoms  available  to  produce  and  maintain 
the  statistical  distribution  of  embryos. 

Experimental  programs  using  either  an  emulsion  of  droplets 
[47-51],  or  single  droplets  [52]  as  a  means  of  isolating  the 
heterogeneous  nucleation  catalysts  in  a  small  fraction  of  the 
liquid,  have  been  aimed  at  determining  the  maximum  undercooling 
ATn  (conceivably  the  onset  of  homogeneous  nucleation)  achievable 
in  a  metal  or  alloy.  Hollomon  and  Turnbull  [48]  found  that  most 
metals  investigated  reach  a  maximum  undercooling  of  about 
0.18  Tm  (e.g.,  they  obtained  a  ATn  for  aluminum  of  195  K  or 
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0-21  Tm).  More  recently  Perepezko  et  al.  [49]  found 
undercoolings  for  Sn  and  Bi  in  the  range  of  0.35  and  therefore 
indicated  that  the  previously  reported  values  (and  perhaps  even 
these  latter  ones)  correspond  to  heterogeneous  and  not 
homogeneous  nucleation.  Furthermore,  it  has  been  shown  that 
reducing  the  droplet  size  markedly  increases  the  achievable 
undercooling  [50].  Droplet  experiments  in  the  Sn-Bi  [49]  and 
Cu-Ni  [51]  systems  revealed  that  undercoolings  comparable  to 
those  of  the  pure  metals  can  be  obtained  in  alloys  throughout  the 
entire  phase  diagram. 

Undercooling  results  have  been  used  in  the  past  to  estimate 
the  values  for  the  liquid-solid  interfacial  energies  from 
classical  nucleation  theory  [48].  However,  the  data  obtained 
will  always  be  underestimated  unless  homogeneous  nucleation 
conditions  are  satisfactorily  established  [49]. 

Christian  [42]  has  examined  the  compositional  variations 
during  the  nucleation  of  solid  from  a  liquid  alloy.  Using  the 
free  energy  (F)  vs.  composition  (X)  diagram  he  indicated  that 
the  maximum  driving  force  for  a  solid  of  composition  X'  to 
nucleate  from  a  liquid  of  composition  Xg  corresponds  to  the  point 
at  which 

(,3FS/3X)X-  =  (<^Fl/3X)  Xq  (2.3) 

where  F^  and  F.  are  the  molar  free  energies  of  the  solid  and  the 
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liquid  respectively.  Therefore,  if  nucleation  occurs  in  the 
liquid-solid  range  the  composition  of  the  solid  X'  differs  from 
that  of  the  liquid  Xq  by  a  greater  amount  than  does  the 
composition  of  the  equilibrium  solid,  XgE.  However,  the  chemical 
term  of  the  surface  free  energy  tends  to  drive  X'  in  the  opposite 
direction.  If  nucleation  occurs  below  the  solidus  temperature  it 
is  not  possible  to  generalize  if  the  first  solid  to  form  will  be 
richer  or  poorer  than  XgE> 

For  alloys  in  the  two-phase  range  of  a  eutectic  binary,  a 
critical  undercooling  must  be  achieved  for  the  nucleation  of  each 
phase,  and  the  second  phase  must  either  nucleate  heterogeneously 
on  the  primary  phase  or  homogeneously  in  the  liquid.  The  general 
conclusion  from  Christianas  analysis  is  that  the  proportion  of 
the  two  phases  in  a  solidified  alloy  within  the  eutectic  range 
will  only  agree  with  that  calculated  from  the  equilibrium  phase 
diagram  if  one  phase  acts  as  a  perfect  catalyst  for  the 
nucleation  of  the  other.  The  discussion  assumes  that  no 
recalescence  takes  place  during  solidification,  and  that  no 
foreign  catalysts  are  present. 

Several  authors  have  utilized  variations  of  the  droplet 
technique  to  study  the  catalyzing  effect  of  the  primary  phase  on 
the  nucleation  of  the  second  phase  in  several  systems  [50,53,54]. 
For  example,  undercoolings  of  68  and  176  K  (with  respect  to  the 
solidus  and  liquidus  corresponding  to  the  composition  of  the 
liquid  phase  at  the  moment  of  nucleation)  were  obtained  for  the 
nucleation  of  Si  on  Al.  Furthe  ,  e,  it  was  found  that  for 
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certain  systems  like  Al-Ge  the  undercoolings  observed  were 
greater  than  those  reported  previously  for  "homogeneous" 
nucleation,  implying  that  the  crystalline  substrate  (Al)  had  no 
influence  on  the  nucleation  of  the  second  phase  from  the  liquid. 
These  results  were  examined  by  Cantor  and  Doherty  [55]  who 
proposed  a  model  to  quantify  the  catalyzing  activity  of  the 
primary  phase,  and  concluded  that  in  some  cases  like  Zn-Bi  the 
solute  (Bi)  had  nucleated  homogeneously  from  the  liquid. 

2.2.2  Growth  kinetics 

Most  of  the  rapid  solidification  processes  are  believed  to 
take  place  in  the  continuous  growth  regime  with  a  diffuse 
interface  [56].  The  kinetic  relationship  suggested  by  Turnbull 
[57]  for  the  steady-state  continuous  growth  of  a  planar  front  is 

R=  (dlm/<5)  exp  ( -ED AT*/T* )  [1-exp  (- AHmAt*/R'TmT*)  ]  (2.4) 

where  R  is  the  interface  velocity,  Dru  is  the  self-diffusion 

IjM 

coefficient  in  the  liquid  at  the  melting  temperature,  ED  is  the 

It  it 

activation  energy  for  diffusion,  AT  and  T  are  the  interface 
undercooling  and  temperature,  respectively,  AH^  is  the  heat  of 
fusion  at  the  melting  temperature  TM,  <5  is  the  length  of  an 
interface  step  and  R'  is  the  gas  constant.  For  moderate 
undercoolings  the  above  equation  reduces  to 

R  =  e(DLM/6)  <ahm/r't2)  at* 


(2.5) 
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6  =  (  Vy)  2  (u*/uL)  (2.6) 

where  H  is  a  factor  inco^r  rated  by  Cahn  et  al.  [56]  including 

★ 

the  ratio  of  the  jumpin  • aquencies  near  the  interface  u  and  in 
the  liquid  u^,  the  average  jumping  distance  for  an  atom  to  join 
the  solid  y  and  the  length  of  an  interface  step  <5.  Values  for  £ 
vary  from  1  to  100  for  symmetric  molecules. 

Equations  (2.4)  and  (2.5)  will  be  henceforth  called  the 
exponential  and  linear  laws  respectively.  For  the  purposes  of 
this  work  the  factor  0  was  also  included  in  Equation  (2.4). 

Shaefer  and  Glicksman  developed  a  mathematical  model  to 
study  the  effects  of  curvature  and  heat  flow  in  the  transient 
growth  of  a  spherical  crystal  nucleus  in  an  infinite  undercooled 
melt  [58] .  Their  analysis  shows  that  the  initial  growth  rate  is 
very  small  due  to  the  large  effect  of  the  surface  tension  right 
after  nucleation.  The  interface  velocity  increases  rapidly  as 
the  curvature  is  reduced,  goes  through  a  maximum  and  then  reaches 
a  steady  state  whose  characteristics  depend  on  the  melt 
undercooling  and  a  kinetic  parameter. 

The  problem  of  stability  of  a  moving  interface  has  been 
extensively  discussed  in  the  literature.  For  most  cases  of  low 
interface  velocities  the  principle  of  constitutional  supercooling 
has  been  widely  used  as  an  adequate  guideline  to  determine  the 
growth  conditions  that  result  in  unstable  interfaces  during  alloy 
solidification  [59].  This  principle  is  basically  an  application 
of  thermodynamics  to  determine  if  the  liquid  ahead  of  a  moving 
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interface  is  undercooled  with  respect  to  its  composition.  The 
stability  parameter  emerging  from  this  theory  is  the  ratio  of  the 
temperature  gradient  in  the  liquid  ahead  of  the  moving  interface 
Gl  to  its  velocity  R.  Greater  stability  is  promoted  by 
increasing  G^/r. 

The  theory  of  morphological  stability  [60,61]  approaches  the 
problem  from  a  kinetic  point  of  view,  solving  the  thermal  and 
diffusion  field  equations  about  the  interface  under  steady  state 
conditions.  It  is  also  assumed  that  local  equilibrium  prevails 
at  the  solidification  front.  The  resulting  criterion  identifies 
three  factors  contributing  to  the  overall  stability  of  the 
interface.  The  solute  field  created  by  segregation  is  always 
destabilizing,  but  the  capillary  forces  always  oppose  the  growth 
of  a  perturbation.  The  thermal  field,  as  determined  by  the 
average  gradient 


GSL  =  2  (kSGS  +  kLGL*^kS  +  kL^ 


(2.7) 


will  be  stabilizing  if  GgL  is  positive,  and  destabilizing  in  the 
opposite  case.  kg  and  kL  are  the  thermal  conductivities  of  the 
solid  and  liquid  respectively.  Note  that  while  the 
constitutional  supercooling  criterion  predicts  instability  for 
all  cases  where  GL  <  0,  the  morphological  theory  may  predict 
stability  if  Gg  is  sufficiently  positive  to  make  the  average 
gradient  term  also  positive. 
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Recent  studies  have  compared  the  predictions  of  both 
theories  for  the  case  of  rapid  solidification.  Cohen  et  al. 

[62]  mapped  the  regimes  of  stability  in  a  GL  vs.  R  plot  and 
indicated  that  sufficiently  increasing  the  cooling  rate  and 
decreasing  composition  will  promote  planar  interfaces.  Other 
studies  [63,64]  have  extended  the  morphological  stability  theory 
and  indicated  that  two  effects  are  important  in  rapid 
solidification  phenomena.  The  first  one  is  known  as  absolute 
stability  and  arises  because  only  short  wavelength  perturbations 
are  important  at  high  velocities,  and  these  are  stabilized  by 
surface  energy.  The  net  result  is  that  given  a  sufficiently  high 
velocity  the  interface  will  be  stable  for  any  value  of  GL 
(assuming  that  GgL  is  positive) .  The  second  effect  is  caused  by 
departure  from  local  equilibrium  at  the  interface.  Coriell  and 
Sekerka  [64]  concluded  that  most  departures  from  local 
equilibrium  enhance  the  interface  stability,  especially  in  the 
case  where  the  partition  coefficient  tends  to  1.  Baker  [o5] 
analyzed  the  behavior  of  the  partition  coefficient  as  a  function 
of  interface  velocity  for  various  models  and  indicated  that  given 
a  sufficiently  high  velocity  the  partition  coefficient  will 
always  tend  to  unity. 

Interface  instability  usually  leads  to  cellular  or  dendritic 
patterns  depending  on  the  growth  conditions.  Analysis  of 
dendritic  solidification  has  been  approached  by  several  authors 
[66-68] .  In  general,  the  growth  of  a  dendrite  tip  into  a 
supercooled  melt  can  be  described  by  an  expression  of  the  type 
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6*  =  -Pe  ePe  Ei(Pe)  (2.8) 

9*  =  Cl(T*  -  Tm)/AHm  =  -CL  AT*/AHm  (2.9) 


Pe  =  (R  rfc/2aL) 


(2.10) 


4r 

where  0  is  the  dimensionless  tip  temperature,  Pe  is  the  Peclet 
number,  relating  the  tip  velocity  R  and  radius  rfc,  t*  and  AT*  are 
the  temperature  and  undercooling  of  the  dendrite  tip 
respectively,  and  aL  and  are  the  thermal  diffusivity  and 
specific  heat  of  the  liquid  respectively. 

Glicksman  and  Huang  [69]  have  developed  a  series  of 
relationships  between  the  microstructural  size  scale  (segregate 
spacing)  and  interface  undercooling,  and  concluded  that  for  small 
AT  the  spacing  should  be  a  fixed  multiple  of  the  critical 
nucleation  radius.  Furthermore,  they  suggested  that  as 
supercooling  increases  towards  unity  the  scale  of  the  morphology 
should  become  rapidly  smaller.  They  also  studied  the  coarsening 
process  during  solidification  and  indicated  that  the  initial 
microstructure  is  produced  far  from  thermodynamic  equilibrium  in 


order  to  comply  with  the  transport  processes  during  growth.  The 
final  microstructure  tends  to  be  substantially  different  from 
that  produced  at  the  dendrite  tip  conditions,  and  can  be 
correlated  to  the  coarsening  time  by  a  scaling  law. 
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2.2.3  Extension  of  solid  solubility 

Massalski  and  Bienvenu  have  discussed  the  limitations 
imposed  by  the  Hume-Rothery  rules  on  the  metastable  extension  in 
solid  solubility  in  the  systems  formed  by  the  quartet  Cu,  Zn,  Ag 
and  Cd  [70] .  They  indicate  that  complete  miscibility  through  the 
phase  diagram  requires  a  certain  undercooling  for  which  the  free 
energy  of  the  liquid  FL  is  higher  than  that  of  the  solid  Fg  for 
all  compositions.  If  the  size  factor  or  a  difference  in  crystal 
structure  produces  a  "hump"  or  a  discontinuity  on  the  Fg  curve, 
complete  miscibility  may  never  be  achieved. 

The  assumption  of  local  equilibrium  at  the  interface  implies 
that  extension  in  solid  solubility  results  from  suppressing 
nucleation  of  the  second  phase  below  the  eutectic  temperature. 

The  system  then  follows  the  metastable  solidus  line  as  it  moves 
into  the  two  phase  region  of  the  phase  diagram.  Therefore,  the 
maximum  solid  solubility  is  limited  to  the  composition  dictated 
by  this  metastable  solidus  at  the  solidification  temperature. 
Experiments  in  Zn-Cd,  a  system  with  a  retrograde  solidus,  have 
shown  this  hypothesis  to  be  incorrect,  and  furthermore  have 
proven  the  occurrence  of  solute  trapping  and  deviation  from  local 
equilibrium  at  the  liquid-solid  interface  [71] . 

Baker  and  Cahn  [72]  have  discussed  the  thermodynamics  of 


nonequilibrium  solidification  in  a  binary  alloy.  With  the  aid  of 
free  energy  (F)  vs.  composition  (X)  curves,  they  defined  the 
domain  of  possible  interface  compositions  for  isothermal  and 
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steady  state  solidification.  They  concluded  that  for  any  liquid 
composition  XL,  there  is  a  range  of  solid  compositions  Xs  which 
are  thermodynamically  viable.  The  maximum  Xg  that  can  be  formed 
at  the  temperature  of  interest  is  determined  by  the  point  at 
which  the  liquid  and  solid  free  energy  curves  intersect.  This  is 
known  as  the  X(Tq)  composition  and  the  locus  of  all  these  points 
forms  a  curve  between  the  liquidus  and  the  solidus  of  the 
equilibrium  phase  diagram.  Therefore,  the  Tq  curve  effectively 
limits  the  achievable  solid  solubility  at  a  given  temperature, 
regardless  of  solidification  rate.  Partitionless  transformations 
can  occur  only  if  XL  <  X(T0)  (or  equivalently  if  T  <  TQ  for  any 
XL) .  This  involves  an  increase  in  the  chemical  potential  of  the 
solute  on  crossing  the  interface  and  therefore  implies  that  the 
two  species  have  to  solidify  in  a  "cooperative"  manner.  The 
solute  is  then  "trapped"  by  a  reaction  which  produces  an  overall 
decrease  in  the  free  energy  of  the  system. 

Cahn  et  al.  [63]  have  indicated  that  for  a  system  with  a 
retrograde  solidus  the  maximum  value  of  X(Tq)  occurs  at  0  K  and 
corresponds  approximately  to  the  value  of  the  liquid  composition 
in  equilibrium  at  the  retrograde  temperature.  They  also  compared 
calculated  maxima  for  Zn-Cd  and  some  Si-base  alloys  with  the 
maximum  solubility  extension  observed  in  those  systems  during 
rapid  solidification,  and  concluded  that  the  estimated  bounds 
were  consistent  with  the  experimental  results. 

Hillert  and  Sundman  have  used  the  solute  drag  treatment  to 
determine  the  limiting  conditions  for  partitionless 
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solidification  [73].  Their  results  indicate  that  the  critical 
temperature  for  this  process  is  lower  than  the  Tq  value  and 
depends  on  the  interface  velocity  and  the  ratio  of  the 
dif fusivities  in  the  solid  and  the  liquid. 

When  nucleation  is  suppressed  below  the  solidus  temperature, 
the  first  solid  to  form  is  of  composition  very  similar  to  the 
liquid.  However,  since  heat  extraction  in  metal  droplets  is 
usually  controlled  at  the  surface,  and  the  growth  rates  are  large 
for  undercoolings  higher  than  a  few  degrees,  the  solidification 
process  after  nucleation  is  considered  essentially  adiabatic. 
Therefore,  the  extent  of  homogeneous  or  partitionless 
solidification  is  limited  by  the  reheating  of  the  metal  back  into 
the  liquid-solid  region  due  to  the  transformation  of  latent  into 
sensible  heat.  As  in  other  metastable  effects,  the  onset  of 
segregation  must  be  characterized  by  a  critical  temperature 
Tc,  whose  value  falls  presumably  between  the  solidus  and  the 
Tq  curve.  The  production  of  a  homogeneous  solid  under  adiabatic 
conditions  requires  then  an  initial  undercooling  dictated  by  the 
following  relationship 


Tc 

AHm  <  /  CL  dT  (2.11) 

TN 

where  TN  is  the  temperature  of  nucleation.  The  maximum  value  of 
Tn  that  satisfies  the  above  relationship  is  the  "hypercooling” 
temperature  Th« 
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The  heat  transfer  to  the  surroundings  is  then  more  important 
prior  to  than  during  solidification.  Hirth  [74]  has  estimated 
the  cooling  rates  necessary  for  a  particle  to  satisfy 
Equation  (2.11)  before  the  number  of  nuclei  given  by 
Equation  (2.2)  becomes  unity.  His  calculations  show  that 
homogeneous  solidification  appears  feasible  for  alloys  of  Fe  and 
Co  but  not  for  those  of  Ni,  Cu  or  Al,  under  the  cooling  rates 
currently  achievable  in  atomization  processes.  However,  these 
conclusions  are  questionable  since  he  used  the  steady  state 
expression  for  J  and  considered  homogeneous  nucle^tion  to  occur 
at  the  maximum  undercoolings  reported  in  the  emulsified  droplet 
exper iments . 


mm 


26 


CHAPTER  3 

HEAT  FLOW  ANALYSIS:  OVERVIEW 


Previous  investigations  have  indicated  that  the  production 
of  metastable  solid  structures  from  the  liquid  depends  on  the 
achievable  interface  velocity  (or  undercooling)  during  the 
solidification  process.  It  is  also  agreed  that  for  most  metals 
the  flow  of  latent  heat  away  from  the  moving  interface  is  the 
kinetically  rate  limiting  process,  or  at  least  a  major  component 
of  it. 

During  solidification  of  metal  powders,  heat  is  extracted 
from  the  droplets  by  both  convection  and  radiation  at  their 
surface.  Most  atomization  techniques  are  characterized  by  Biot 
numbers  lower  than  0.1,  where 

Bi  =  hr0/kL  O-U 

is  the  ratio  of  the  surface  conductance,  or  heat  transfer 
coefficient  h,  to  the  metal  conductance  ;  kL  is  the 

thermal  conductivity  of  the  liquid  metal  and  Kq  is  the  droplet 
radius.  It  is  then  concluded  that  the  dominant  resistance  to 
heat  flow  occurs  at  the  surface,  and  that  no  significant 
temperature  gradients  can  develop  inside  the  droplet.  This 
assumption  is  the  basis  of  the  Newtonian  model,  which  has  been 
used  in  previous  studies  on  this  topic  for  single-phase  cooling 


and  solidification  without  undercooling  [75,76]. 

The  Newtonian  approximation  will  become  less  valid  as  the 
droplet  diameter  or  the  heat  transfer  coefficient  increase,  and 
some  numerical  simulations  have  been  made  in  the  past  to  study 
this  effect  [77] .  Furthermore,  its  basic  assumption  of  uniform 
temperature  distribution  will  also  lose  validity  if  any 
substantial  undercooling  prior  to  nucleation  is  achieved  in  tha 
powders,  since  the  local  recalescence  of  the  interface  is  likely 
to  be  faster  than  the  overall  recalescence  of  the  droplet. 

In  dealing  with  rapid  solidification  processes  it  is 
commonly  accepted  that  two  extreme  heat  flow  situations  may  arise 
[69] .  In  the  first  one  the  solidification  front  moves  into  a 
superheated  liquid,  and  the  heat  of  fusion  must  be  withdrawn 
through  the  solid.  Such  are  the  cases  of  surface  melting  and 
splat  cooling,  where  the  achievable  interface  velocities  are 
limited  by  the  rate  of  heat  extraction  to  the  substrate  and  the 
thickness  of  material  being  processed.  The  second  situation 
involves  solidification  into  a  supercooled  melt,  where  most  of 
the  heat  of  fusion  is  absorbed  by  the  liquid  itself,  thereby 
raising  its  temperature  (recalescing) .  The  external  role  of  heat 
extraction  and  the  volume  of  the  casting  become  less  relevant 
during  growth,  as  the  interface  velocities  are  now  mainly 
controlled  by  the  melt  supercooling. 

Since  atomization  is  more  limited  than  surface  melting  or 
splat  cooling  in  the  achievable  rates  of  heat  extraction,  it  is 
desirable  to  reach  substantial  undercoolings  prior  to  nucleation 
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in  order  to  enhance  the  rapid  solidification  effects.  For 
example,  it  is  considered  that  homogeneous  solidification  can 
only  be  obtained  if  an  amount  of  enthalpy  equivalent  to  the  heat 
of  fusion  is  withdrawn  from  the  liquid  before  solidification 
starts  [74] .  Experimental  evidence  from  this  and  other 
investigations  [22,49,50]  has  shown  that  large  metastable  effects 
can  indeed  be  obtained  under  slow  cooling  conditions  by 
suppressing  the  onset  of  nucleation. 

It  has  also  been  contended  that  the  external  cooling  during 
atomization  may  be  sufficiently  high  to  partially  offset  the 
recalescence  of  the  droplet  and  extend  the  regime  of  rapid 
solidification.  Proper  assessment  of  this  hypothesis  requires  a 
clear  knowledge  of  the  important  process  variables  and  their 
effects  on  the  thermal  history  of  freezing  undercooled  droplets. 
Unfortunately,  previous  studies  relating  heat  flow  and  interface 
kinetics  have  been  mainly  focused  on  the  growth  of  dendrite  tips 
or  other  curved  interfaces  into  large  supercooled  melts 
[58,66-69]  and  are  difficult  to  relate  to  the  transient  heat  flow 
conditions  inside  a  powder. 

It  is  then  the  purpose  of  this  investigation  to  further 
develop  our  understanding  of  the  thermal  history  during 
solidification  of  liquid  metal  droplets.  Most  of  the  analysis 
will  be  presented  in  dimensionless  terms,  since  it  tries  to 
encompass  as  wide  a  range  of  conditions  as  possible. 

Furthermore,  general  relationships  will  be  developed  between  the 
process  variables  and  the  main  solidification  parameters. 
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3.1  Parameters 

As  will  be  evident  later  in  this  work,  the  most  important 
process  variable  is  the  droplet  size;  values  from  0.1  to  1000  ym 
were  considered  to  cover  the  typical  range  of  most  current 
laboratory  and  commercial  atomization  techniques. 

The  second  important  parameter  is  the  combined  radiative  and 
convective  heat  transfer  coefficient  at  the  droplet  surface,  h. 
Values  for  h  during  atomization  have  been  for  sometime  a  matter 
of  uncertainty,  since  its  experimental  determination  from  direct 
measurements  of  the  heat  flow  in  a  moving  droplet  will  be 
extremely  difficult,  if  not  impossible.  Nevertheless,  it  is 
possible  to  estimate  h  from  some  analytical  expressions  and  other 
indirect  methods. 

In  gas  atomization  the  convective  heat  transfer  coefficient 
h  is  overriding  and  can  be  related  to  other  process  parameters 

V 

by  the  following  expression  [78] 

Nu  »  2  +  0.6  ?¥r  (3.2) 

where  Nu,  Re  and  Pr  are  the  Nusselt,  Reynolds  and  Prandtl 
dimensionless  numbers  respectively,  given  by 

Nu  =  (2hrQ)/kG  (3.3) 

Re  =  <2rovJ/vG  <3*4) 


t 
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Pr  =  VG/UG  (  3  •  5 ) 

kG,  <1q  and  vG  are  the  thermal  conductivity,  thermal  diffusivity 
and  kinematic  viscosity  of  the  gas,  respectively,  and  V  is  the 
particle  velocity  relative  to  the  gas.  Rearranging 
Equation  (3.2)  we  obtain  an  expression  for  the  heat  transfer 
coefficient  as  a  function  of  particle  size  and  velocity 

hc  =  (kGAo)  +  0.42  [(kG/  ^  ^)  ]  (3.6) 

The  first  term  in  this  expression  represents  the  contribution 
from  conductive  heat  transfer  when  the  droplet  is  cooling  in  a 
stagnant  gas.  Values  determined  from  Equation  (3.6)  for  helium 
atomization  are  shown  in  Table  2,  where  is  also  expressed  as 
Mach  number,  or  fractions  of  the  velocity  of  sound  (965  m  s-*- 
in  He)  . 


Table  2 


Calculated  Heat 

Transfer 

Coefficients 

for  Helium 

Atomization 

Mach 

Vco 

Droplet  Diameter  (pm) 

(m  s-1 ) 

1 

10 

100 

1000 

0 

0 

5.1x10;? 

5.1x104 

5.1x10;* 

5.1x10:? 

0.1 

97 

5.7x10;? 

7.1x10’ 

1.1x10’ 

2.5x10;* 

0.3 

290 

6.2x10;? 

8.6x10’ 

1.6x10’ 

4 . 0x10  ^ 

0.  5 

483 

6.5x10^ 

9.6x10’ 

1.9x10’ 

5 . 0x10  ^ 

1.0 

965 

7 . 1x10  5 

l.lxlO5 

2. 5xl04 

6 . 9xl03 

+  from  Equation  (3.6),  in  WnT^K  ^ 

droplet  is  Aluminum  at  TM  =  933  K,  He  at  300  K 
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Analysis  of  the  hc  values  in  this  table  reveals  that  the 
cooling  efficiency  is  enhanced  by  reducing  the  particle  size 
and/or  increasing  its  velocity  relative  to  the  gas.  However,  the 
contribution  of  the  latter  decreases  with  increasing  velocity  and 
is  more  pronounced  in  the  larger  particles.  An  additional 
conclusion  from  this  table  is  that  as  the  particle  radius  tends 
to  zero  the  product  hrgf  and  hence  the  Biot  number,  becomes  a 
constant,  since  the  Nusselt  number  is  dominated  by  the  first  term 
of  Equation  (3.2).  It  should  also  be  noted  that  this  equation  was 
developed  for  laminar  flow  past  a  sphere  and  its  accuracy  when 
applied  to  gas  atomization  has  not  been  tested. 

Indirect  estimates  of  heat  transfer  coefficients  in  various 
atomization  processes  have  also  been  made  by  comparison  of 
measured  segregate  (dendrite  arm)  spacings  in  crystalline  alloy 
powders  with  predetermined  relationships  between  these  spacings 
and  average  cooling  rates  during  solidification.  Equation  (2.1). 
The  heat  transfer  coefficients  for  gas  atomization  deduced  in 
this  way  are  of  the  same  order  of  magnitude  as  those  calculated 
from  mathematical  expressions  [79] . 

In  vacuum  (e.g.  Electrohydrodynamic)  atomization  the 
radiative  heat  transfer  coefficient  hr,  is  dominant  and  can  be 
estimated  from  the  following  relationship 

hr  =  e  0  Tg/(TS  -  TG)  (3.7) 


f 
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where  e  is  the  surface  emissivity,  o  is  the  Stefan  Boltzmann 
constant,  Tg  and  TG  are  the  droplet  and  cooling  environment 
temperatures  respectively.  The  cause  for  uncertainty  in 
Equation  (3.7)  is  the  emissivity  factor,  since  values  for  liquids 
are  usually  unknown  and  values  from  0.05  to  0.3  are  reported  for 
solid  aluminum  surfaces  in  the  literature  [80] .  An  hr  of  about 
80e  W  m-^  K-^  can  be  calculated  for  for  Tg  =  933  k  and 
Tg  =  300  K. 

In  general  then,  the  combined  radiative  and  convective  heat 

—  2  —I 

transfer  coefficient  may  vary  between  10  W  m  K  for  the  pure 

radiation  case  to  less  than  10°  W  m  ^  K  for  convection  under 

the  most  favorable  experimental  conditions  [81]  .  The  limitations 

on  h  translate  into  Biot  numbers  lower  than  1  for  most  metals  and 

particle  sizes  of  interest.  For  example,  500  um  droplets  of  Al, 

Fe  and  Ni  result  in  Biot  numbers  of  0.28,  0.62  and  0.64  for  a 

S  —2  -1 

high  heat  transfer  coefficient  of  10  W  m  K  . 

The  third  major  parameter  is  the  kinetic  relationship 
between  the  interface  undercooling  and  velocity.  Although  models 
for  the  different  growth  mechanisms  have  been  fairly  well 
discussed  in  the  literature  [56],  their  experimental  verification 
has  been  limited  to  non-metallic  substances  [82,83]  and  a  few  low 
melting  metals  like  Ga  [84] . 

The  application  of  the  kinetic  models  discussed  in  section 
2.2.2  involves  an  assumption  of  the  factor  g,  whose  values  may 
range  over  two  orders  of  magnitude.  Furthermore,  data  for  the 
self-diffusion  coefficient  in  supercooled  liquid  metals  DTM 

LjM 


are 
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generally  unknown,  and  even  values  for  the  superheated  liquid  are 
unavailable  for  Al.  For  the  purposes  of  this  work  the  following 
expression  will  be  used 

Dl  =  3.3  x  10~9  exp ( -2 . 3AT/T)  m2  s'1  (3.8) 

Equation  (3.8)  was  estimated  from  the  theory  of  corresponding 
states  [85]  using  available  data  on  liquid  viscosity  above  the 
melting  point  [86] . 

The  exponential  and  linear  laws,  Equations  (2.4)  and  (2.5), 
have  been  plotted  as  R/8  versus  dimensionless  undercooling  for  Al 
in  Figure  1,  where  it  is  evident  that  the  latter  always  dictates 
higher  velocities.  The  differences  between  the  two  models  become 
more  pronounced  for  larger  undercoolings,  and  for  0*  <  -0.65 

fc  jc 

(AT  >  240  K)  they  predict  opposite  trends  with  increasing  AT  . 

Despite  the  relative  uncertainty  in  the  variables  discussed 
above,  it  will  be  shown  in  the  course  of  this  work  that  the 
conclusions  regarding  the  general  trends  and  limitations  to  the 
effect  of  heat  flow  on  the  solidification  process  are  fairly 
independent  of  the  numerical  values  of  the  individual  parameters. 

3.2  The  Enthalpy-Temperature  Diagram 

The  thermal  history  of  powder  solidification  can  be  readily 
described  on  the  enthalpy  temperature  diagram  shown  in  Figure  2. 
This  diagram  represents  graphically  the  following  equation 
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DIMENSIONLESS  INTERFACE  TEMPERATURE 

e*  =  cl(t*-tm)/ahm 


Figure  1. 


Interface  velocity  as  a  function  of  temperature  for 
the  exponential  kinetic  model  and  its  linear 
approximation. 


RELATIVE  ENTHALPY  (H-HCJ,  kJ  mole 
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(H~HSM>  =  ^AhM  +  CL(T-TMn  (1_g)  +  CS(T-TM)g  (3‘9) 

where  H  is  the  enthalpy  of  the  system  at  temperature  T  and 
fraction  solid  g,  above  the  enthalpy  of  the  solid  at  the  melting 
temperature  HgM,  AHM  is  the  heat  of  fusion  at  the  melting 
temperature  TM,  and  Cg  are  the  specific  heats  of  the  liquid 
and  solid  respectively.  Equation  (3.9)  can  be  written  in  terms 
of  the  dimensionless  enthalpy  >p  and  temperature  0  as 

ip  =  (1+0)  (1-g)  +  (Cg/cL)0g  (3.10) 

where 


=  (H"HSM>/AHM 

(3.11) 

=  cl(t-tm)/ahm 

(3.12) 

The  two  terms  on  the  right  hand  side  of  Equation  (3.10)  represent 
the  enthalpies  of  the  liquid  and  solid  fractions  respectively. 

For  a  solid  volume  (g  =  1)  the  dimensionless  nodal  enthalpy  y  is 
always  negative  or  zero  since  0  £  0,  whereas  for  a  volume  of 
superheated  liquid  (g  =  0,  B>  0)  4'  is  larger  than  unity.  Since 
constant  fraction  solid  lines  in  Figure  2  intersect  the  melting 
temperature  axis  (0  =  0)  at  4’  =  1-g,  a  volume  of  partially  solid 
material  at  TM  will  have  an  enthalpy  equal  to  its  fraction  of 
liquid.  Note  also  that  the  lines  of  constant  g  change  gradually 
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in  slope  (do/dO)  from  1  at  g  =  0  to  Cg/cL  at  g  =  1  (solid)  . 

The  system  can  conceivably  follow  any  path  of  decreasing 
enthalpy  in  the  supercooled  region  (41  <  1)  starting  at  the 
nucleation  temperature  (eN  <  0)  on  the  liquid  line,  and  ending  at 
some  point  on  the  solid  line.  A  vertical  path  represents 
isothermal  solidification,  of  which  a  limiting  case  is  the 
no-undercooling  situation  (path  1  in  Figure  2).  A  horizontal 
path  represents  isenthalpic  (or  adiabatic)  solidification.  If 
nucleation  can  be  suppressed  until  an  amount  of  enthalpy  at  least 
equivalent  to  the  heat  of  fusion  has  been  removed  (i.e.  tp  <  0) 
the  system  can  undergo  complete  solidification  without  further 
heat  extraction  (path  2  in  Figure  2) .  The  region  of  the  diagram 
for  which  |i  <  0  is  known  as  "hypercooled" .  A  more  realistic  path 
such  as  that  described  by  curve  3  will  be  discussed  later. 

The  present  investigation  will  address  both  undercooling  and 
no-undercooling  situations,  and  numerical  models  will  be 
developed  for  each  case.  The  Newtonian  model  will  also  be 
reviewed  and  extended  to  use  it  as  a  basis  of  comparison  and 
develop  trends  for  the  limiting  case  where  the  resistance  to  heat 
flow  inside  the  droplet  is  negligible.  Although  most  of  the  work 
is  focused  on  powders  coming  from  a  single  nucleus,  the 
implications  of  multiple  nucleation,  especially  when  leading  to 
microcrystalline  structures,  will  also  be  discussed. 
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CHAPTER  4 

PROBLEM  STATEMENT  AND  SOLUTION  APPROACH 

We  consider  the  cooling  of  a  superheated  pure  liquid  metal 
droplet  and  its  subsequent  solidification  due  to  radially  uniform 
heat  loss  by  radiation  and  convection  from  its  surface.  Surface 
conductance — heat  transfer  coef f icient--and  the  temperature  of 
the  environment  are  assumed  to  remain  constant.  The  thermal 
properties  of  the  solid  and  the  liquid  are  considered  to  be 
different,  but  independent  of  temperature.  Heat  transfer  inside 
the  droplet  occurs  only  by  conduction. 

First,  a  complete  analysis  based  on  the  Newtonian  model  will 
be  used  to  develop  some  general  relationships  and  their  trends 
with  changing  particle  size,  heat  transfer  coefficient,  etc.  The 
liquid  cooling  stage  and  solidification  with  and  without 
undercooling  will  be  discussed,  both  for  one  nucleation  event  and 
massive  nucleation  cases. 

Next,  a  dimensionless  form  of  the  analytical  solution  to  the 
single-phase,  one-dimensional  heat  conduction  equation  in 
spherical  coordinates  will  be  used  to  study  the  liquid  cooling 
stage  and  its  results  compared  with  the  previous  Newtonian 
relationship  and  trends.  It  will  also  be  used  to  determine  the 
temperature  profile  in  the  droplet  when  its  surface  reaches  the 
prescribed  nucleation  temperature.  This  profile  is  then  used  as 
the  initial  boundary  condition  for  the  subsequent  solidification 
problem. 
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The  next  chapter  involves  the  numerical  analysis  of  the 
solidification  stage  with  no  undercooling  required  for  nucleation 
or  growth.  The  concepts  of  the  enthalpy  model  proposed  by 
Shamsundar  and  Sparrow  [87]  are  used  to  develop  the  solutions  for 
this  portion  of  the  problem.  In  this  model  the  enthalpy  is  used 
as  a  dependent  variable  in  addition  to  the  temperature,  and  the 
finite  difference  equations  are  written  in  their  fully  implicit 
form.  The  important  variables  during  cooling  and  solidification 
of  the  droplets  are  determined  and  expressed  in  their  most 
general  dimensionless  form  over  a  range  of  Biot  numbers  of 
practical  interest  for  droplets  of  three  elements:  aluminum, 
nickel  and  iron. 

Finally  the  numerical  modelling  of  the  solidification  of  an 
undercooled  liquid  droplet  from  a  single  nucleation  event 
occurring  at  the  surface  will  be  described.  The  problem  is  two- 
dimensional  in  nature  and  has  a  symmetry  axis  going  from  the 
nucleation  point  on  the  surface  to  its  antipode.  A  modified 
version  of  the  enthalpy  model  is  developed  and  the  mesh  is 
defined  on  a  superimposed  bispherical  (rotational  bipolar) 
coordinate  system.  Since  the  temperature  of  each  solidifying 
node,  i.e.  that  containing  a  portion  of  the  liquid-solid 
interface,  is  allowed  to  depart  from  the  melting  point,  a  kinetic 
relationship  must  be  added  to  compensate  for  this  extra  degree  of 
freedom  in  the  nodal  enthalpy  equations.  As  before,  the  finite 
difference  equations  are  written  in  their  fully  implicit  form  and 
solved  by  a  Gauss  Seidel  iteration  routine.  The  results  are 
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compared  with  the  equivalent  Newtonian  case  and  approximate 
interface  velocities  and  temperature  gradients  in  front  of  the 
liquid  solid  front  are  calculated.  These  last  parameters  are 
later  discussed  in  terms  of  the  interface  stability  theory. 

For  mathematical  consistency,  all  the  dimensionless  numbers 
in  the  following  formulation  are  written  in  terms  of  the 
properties  of  the  liquid.  The  values  of  the  thermophysical 
properties  used  in  the  calculation  are  given  in  Appendix  A. 
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CHAPTER  5 

NEWTONIAN  MODEL 

The  Newtonian  or  "point  source"  model  assumes  that  all  the 
resistance  to  heat  flow  is  concentrated  at  the  surface  and 
therefore  no  temperature  gradients  can  develop  inside  the 
droplet.  It  is  commonly  accepted  that  this  model  may  be  applied 
when  the  Biot  number  is  lower  than  0.1  [78];  however,  it  will  be 
later  shown  that  this  criterion  may  not  be  always  satisfactory  in 
powder  solidification. 

According  to  the  Newtonian  model,  the  thermal  history  of  a 
cooling  droplet  can  be  described  by 

(dii>/ dFo)  *  -3  Bi  (0  +  Ste)  (5.1) 

where  \p  and  0  are  the  dimensionless  enthalpy  and  temperature 
respectively,  Fo  is  the  Fourier  number  or  dimensionless  time, 

Fo  =  aLt/rQ  (5.2) 

and  Ste  is  the  Stefan  number, 

Ste  =  CL  (TM  "  V/^M  *5-3) 

which  is  equivalent  to  the  dimensionless  temperature  of  the 
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cooling  environment  with  opposite  sign. 

5.1  Liquid  Cooling  Stage 

For  single  phase  cooling  the  rate  of  enthalpy  change  in  the 
liquid  is  given  by 


(dt^/dFo)  =  (dO/dFo)  (5.4) 

(for  the  solid  phase  cooling  the  right  hand  side  is  multiplied  by 
the  ratio  of  solid  to  liquid  specific  heats  CS/CL) .  Combining 
Equations  (5.1)  and  (5.4)  produces  an  expression  for  the  actual 
instantaneous  cooling  rate 

e  =  3  ((VkL>  (h/r0)  <T  ~  tg>  (5.5) 

where  r^  is  the  droplet  radius,  T  and  Tq  are  the  liquid  and 
cooling  environment  temperatures  respectively.  The  sign  of  the 
temperature  derivative  has  been  changed  so  that  c  is  positive  for 
T  >  Tg. 

Equation  (5.5)  can  be  now  readily  integrated  to  yield  an 
expression  for  the  cooling  time,  i.e.  the  time  it  takes  for  the 
droplet  to  cool  from  an  initial  temperature  Tq  (Oq  >  0)  to  the 
prescribed  nucleation  temperature  (0N  <  0) 

Fon  =  ln[ ( 6 0  +  Ste)/(0N  +  Ste) ]  /  3Bi 


(5.6) 
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fcN  =  (kL/3aL)  (rQ/h)  ln{(TQ  -  TQ)/(TN  -  TQ) J  (5.7) 

The  above  equations  give  the  cooling  time  dependence  on  the 
process  parameters  in  dimensionless  and  dimensional  forms 
respectively. 

In  general,  Equations  (5.5)  to  (5.7)  indicate  that 
decreasing  the  droplet  radius  by  one  order  of  magnitude  will 
proportionally  increase  the  cooling  rate  and  decrease  the  cooling 
time.  Increasing  the  heat  transfer  coefficient  h  will  have  the 
same  effect. 

5.2  Solidification  without  Undercooling 

For  isothermal  solidification  at  the  melting  point  (6=0) 
the  rate  of  enthalpy  change  is 

- (dtp/dFo)  =  (dg/dFo)  =  3  Bi  Ste  (5.8) 

where  (dg/dFo)  is  the  dimensionless  solidification  rate.  This 
equation  can  be  readily  integrated  to  yield  the  expressions  for 
the  dimensionless  and  dimensional  solidification  times 

Fof  =  (3  Bi  Ste)-1  (5.9) 

tf  =  [  AH^ p/3  (Tj,|  -  Tq)  ]  (fg/h) 


(5.10) 
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where  TM  is  the  melting  temperature,  AHM  and  p  are  the  heat  of 
fusion  and  density  of  the  material  respectively. 

Equations  (5.8)  to  (5.10)  indicate  that  the  solidification 
rate  and  time  depend  on  r^  and  h  in  the  same  way  as  the  cooling 
rate  and  time  discussed  above. 

The  Newtonian  model  per  se  does  not  imply  any  particular 
geometry  for  the  liquid  solid  interface.  Nevertheless,  one  could 
conceivably  calculate  interface  velocities  from  this  model  by 
prescribing  such  a  geometry.  Manipulation  of  Equation  (5.8) 
shows  that,  in  general,  the  dimensionless  (R)  and  dimensional  (R) 
interface  velocities  can  be  written  as 

5  =  (ro/°lL)  R  =  Bi  Ste  g'  (5.11) 

R  =  I<Tm  -  TG)/pAHM]  h  g'  (5.12) 

where  g'  is  a  geometric  factor  solely  dependent  on  the  fraction 
solid.  Therefore,  given  a  fixed  growth  geometry,  the  Newtonian 
interface  velocity  at  a  certain  fraction  solid  is  insensitive  to 
changes  in  the  droplet  radius  and  can  only  be  increased  by 
increasing  the  heat  transfer  coefficient  or  lowering  the  cooling 
environment  temperature. 

It  is  evident  from  Equation  (5.11)  that  the  calculated 
velocities  will  be  strongly  dependent  on  the  shape  asigned  to  the 
interface.  For  example,  if  the  interface  is  concentric  with  the 
droplet  surface  (concentric  geometry  in  Figure  3)  the 
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CONCENTRIC  GEOMETRY 
INTERFACE  MOVES 
TOWARD  CENTER  4>  =  0 


AXISYMMETRIC  GEOMETRY 
INTERFACE  MOVES  FROM 
Z*  =  -1  TO  Z*  =  1 


Interface  geometries  used  for  the  cases  of 
no-undercooled  (concentric)  and  undercooled 
(axi symmetric)  solidification. 


dimensionless  interface  velocity  R  is  given  by 


R  =  (d4>*/dFo)  =  Bi  Ste  {(4>*)-2} 


(5.13) 


where  <f>*  is  the  dimensionless  radial  position  of  the  liquid  solid 

it 

front  (0  <  4>  <  1)  .  However, if  the  interface  is  spherical  in 

shape,  with  its  center  of  curvature  on  the  droplet  surface 
(axisymmetr ic  geometry  in  Figure  3)  the  dimensionless  interface 
velocity  is  now 


R  =  (dZ*/dFo)  =  Bi  Ste  { 4/ [ 2 (1+Z*) 2- (1+Z*) 3] }  (5.14) 

it 

where  Z  is  the  dimensionless  position  of  the  solid  front  along 
the  axis  of  symmetry  for  both  spheres,  henceforth  called  the 
growth  axis  (-1  <  Z*  <  1) . 

Normalized  forms  of  Equations  (5.13)  and  (5.14)  are  plotted 
in  Figure  4,  showing  that  quite  different  results  may  be  obtained 
simply  by  assuming  two  different  growth  geometries.  For  example, 
at  a  fraction  solid  g  =  0.1  the  velocity  predicted  from  the 
concentric  pattern  is  3.5  cm  s~^  whereas  that  obtained  from  the 
axisymmetr ic  geometry  is  23.3  cm  s-^  for  A1  droplets  of  any  size 
(h  =  5  x  10^  W  m  2  K  =  300  K)  .  It  should  also  be  noted 

that  this  figure  is  independent  of  the  metal  properties. 


47 


Figure  4.  Normalized  interface  velocity  as  a  function  of 

fraction  solid  for  Newtonian  solidification  without 
undercooling.  Note  that  the  values  in  the  ordinate 
are  equivalent  to  g"  in  Equation  (5.11)  and  therefore 
independent  of  the  material. 
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5.3  Solidification  with  Undercooling 

When  a  liquid  droplet  is  substantially  undercooled  prior  to 
nucleation,  the  solidification  process  cannot  be  assumed 
isothermal  and  the  system  acquires  one  more  degree  of  freedom. 
Although  the  thermal  history  must  still  follow  Equation  (5.1), 
the  problem  cannot  be  easily  solved  analytically  and  a  numerical 
approach  is  required. 

5.3.1  Mathematical  formulation 

Equation  (5.1)  can  be  written  in  implicit  finite  difference 
form  as  follows 

^k+1  +  ^BiAFo  0|<+i  *  -  3BiAFo  Ste  (5.15) 

where  the  subscripts  k  and  k+1  indicate  the  previous  and  current 
time  steps  respectively,  separated  by  an  interval  AFo.  From  the 
definition  of  enthalpy  given  by  Equation  (3.10)  we  obtain 

^k+1  "  (1+Ac9k+1)  ek+l  =  1  "  9k+i  (5.16) 

where 

AC  =  (Cg  -  CL)/CL  (5.17) 

Combining  Equations  (5.15)  and  (5.16)  produces  the  basic 
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relationship  of  the  Newtonian  model  for  undercooled 
solidification: 


(l+3BiAFo)0k+1  +  (ACek+1-l)gk+1  -  (^k-l-3BiAFoSte)  =  0  (5.18) 


In  order  to  solve  this  equation  as  a  function  of  time  we 
require  a  kinetic  relationship  between  g  and  0  to  account  for  the 
extra  degree  of  freedom  introduced  in  non-isothermal 
solidification.  Since  the  kinetic  models  relate  the  undercooling 
to  linear  displacements,  it  is  necessary  to  specify  the  number  of 
nucleation  events  and  the  liquid-solid  interfacial  geometry  to 
translate  these  displacements  into  changes  of  the  fraction  solid. 
If  we  assume  that  no  temperature  gradients  exist  in  the  droplet, 
and  neglect  any  orientation  dependence  of  the  growth  rate,  the 
interface  should  preserve  the  original  nucleus  geometry. 

In  this  section  we  will  first  develop  the  model  for  a  single 
nucleation  event  and  later  approach  the  problem  of  massive 
nucleation  that  leads  to  a  microcrystalline  structure. 

5 . 3 . 1 . 1  Single  nucleation  event 

We  now  consider  a  nucleation  event  in  the  form  of  a 
spherical  cap  occurring  at  the  droplet  surface.  The  interface 
follows  the  axisymmetric  geometry  depicted  in  Figure  3.  The 
fraction  solidified  is  given  by 
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(5.19) 


g  =  0 . 5 ( 1+Z*) 3  -  0.1875(1+Z*)4 

where  (1+Z  )  is  the  distance  solidified  along  the  growth  axis. 

For  any  kinetic  expression,  the  interface  position  is 

Zk+1  =  zk  +  B(0)  AFo  (5.20) 

where  R(0)  is  the  temperature  dependent  dimensionless  interface 
velocity.  Substituting  the  above  equation  into  (5.19)  and 
rearranging,  we  obtain  the  relationship  between  g  and  0  which  is 
of  the  form 


p=4 

gk+l  =  l  Yn  t£(e)JP  (5.21) 

p=0 

where  the  coefficients  Yp  are  only  functions  of  the  time  interval 
AFo  and  the  interface  position  in  the  previous  time  step  z£. 

5 . 3 , 1 . 2  Massive  nucleation 

For  the  simulation  of  the  microcrystalline  solidification  we 
choose  an  array  of  close  packed  dodecahedral  grains,  each  one 
originating  from  a  spherical  nucleus.  If  the  spacing  between  two 
nearest  neighbors  is  d^x,  the  approximate  number  of  grains  in  a 
droplet  of  radius  rQ  is 

N  =  (Tt/3rZ)  (2rQ/dMX)3 


(5.22) 
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Assuming  that  the  spherical  solid  grows  concentrically  within  the 
dodecahedral  liquid  volume,  the  fraction  solid  is  given  by 

g  =  (tt/^2)  (Z*/3)  0  <  Z*  <  1  (5.23) 

g  =  (tt/^2)  [-(5/3)  (Z*)  3  +  3(Z*)2  -  1]  1  <  Z*  <  2/^f  (5.24) 

where  Z*  is  now  the  distance  between  the  interface  and  the 
nucleation  point  normalized  by  half  the  grain  spacing  dMx. 

Equation  (5.23)  applies  when  the  spheres  are  growing  freely 
inside  the  dodecahedra.  Equation  (5.24)  applies  from  the  moment 
when  two  neighbor  spheres  touch  (Z*  =  1)  until  3  spheres  touch 
and  start  forming  a  corner  (Z*  =  2/^3) .  At  this  moment  the  value 
of  g  is  approximately  0.96.  Since  the  mathematical  expressions 
for  the  rest  of  the  solidification  problem  are  rather  complex, 
the  analysis  was  terminated  at  this  point. 

The  interfacial  displacements  are  given  in  this  case  by 

<+l  -  zjj  +  (2r0/dMX)  R(9)  AFo  (5.25) 

which  combined  with  Equations  (5.23)  or  (5.24)  yields  an 
expression  similar  to  (5.21) 

p=3 

“  [<2VdMX>  YP]  E£(0nP 
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(5.26) 


where  the  coefficients  Y'  are  only  functions  of  Fo  and  the  value 
of  Z*. 

Equation  (5.18)  can  be  solved  with  the  aid  of 
Equations  (^.21)  or  (5.26)  for  the  single  and  massive  nucleation 
cases  respectively,  using  a  Newton-Raphson  technique  at  each  time 
step,  starting  with  as  the  initial  condition.  Calculations 
were  performed  with  both  the  linear  and  exponential  laws  defined 
in  Chapter  2.  The  solution  procedure  was  set  for  an  adjustable 
time  increment  to  produce  roughly  uniform  displacements  of  the 
liquid-solid  interface. 

Before  entering  the  description  of  the  results,  it  should  be 
emphasized  that  the  inherent  assumption  of  Newtonian  conditions 
during  solidification  precludes  temperature  gradients  in  the 
droplet.  Therefore,  no  distinction  can  be  made  between  droplet 
and  inter  facial  temperature,  and  the  motion  of  the  solid  front 
will  be  dictated  by  the  overall  undercooling. 

5.3.2  Results  and  discussion 

The  general  shape  of  the  enthalpy-temperature  (<|>-8)  curve 
has  been  sketched  in  Figure  2  (case  3)  where  two  distinct 
solidification  regimes  can  be  identified.  The  first  one  is  the 
recalescence  or  "rapid  solidification"  stage,  where  the 
supercooled  droplet  is  absorbing  most  of  the  liberated  heat  of 
fusion  and  the  heat  loss  to  the  surroundings  is  more  or  less 
irrelevant.  The  second  regime  is  one  of  a  slower  growth,  and 

i 

I 

i 
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develops  after  the  droplet  undercooling  has  been  largely 
relieved,  therby  limiting  the  progress  of  solidification  to  the 
conditions  imposed  by  the  external  heat  flow.  This  regime 
usually  involves  temperatures  close  to  the  melting  point  of  the 
material  and  is  roughly  isothermal. 

The  thermal  history  reflects  the  competition  between  the 
external  rate  of  heat  extraction  and  the  recalescence  rate 
(release  of  latent  heat)  which  is  proportional  to  the 
solidification  rate  (change  in  fraction  solid  with  time,  dg/dt) . 
The  latter  is  in  turn  a  product  of  the  liquid  solid  interface 
velocity  times  its  area.  Predominance  of  the  recalescence  rate 
translates  into  a  shallow  negative  slope.  As  the  heat  extraction 
process  is  able  to  relieve  some  of  the  heat  of  fusion  being 
transformed  into  sensible  heat,  the  slope  of  the  curve  becomes 
steeper  (more  negative)  and  eventually  positive  if  the  heat  flow 
is  dominant. 

In  addition  to  the  two  regimes  noted  above,  there  is  usually 
an  initial  transient  in  the  ip  — 0  curve,  during  which  time  the 
cooling  rate  decreases  rapidly  from  the  liquid  value  at  the 
moment  of  nucleation  until  it  changes  sign  and  the  overall 
recalescence  begins.  At  the  start  of  solidification  the  velocity 
is  finite  but  the  interfacial  area  is  very  small,  and  so  is  the 
solidification  rate.  Therefore,  the  rate  of  heat  extraction 
dominates  the  rate  of  recalescence  and  the  droplet  temperature 
drops  further  below  0N>  There  is  also  a  final  transient,  of 
increasing  cooling  rate,  due  to  the  same  competitive  process. 


If  we  consider  that  the  enthalpy  is  only  a  function  of  the 


temperature  and  fraction  solid,  then 

(d^/d6)  =  (3^/90)g  +  ( 9^/9g) q  (dg/d9)  (5.27) 

and  from  Equation  (3.10) 

(9^/96)  =  1  +  AC  g  (5.28) 

(3v/3g)Q  =  Ac  6  -  1  (5.29) 

Equation  (5.27)  can  also  be  written  as 

(di^/de)  [l-(9^/9g)0(dg/dZ*)  (dZ*/d^)]  =  { 3^/39)  (5.30) 

dZ*  and  d^  are  given  respectively  by 

dZ*  =  R  dFo  (5.31) 

d^  =  -3  Bi  (0  +  Ste)  dFo  (5.32) 

Combining  Equations  (5.28)  to  (5.32)  we  obtain 

(dt^/d0 )  { 1+  [  ( AC0+1)  /3  (0+Ste)  ]  (dg/dZ*)  (R/Bi)  } -1-ACg  =  0  (5.33) 


Since  the  interface  velocity  R  can  be  written  in  its  general 
form  as 

R  =  K  4>(6)  (5.34) 

where  K  is  a  dimensionless  kinetic  coefficient  at  the  melting 
point,  and  4>(0)  is  a  function  of  temperature.  K  is  related  to 
the  actual  kinetic  parameter  KM  (i.e.  the  constant  relating 
velocity  and  undercooling  in  the  linear  law,  Equation  (2.5))  by 

K  =  (r0/aL>  <ahm/cL>  kM  (5.35) 

KM  =  StD^/fi)  (AHm/R'T2)  (5.36) 

Equation  (5.33)  can  now  be  written  as 


(di^/d© )  { 1+  [  ( AC0+1) /3 ( 0+Ste)  ]  (dg/dZ*)  (K/Bi)  0(0)} -1- ACg  =  0  (5.37) 

Since  both  K  and  Bi  contain  an  rQ  factor,  their  ratio  is 
independent  of  particle  radius.  In  addition,  (dg/dZ*)  is  only  a 
function  of  g  for  a  fixed  interfacial  geometry.  Therefore,  the 
ip—  0  path  described  by  the  integration  of  Equation  (5.37)  is  also 
independent  of  particle  size.  The  only  effect  of  Tq  is  to 
increase  or  decrease  proportionally  the  time  scale  of  the 
process . 
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For  a  particular  solidification  geometry  (e.g.  concentric), 
the  thermal  history  and  the  interface  velocity  at  a  given 
fraction  solid  are  only  functions  of  the  initial  temperature 
0N,  and  the  dimensionless  ratio  K/Bi.  This  latter  parameter 
measures  the  degree  of  predominance  of  the  recalescence  rate  over 
the  rate  of  heat  loss  to  the  environment.  Appropriate 
substitution  of  the  variables  in  Equation  (5.36)  yields  an 
expression  for  the  K/Bi  ratio  in  terms  of  0  and  h,  which  is  used 
throughout  the  calculations 

(K/Bi)  =  2  x  107  (iS/h)  (5.38) 

Kinetic  considerations  generally  limit  8  to  within  1  and  10*, 
while  the  heat  transfer  coefficients  during  atomization  range 

c.  _o  —1 

from  a  minimum  of  10  to  about  1GU  W  m  K  .  Consequently,  the 

o 

parameter  K/Bi  can  vary  from  about  20  to  2  x  10  . 

For  the  case  of  massive  nucleation  a  third  parameter,  the 
ratio  2rQ/dMX  (or  the  number  of  grains) ,  will  appear  in 
Equation  (5.33)  multiplying  R. 

5 . 3 . 2 , 1  Solidification  from  a  single  nucleation  event 

Figure  5  shows  the  effect  of  the  initial  undercooling  on  the 
enthalpy-temperature  curve  for  a  constant  K/Bi  =  100,  considering 
both  linear  and  exponential  kinetics.  It  is  evident  that 
increasing  the  undercooling  prior  to  nucleation  results  in  a 
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corresponding  increase  in  the  solidification  rate,  thus  reducing 
the  initial  transient  and  the  relative  contribution  of  the 
external  heat  extraction.  The  latter  is  also  diminished  due  to 
the  smaller  temperature  differences  between  the  droplet  and  its 
environment.  In  addition,  this  figure  also  shows  that  the  "slow 
growth"  regime  is  gradually  suppressed  as  r’N  decreases  (or  as 
fraction  solid  at  the  end  of  recalescence  increases  towards  one) . 

Changing  the  kinetic  model  from  linear  to  exponential 
(slower) ,  has  the  effect  of  slightly  extending  the  "rapid 
solidification"  regime,  as  shown  in  Figure  5.  This  follows  from  a 
reduction  in  the  solidification  and  recalescence  rates  due  to  the 
lower  interface  velocities  predicted  by  the  exponential  law, 
therby  enhancing  the  role  of  heat  loss  to  the  environment.  The 
differences  between  the  two  models  increase  with  undercooling,  as 
expected  from  Figure  1, 

The  dependence  of  the  thermal  history  on  the  ratio  K/Bi  for 
a  constant  initial  temperature  0 ^  =  -0.5  ( ATN  =  182  K)  is  shown 
in  Figure  6.  Curves  for  K/Bi  >  10  closely  approach  the  limiting 
case  of  adiabatic  ( isenthalpic)  rapid  solidification  followed  by 
an  isothermal  slow  growth  regime  at  the  melting  temperature.  As 
K/Bi  decreases  (slower  kinetics  or  more  efficient  cooling) ,  the 
curves  show  an  increasing  contribution  of  the  heat  extraction 
rate  during  recalescence.  At  the  extreme,  e.g.  K/Bi=10,  the 
regime  of  "rapid  solidification"  is  extended  to  a  point  where  the 
droplet  temperature  never  goes  beyond  its  initial  value  eN« 
However,  even  for  a  low  kinetic  coefficient  KM  =  2  cm  s-^  K~* 


TEMPERATURE  (K) 


700  800  900 


DIMENSIONLESS  TEMPERATURE  0  =  CL(T-T^)/AH^ 


Figure  6.  Effect  of  the  ratio  K/Bi  on  the  thermal  history  of 
undercooled  droplets  for  a  fixed  nucleation 
temperature.  Note  that  for  typical  atomization 
conditions  (K/Bi>100)  the  external  cooling  is  not  able 
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(3=1),  the  heat  transfer  coefficient  would  have  to  be  about 

2  x  10^  w  m- ^  k-^  in  order  to  produce  this  effect.  For  a  higher 

6  such  as  the  value  of  100  suggested  by  Cahn  et  al.  [56],  we 

8  —2  —1 

would  require  an  h  =  2  x  10  W  m  K  ,  which  is  clearly  above 
the  typical  values  encountered  in  atomization. 

Three  characteristic  undercooling  values  can  be 
distinguished  in  the  ip—  0  curves  (see  Figure  6) 

a)  The  maximum  undercooling  at  the  end  of  the  initial 
transient,  ATmax  (or  9max) .  This  parameter  is  of  small 
importance  since  6m3X  is  close  to  0N  for  typical 
atomization  conditions,  i.e.  K/Bi  >  100. 

b)  The  minimum  undercooling,  ATmin  (or  6min) ,  corresponding 
to  the  "nose"  that  every  ^-0  curve  exhibits  after  the 
recalescence  regime.  This  point  represents  the  maximum 
temperature  reached  by  the  droplet  during  solidification, 
where  the  external  rate  of  heat  extraction  becomes  again 
the  dominant  process,  and  changes  the  slope  of  the  curve 
from  negative  to  positive.  If  a  certain  metastable 
effect  requires  a  critical  undercooling  ATC  during 
growth,  the  "nose"  of  the  recalescence  curve  should 

always  have  a  AT  .  >  ATr,. 

J  min  t 

c)  The  undercooling  in  the  droplet  at  a  given  fraction  solid 
AT(g)  or  0(g).  In  the  limiting  case  of  adiabatic 
solidification  (K/Bi-^<»)  ,  the  undercooling  at  the  end  of 
recalescence  is  zero  and  the  fraction  solidified  is 

ga  =  ~aN.  The  effect  of  heat  extraction  during 
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recalescence  (finite  K/Bi)  can  then  be  assessed  by 
comparing  the  values  of  AT(ga) .  The  ratio  of 
AT(ga)/ATN  indicates  the  fraction  of  the  initial 
undercooling  still  present  in  the  droplet  when  g  =  ga, 
and  therefore  is  a  measure  of  how  efficiently  the 
external  cooling  is  able  to  slow  down  the  reheating 
process. 


Table  3 

Effect  of  the  K/Bi  Ratio  on  the  Characteristic 
Undercoolings  of  the  Enthalpy-Temperature  Curve  (K) ' 


Linear 

Exponential 

K/Bi 

AT(ga) 

AT  . 
min 

AT  (gd) 

AT  . 
min 

10  2 

51 

23 

65 

26 

200 

31 

11 

38 

12 

400 

18 

5.4 

22 

5.4 

1600 

6.3 

1.3 

7.3 

1.3 

10  4 

1.5 

0.21 

1.6 

0.21 

10  6 

0.004 

0.002 

0.078 

0.002 

N  -  182 

K,  h  =  5xl04  Wm'2K~1, 

KM  "  2 

cms-1K~1 

Table  3  shows  the  characteristic  AT  values  as  a  function  of 
K/Bi  for  both  the  linear  and  exponential  kinetic  models.  The 
results  indicate  that  even  for  a  very  low  K/Bi  like  100,  the  heat 
flow  can  only  produce  an  undercooling  AT(ga)  =  65  K.  This  means 
that  about  one  third  of  the  liberated  heat  of  fusion  may  be  lost 
to  the  environment  before  the  droplet  reaches  a  fraction  solid  of 
ga  =  0.5.  One  can  then  conclude  from  the  Newtonian  analysis  that 
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the  external  cooling  has  little  influence  on  the  thermal  history 
of  the  droplet  during  recalescence . 

A  further  conclusion  from  Figures  5  and  6  is  that  the  K/Bi 
ratio  is  more  influential  in  shifting  the  "nose"  of  the  ^-0  curve 
than  the  initial  undercooling  0  if  atc  is  relatively  small 
rapid  solidification  effects  may  be  obtained  throughout  the 
droplet  with  moderate  AT^,  provided  the  value  of  K/Bi  is 
sufficiently  small  to  keep  ATmin  above  ATC.  However,  if  ATC  is 
large,  the  extent  of  rapid  solidification  will  be  essentially 
limited  by  the  undercooling  prior  to  nucleation. 

Although  the  thermal  history  is  uniquely  defined  for  a  given 

0N  and  K/Bi,  (i.e.,  the  undercooling  is  fixed  for  every  fraction 
solid) ,  the  Newtonian  interface  velocities  depend  on  the  specific 
values  of  kM'  h  and  the  prescribed  kinetic  model.  This  is  due  to 
the  assumption  of  infinite  heat  conductance  inside  the  droplet, 
making  the  interface  motion  independent  of  the  thermal  field  in 
its  vicinity.  For  the  Newtonian  case,  then,  the  interface  will 
move  as  fast  as  the  overall  undercooling  of  the  droplet  allows  it 
to  go,  rather  than  adjust  its  undercooling  to  comply  with  the 
heat  flow  conditions  around  it. 

The  assumption  5f  an  axisymmetric  solidification  geometry 
influences  the  calculated  growth  rates,  as  previously  discussed 
for  the  case  of  no  undercooling.  In  order  to  isolate  the 
geometric  effects  from  those  imposed  by  the  kinetics  and  heat 
flow  conditions,  a  normalized  interface  velocity  is  defined  as 
the  ratio  of  undercooled  (Ry)  to  the  no  undercooled  (RNU) 
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interface  velocities  for  the  same  fraction  solid  g.  This  ratio 
is  independent  of  particle  size,  as  shown  before,  as  well  as  of 
the  kinetic  coefficient  KM,  making  it  a  single  valued  function  of 
g  for  constant  and  K/Bi. 

Figure  7  shows  the  effect  of  fiN  on  the  normalized  interface 
velocity  for  a  fixed  K/Bi  =  10^,  whereas  Figure  8  describes  the 
effect  of  K/Bi  for  a  constant  0N  =  -0.5.  The  curves  in  these 
figures  consist  of  a  "hump"  where  Ry  >>  rnu  reflecting  the 
undercooling  driven  velocity,  and  the  flat  portion  where 
RU  a  rNU  characteristic  of  the  heat  flow  driven  velocity. 

The  initial  and  final  transients  are  geometric  effects  since  the 
no-undercooling  model  predicts  infinite  velocities  at  the 
beginning  and  at  the  end  of  solidification  (see  Figure  4)  but  the 
kinetics  will  limit  the  interfacial  motion.  It  should  be 
emphasized  that  these  velocities  will  only  be  meaningful  if 
Newtonian  conditions  are  fulfilled.  Results  of  the  Enthalpy 
model  for  undercooled  solidification,  discussed  in  Chapter  8, 
reveal  that  substantial  temperature  differences  during 
recalescence  may  arise  in  the  droplet  even  for  Bi  numbers  as  low 
as  10”^.  Therefore,  care  should  be  exercised  in  using  the 
Newtonian  model  for  rapid  solidification  of  undercooled  metals. 
Better  guidelines  are  necessary  to  determine  its  applicability, 
since  compliance  with  the  condition  of  Bi  <  0.1  does  not  seem  to 
be  a  satisfactory  criterion. 

The  dependence  of  the  solidification  times  on  the  K/Bi  ratio 
for  a  linear  kinetic  model  and  different  nucleation  temperatures 


NORMALIZED  INTERFACE  VELOCITY  (R../R, 


|  ALUMINUM  | 
NEWTONIAN  MODEL 
LINEAR  KINETICS 
K/Bi  =  103 


=  -1.0 


0.4  0.6 

FRACTION  SOLID  g 


Figure  7.  Normalized  interface  velocity  for  undercooled 

Newtonian  solidification,  as  a  function  of  fraction 
solidified  and  nucleation  temperature.  When 
Ry/R-jU  =  1  the  interface  motion  is  solely  controlled 
bv  tne  heat  flow. 
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Figure  8.  Normalized  interface  velocity  for  undercooled 

Newtonian  solidification,  as  a  function  of  fraction 
solidified  and  (K/Bi) .  Note  that  for  K/Bi>103  the 
interface  falls  under  pure  heat  flow  control  shortly 
after  g  surpases  the  adiabatic  value  ga  =  -0M. 


NORMALIZED  SOLIDIFICATION  TIME  [(3Bi  Ste)(a,  /rn2)] t 


Figure  9 


Normalized  solidification  time  as  a  function  of  K/Bi 
and  initial  undercooling.  This  plot  is  independent  of 
particle  radius  and  the  particular  values  of  K«  and  h. 
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is  presented  in  Figure  9.  The  calculated  tim^s  have  been 
normalized  by  the  solidification  time  for  a  no-undercooled 
situation,  (t^)NU  given  by  Equation  (5.10),  making  the  plotted 
values  independent  of  particle  size  and  the  individual  values  of 
K  and  Bi.  The  effect  of  an  exponential  kinetic  model  would  be  to 
shift  the  curves  to  the  right  (higher  time  values)  but  it  is  only 
noticeable  in  the  low  range  of  K/Bi  (<10  )  and  high 
undercoolings . 

Figure  9  shows  that  as  K/Bi  increases  the  curves  approach 
asymptotically  the  value  predicted  from  the  limiting  case  of 
adiabatic  solidification 

(^f )  u  =  (1  _  ^f^NU  (5.39) 

which  implies  an  instantaneous  recalescence  after  nucleation.  As 
the  rate  of  the  heat  flow  increases  (K/Bi  is  reduced)  the  time 
values  become  larger  than  those  predicted  by  Equation  (5.39). 
Examples  of  calculated  solidification  times  for  5  yim  A1  droplets 
are  reported  in  Table  4.  Increasing  rQ  will  proportionally 
increase  the  values  in  this  table. 

Table  5  shows  the  effect  of  the  nucleation  temperature  on 
the  overall  time  scale.  Although  the  solidification  time 
decreases  with  larger  undercoolings,  this  reduction  is  offset  by 
a  corresponding  increase  of  the  cooling  time  (i.e.  the  time 
spent  in  the  liquid  state  before  reaching  the  nucleation 
temperature) .  Therefore,  the  time  it  takes  to  form  a  solid 
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powder  from  a  liquid  droplet  at  or  above  its  melting  temperature 
increases  with  increasing  undercooling  prior  to  nucleation.  This 
is  a  significant  observation  in  determining  whether  the  flight 
time  available  in  a  given  atomizer  will  be  sufficient  for 
complete  solidification. 

Table  4 

Solidification  Times  for  10  yin  Aluminum  Droplets  (ys) 


km 

h 

Nucleation 

Temperature  0N 

(cms~^K“l ) 

(WnT2K-1> 

(ys) 

-0.01 

-0.5 

o 

o 

• 

H 

/ 

2 

2x10s 

12 .8 

15.2 

8.9 

3.2 

2 

10  5 

25.6 

28.3 

15.7 

3.9 

2 

5x10^ 

51.2 

54.2 

29.1 

4.8 

8 

5x10* 

51.2 

51.9 

26.9 

1.8 

50 

5x10* 

51.2 

51.0 

25.9 

0.5 

Table  5 

Effect  of  the  Nucleation  Temperature 
on  the  Time  Scale  ( ys) ^ 


0N 

tt 

fcN 

tN  + 

0 

51.2 

— 

51.2 

-0.01 

54.2 

0.5 

54.7 

-0.25 

41.9 

13.8 

55.7 

-0.50 

29.1 

30.2 

59.3 

-0.75 

16.3 

50.3 

66.6 

-1.0Q 

4.8 

103.7 

108.5 

2rQ  =  10 

ym,  tq  = 

h  =  5xl04 

Wm'^'K'1 , 

km 

=  2  cms~ 

V1 
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5 . 3 . 2 . 2  Solidification  from  massive  nucleation 

The  net  effect  of  increasing  the  number  of  nucleation  events 
(or  2rg/dMX)  manifests  itself  as  an  apparently  higher  K/Bi  due  to 
the  increase  in  the  interfacial  area,  which  in  turn  increases  the 
solidification  rate  (dg/dt) .  Therefore,  the  thermal  history  is 
now  dependent  on  0N  and  the  product  (2r q/<5^)  (K/Bi)  and  must  be 
unique  as  long  as  these  parameters  remain  constant.  It  should  be 
noted  that  every  twofold  increase  in  2rg/dMX  translates  into  an 
eight-fold  increase  in  the  number  of  grains. 

Figure  10  shows  the  enthalpy-temperature  curves  for  a  fixed 
0N  =  -0.5  and  several  values  of  the  product  (2r0/dMX) (K/Bi) .  If 
one  takes  a  low  K/Bi  value  of  50  (e.g.  KM  =  2  cm  s'1  K_1  and 
h  =  4  x  105  W  m-2  K-1)  these  curves  will  represent  the  thermal 
histories  of  droplets  with  1,  8,  64  and  1000  nucleation  events 
respectively.  As  expected,  increasing  the  number  of  grains 
reduces  the  role  of  the  heat  transfer  during  recalescence ,  and 
consequently  the  extent  of  the  "rapid  solidification"  regime. 
However,  the  occurrence  of  multiple  nucleation  does  not  have  a 
critical  effect  at  values  of  K/Bi  >  10  because  the  role  of  the 
heat  flow  is  already  small  even  for  solidification  from  a  single 
nucleus. 

Although  numerical  calculations  for  this  case  were  not 
performed  due  to  the  geometrical  complexity  of  the  problem,  it  is 
anticipated  that  the  validity  of  the  Newtonian  approximation  is 
enhanced  as  the  number  of  nuclei  increases.  This  is  due  to  the 


DIMENSIONLESS  ENTHALPY  ip  =  (H-HCM)/AH 


TEMPERATURE  (K) 

700  800  900 


DIMENSIONLESS  TEMPERATURE  0  =  CL(T-TM)/AHM 


Figure  10.  Effect  of  the  number  of  nucleation  events  on  the 

thermal  history  of  undercooled  droplets  solidifying 
under  Newtonian  conditions. 


reduction  in  the  scale  of  the  heat  transfer  problem  (from  r^  to 
^mx/2)  which  hinders  the  development  of  temperature  differences 
inside  the  droplet.  Furthermore  it  is  expected  that  other 
solidification  morphologies  which  increase  the  liquid  solid 
interfacial  area  (and  therefore  dg/dt)  such  as  dendritic  vs. 
plane  front  structures,  should  have  an  effect  on  the  thermal 
history  similar  to  the  occurrence  of  multiple  nucleation. 
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CHAPTER  6 

LIQUID  COOLING  STAGE:  NON-NEWTONIAN  SOLUTION 


The  following  section  is  devoted  to  the  analysis  of  the 
pre-solidification  period  allowing  for  a  non-uniform  temperature 
distribution  to  develop  inside  the  droplet.  We  consider  the 
cooling  of  a  superheated  liquid  sphere  from  an  initially  uniform 
temperature  0q  until  its  surface  reaches  some  prescribed 
nucleation  temperature  0N.  Assuming  that  only  conductive  heat 
transfer  takes  place  in  the  metal  and  that  the  surface 
conductance  is  constant  in  space  and  time,  we  obtain  a 
concentric-isotherm  temperature  distribution  given  by  [88] 

(0+Ste)  / (  0Q+ste)  =  2Bi  £  Lp [sin  ( Ap<J>) /<j>]  exp(-A2Fo)  (6.1) 

P 

Lp  =  sin(Ap)/[A^  -  (1-Bi) sin2 ( Ap) ]  (6.2) 

where  4>=  (r/rQ)  is  the  dimensionless  radial  position  and  A  are 
the  eigenvalues  defined  by 

Ap  cos(Ap)  -  (1-Bi)  sin ( A  )  =  0  (6.3) 


The  analytical  solution  given  above  can  now  be  used  to  evaluate 
the  characteristic  parameters  of  this  period. 
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It  has  been  shown  elsewhere  that  the  achievable  supercooling 

in  a  liquid  is  determined  by  the  activity  of  the  heterogeneous 

nucleants  147-55]  and  the  cooling  rate  [46,74].  The  expression 

for  the  instantaneous  average  cooling  rate  obtained  from 

Equation  (6.1)  was  found  to  be  identical  to  that  obtained  from 

the  Newtonian  model  when  the  surface  temperature  T  is 

s 

substituted  for  T  in  Equation  (5.5) 

eL  =  3  (ctL/kL)  (h/r0)  (Ts  -  Tq)  (6.4) 

Examples  of  calculated  cooling  rates  for  50  ym  droplets  of  Al,  Fe 
and  Ni,  assuming  a  high  heat  transfer  coefficient  of 
5  x  104  W  m-2  K-1  are  given  in  Table  6. 


Table  6 

-l  + 

Average  Cooling  Rates  in  50  ym  Liquid  Metal  Droplets  (MK  s  ) 

Al  Fe  Ni 

Biot  Number,  hrQ/kL  0.014  0.031  0.032 

Initial  cooling  rate 

(0S  =  0O  =  0. 2Ste)  1.7  1.9  1.7 

at  the  melting  point 

(6S  =  0M  =  0)  1.4  1.6  1.4 

at  the  onset  of  hypercooling 
(%  =  0) 

t  h  =  5xl04  Wm”2K~1 


* ' 


0.6 


1.2 


1.0 


I 
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Cooling  times  can  also  be  estimated  from  Equation  (6.1) 
setting  6  =  0N,  <J>  =  1,  and  looking  for  the  corresponding  Fo  value 
via  Newton-Raphson  or  any  other  root  searching  technique. 

Figure  11  shows  the  dependence  of  the  cooling  time  for  a  droplet 
nucleating  at  its  melting  point,  tM,  on  the  Biot  number  and 
initial  superheat.  The  time  calculated  from  the  analytical 
solution  has  been  normalized  for  purposes  of  representation, 
dividing  it  by  the  Newtonian  cooling  time  obtained  from 
Equation  (5.7).  The  values  in  the  ordinate  of  Figure  11  indicate 
the  divergence  between  the  two  models  and  may  be  interpreted  as 
correction  factors  to  the  Newtonian  equation. 


Table  7 

Cooling  Times  for  Atomized  Metal  Droplets  (us) 


0  0/Ste 

2r0  (ym) 

Al 

Fe 

Ni 

0.3 

50 

41 

84 

87 

(42) 

(90) 

(93) 

500 

304 

445 

453 

(417) 

(900) 

(935) 

0.3 

50 

114 

243 

252 

(115) 

(248) 

(257) 

500 

1053 

2002 

2063 

(1148) 

(2476) 

(2572) 

t  h  =  5xl04 

Wm“2K_1 

Values  in 

parenthesis 

are  from 

the  Newtonian  model. 

Examples  of 

cooling  times 

calculated  from 

Figure  11  for 

droplets  of  Al,  Fe  and  Ni  are 

given  in 

Table  7. 

As  expected, 

f 


NORMALIZED  COOLING  TIME  [3Bi / 1 n ( l+0n/Ste )]  • (a,  /rn2)t 
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Figure  11.  Normalized  time  for  the  surface  of  a  liquid  droplet  to 
cool  from  an  initial  temperature  TQ  to  its  melting 
temperature  TM  by  convection  to  an  environment  with 
constant  temperature  TG,  as  a  function  of  Biot  number 
and  initial  superheat. 


J 
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equivalent  initial  superheats  the  time  for  non-Newtonian  cooling 
is  shorter  because  not  all  the  sensible  heat  in  the  droplet  is 
removed  before  its  surface  reaches  the  nucleation  temperature. 
Note  also  that  reducing  the  droplet  diameter  will  result  in  a 
smaller  reduction  in  tM  than  that  predicted  by  the  Newtonian 
model.  On  the  other  hand,  increasing  the  heat  transfer 
coefficient  will  have  a  larger  effect  in  reducing  the  cooling 
time  than  the  inverse  proportionality  expected  from 
Equation  (5.7). 

Normalized  temperature  distributions  in  liquid  droplets  when 
their  surfaces  reach  the  melting  point,  for  different  Biot 
numbers  and  initial  temperatures  are  shown  in  Figure  12.  The 
largest  gradient  occurs  at  the  surface  and  is  given  by  the 
boundary  condition 
) 

OT/3r)s  =  (h/kL)  (Ts  -  Tg)  (6.5) 

Appropriate  substitution  of  the  variables  in  this  equation  reveal 

that  significant  gradients  may  exist  within  the  droplet  even  for 

Biot  numbers  of  0.01  or  lower.  For  example,  a  heat  transfer 

coefficient  of  5  x  10^  W  m-^  (Bi=0.014  for  2r q  =  50  ym) 

5  -1 

produces  a  maximum  gradient  of  3.5  x  10  Km  in  A1  droplets 
when  Ts  =  Tm. 

Although  the  surface  gradient  is  independent  of  superheat 
and  particle  size,  as  concluded  from  Equation  (6.5),  these 
variables  have  an  important  effect  on  the  internal  temperature 


Figure  12.  Normalized  temperature  distributions  in  a  liquid 
droplet  when  its  surface  reaches  the  melting 
temperature  TM,  for  different  Biot  numbers  and  initial 
superheats. 
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profiles.  For  example,  maximum  temperature  differences  A(T)  of 
4.4  and  45  K  between  the  center  and  the  surface  were  estimated 
from  Equation  (6.1)  for  50  and  500  pm  Aluminum  droplets 
respectively,  using  the  h  and  T  values  given  above  and  an 
initial  superheat  of  190  K.  Thus,  the  basic  assumption  that 
temperature  differences  inside  the  body  are  negligible  (and  hence 
that  the  Newtonian  model  is  a  valid  approximation)  for  Bi  <  0.1 
may  not  always  be  justified.  The  original  superheat  was  found  to 
influence  only  slightly  the  temperature  distributions  for 
Bi  <  0.01. 

During  metal  atomization  one  wishes  to  suppress  nucleation 
until  a  certain  level  of  enthalpy  is  reached  in  the  liquid,  since 
this  enthalpy  is  essentially  a  measure  of  the  volume  fraction 
that  can  be  solidified  before  the  droplet  recalesces  to  the  "slow 
growth"  regime,  as  indicated  in  Chapter  5.  Because  the  surface 
is  always  cooler  than  the  rest  of  the  droplet,  the  volume 
averaged  enthalpy  of  the  latter,  <jj ,  is  always  larger  than  that 
corresponding  to  the  surface  temperature.  ijJ  is  related  to  the 
average  temperature  0  through  Equation  (3.10)  for  g  =  0, 

(ip  =  1+0)  ,  and  can  be  calculated  from  the  following  expression 

(ijj-l+Ste)/ (eQ+Ste)  =  6Bi2  £  L'  exp(-ApFo)  (6.6) 

P 

L'  =  [sin ( Ap) /A2]  Lp  (6.7) 

We  now  define  <J*N  as  the  average  enthalpy  of  the  droplet  at  the 
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moment  of  nucleation;  represents  the  minimum  amount  of  heat 
that  must  be  withdrawn  by  the  cooling  environment  during 
solidification.  Heat  extraction  after  nucleation  will  then  be 
necessary  in  all  cases  where  >  0. 


Table  8 

Undercooling  Required  at  the  Surface  for  the^. 
Droplet  to  Reach  the  Hypercooled  Regime,  (K) 


Diameter  (urn) 

Al 

Fe 

Ni 

500 

372 

415 

474 

50 

365 

353 

418 

Newtonian 

364 

345 

411 

t  h  =  5xl04  Wm  2K-1,  9q  =  0.2Ste 


Combining  Equations  (6.1)  and  (6.6)  yields  a  relationship 
between  the  desired  value  and  the  required  0N  at  the  surface, 
(<J>  =  1)  ,  which  represents  the  maximum  allowable  nucleation 
temperature.  This  relationship  is  plotted  in  Figure  13  as  a 
function  of  the  Biot  number  for  two  levels  of  enthalpy.  =  1 

is  the  point  where  all  the  superheat  has  been  removed  from  the 
droplet,  and  =  0  is  the  enthalpy  where  the  droplet  reaches  the 
hypercooled  regime.  Values  of  the  undercooling  required  at  the 
surface  when  =  0  for  droplets  of  Al,  Fe  and  Ni  are  reported  in 
Table  8.  Note  that  as  r^  increases,  nucleation  must  be  suppressed 
to  lower  temperatures  in  order  to  reach  the  hypercooled  regime. 


r 
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but  the  cooling  rates  (and  the  likeliness  to  achieve  larger 
undercoolings  prior  to  solidification)  decrease  in  an  inverse 
proportion,  as  indicated  by  Equation  (6.4).  One  can  then  conclude 
that  the  achievement  of  large  undercoolings  is  not  only  hindered 
kinetically  but  also  thermally  when  the  particle  size  increases. 
The  effects  are  more  noticeable  in  Fe  and  Ni  due  to  their  lower 
thermal  conductivity. 

The  time  required  to  reach  an  enthalpy  ij]  can  be  calculated 
from  Equation  (6.6)  using  the  Newton-Raphson  formulation  as 
before.  Results  for  A1  are  reported  in  Figure  14,  where  the 
values  have  been  normalized  by  their  Newtonian  equivalent,  (see 
Equation  (5.7)  for  0N  =  i(/N-l).  The  times  calculated  from  the 
series  solution  are  now  higher  than  the  corresponding  Newtonian 
values  since  the  ratio  of  the  cooling  rates  at  the  same  average 
enthalpy  is  given  by 

e (Newtonian) /e (non-Newtonian)  =  (Q  +  Ste)/(6  +  Ste)  (6.8) 

s 

and  this  ratio  is  always  larger  than  one  for  0  =  0  .  For  example, 

the  non-Newtonian  and  Newtonian  times  to  reach  the  hypercooled 

.  4  —2  —1 

regime  in  a  50  pm  droplet  of  A1  with  h  =  5  x  10  WmK  and  an 
initial  superheat  of  190  K  are  490  and  455  ;is  respectively. 


0  0.2  0.4  0.6  0.8  1.0 

BIOT  NUMBER  (Bi  =  hrQ/kL) 


Figure  14.  Normalized  cooling  time,  for  the  droplet  to  reach  an 
average  enthalpy  content  of  as  a  function  of  Biot 
number . 
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CHAPTER  7 


83 


SOLIDIFICATION  WITHOUT  UNDERCOOLING:  ENTHALPY  MODEL 


We  now  address  the  solidification  problem  where  no 
significant  undercooling  is  required  for  either  nucleation  or 
growth.  This  situation  resembles  the  case  where  active  catalysts 
limit  the  achievable  undercooling  to  a  few  degrees.  As  in  the 
previous  chapter,  we  assume  that  only  conductive  heat  flow  takes 
place  inside  the  droplet  and  that  the  heat  transfer  coefficient 
is  constant  in  time  and  space,  resulting  in  a  concentric  isotherm 
geometry.  Since  all  the  points  on  the  surface  have  the  same 
enthalpy  (or  temperature) ,  solidification  must  sta  t 
simultaneously  as  a  thin  shell,  and  the  interface  fo) lows  at  any 
time  the  spherical  shape  of  the  isotherms.  We  further  assume 
that  the  interface  configuration  is  stable  through  the  entire 
solidification  process  (no  breakdown  to  cells  or  dendrites). 


7.1  Mathematical  Formulation 

For  this  part  of  the  problem  the  spherical  droplet  is 
divided  into  a  fixed  number  of  concentric  shell  domains  of 
thickness  Ar  for  the  internal  nodes  (1  <  i  <  n)  and  Ar/2  for  the 
center  (i  =  1)  and  the  surface  (i  =  n)  nodes,  as  shown  in 
Figure  15.  An  enthalpy  balance  for  an  element  of  volume  and 

surface  area  A^  can  be  written  as 


i 
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(d/dt)  /  p  H  dV  =  /  (k  VT)  •  (e  dA)  (7.1) 

Vi  Ai 

where  e  is  the  unit  vector  normal  to  the  bounding  surface  of  the 
element.  The  resulting  expression  for  an  internal  node  is 

(dH./dt)  =  k+ (At/vp  (3T/3r)t  -  k“  (A'/VT)  ( 3T/3r )  7  (7.2) 

where  a7,  a+  and  denote  the  inner  and  outer  surface  areas  and 
the  volume  of  the  spherical  shell  i  respectively.  The 
superscripts  (+)  or  (-)  refer  to  the  positive  and  negative  r 
directions . 

In  general,  in  Equation  (7.2)  refers  to  the  specific 
enthalpy  of  the  volume  element  i  and  assumes  different  forms  when 
the  node  is  in  the  liquid,  the  solid,  or  contains  the 
liquid-solid  interface.  k  is  the  thermal  conductivity  of  the 
liquid  or  the  solid;  k+  =  k-  for  all  nodes  except  that  containing 
the  interface,  i  =  2,  where  k+  =  kg  and  k  =  k^  (see  Figure  15). 

Equation  (7.2)  can  now  be  put  in  terms  of  the  dimensionless 
nodal  enthalpy  ij^,  temperature  fK,  distance  <p  and  time  Fo,  and 
becomes 


(d-lj/dFo)  =  ( r q /V A )  (Qj  -  Qt) 


(7.3) 


where 


Q\  =  ~  (k+/kL) (At/rJ)  (30/30)^ 


i  /  n  (7.4) 


/ 
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Q~  =  -  (k  /kL)  (Ai/rjj)  (30/;>4>)"  i?l  (7.5) 

For  the  center  node,  i  =  1,  the  term  Q”  on  the  right  hand  side  of 
Equation  (7.3)  becomes  zero.  For  the  surface  node,  i  =  n,  an 
energy  balance  similar  to  expression  (7.2)  is  written,  modifying 
the  (+)  term  to  account  for  the  heat  loss  to  the  environment. 
Appropriate  manipulation  of  this  equation  using  the  definitions 
of  dimensionless  nodal  enthalpy  and  temperature  yields  the 
following  expression  for  Q* 

Qn  =  (A^Ao>  Bi  ( en  +  Ste)  (7.6) 

As  concluded  from  the  discussion  of  the  enthalpy-temperature 
diagram.  Figure  2,  the  dimensionless  nodal  enthalpy  is  always 
negative  in  the  solid  region  (g^^  =  1,  0^  <  0)  ,  whereas  in  the 
superheated  liquid  region  iJj^  is  larger  than  unity,  (g^  =  0, 

>  0) .  A  volume  element  containing  the  liquid  solid  interface 
(i  =  f.)  is  at  the  melting  point  of  the  metal;  therefore,  the 
value  of  C-'j,  is  equal  to  the  volume  fraction  of  the  element 
which  is  in  the  liquid  state,  1-g^. 

Equations  (7.3)  to  (7.6)  were  put  into  their  equivalent 
fully-implicit  finite  difference  forms.  The  resulting  system  of 
simultaneous  equations  for  the  different  nodes  was  solved  using 
the  Gauss-Seidel  iterative  technique  starting  with  the 
temperature  distribution  obtained  from  the  analytical  solution 


t 


for  the  pre-solidification  stage.  Equation  (6.1),  when  the 


surface  reaches  the  melting  temperature  (0^  =  0) .  The 
convergence  criterion  used  for  each  iteration  cycle  was 


i=n+l 

I 


i=l 


('P 


t 

i ,  k+1 


/;'l,k+l> 


1  <  10 


(7.7) 


where  the  superscripts  t  and  c  denote  the  trial  and  calculated 
values  respectively. 

The  dimensionless  time  interval,  AFo,  was  adjusted  in  every 

step  k  to  give  roughly  uniform  displacements  of  the  liquid-solid 

interface.  Trial  runs  were  performed  changing  the  number  of 

nodes  N,  and  the  length  of  the  interface  displacements  AcJ> * ,  (i.e. 

the  length  of  the  time  steps) .  It  was  observed  that  for  the 

metals  investigated,  increasing  the  number  of  nodes  above  21  does 

not  significantly  affect  the  results.  On  the  other  hand,  some 

parameters  like  Gl/r,  cooling  rates  in  the  liquid  and  interface 

★  * 

velocities  at  the  beginning  ($  -»•  1)  and  at  the  end  ($  -+  0)  of 

solidification  are  particularly  sensitive  to  the  time  step  size, 
especially  in  the  low  range  of  Biot  numbers.  For  these 
calculations  the  time  step  sizes  were  appropriately  selected 
until  convergence  was  reached.  Calculations  were  performed  on  a 
CDC  CYBER  175  computer. 


7.2  Results  and  Discussion 

The  computer  heat  flow  model  was  used  to  calculate  the 
important  cooling  and  solidification  parameters  (e.g.  cooling 
rates,  interface  velocities,  temperature  gradients,  etc.)  as  a 
function  of  time  and  the  dimensionless  var iables--mainly  Biot 
number--governing  the  rate  of  heat  extraction  from  the  droplets. 

Curves  relating  the  net  solidification  time  for  droplets  of 
Al,  Fe  and  Ni  to  the  Biot  number  and  dimensionless  initial 
superheat  are  presented  in  Figure  16.  The  time  resulting  from  the 
numerical  solution  of  the  Enthalpy  model  has  been  normalized 
dividing  it  by  the  time  calculated  for  solidification  in  a 
Newtonian  regime,  given  by  Equation  (5.10). 

Table  9 

Net  Solidification  Times  for  Atomized^ 

Metal  Droplets:  no  Undercooling,  (ys) 


2r0  (urn) 

90/Ste 

Al 

Fe 

Ni 

50 

0.3 

263 

274 

317 

0 

257 

265 

305 

Newtonian 

252 

216 

283 

500 

0.3 

3172 

5405 

4898 

0 

3047 

4795 

4267 

Newtonian 

2518 

2160 

2826 

f  h  =  5x10^  Wm-^K-^ 

Examples  of  calculated  solidification  times  for  droplets  of 
the  three  metals  investigated  are  shown  in  Table  9  and  compared 


BIOT  NUMBER  (Bi  =  hrQ/kL) 


Normalized  solidification  time  for  liquid  droplets  of 
aluminum,  iron  and  nickel,  as  a  function  of  Biot 
number  and  dimensionless  initial  superheat  Oq.  The 
calculated  curve  for  nickel  with  no  superheat  closely 
follows  that  of  aluminum  for  0  =  0.3Ste. 


to  the  corresponding  Newtonian  predictions.  As  anticipated, 
increasing  initial  superheat  prolongs  solidification;  the 
sensible  heat  retained  in  the  liquid  portion  of  the  droplet 
increases  resulting  in  longer  times  from  initiation  to  completion 
of  solidification.  Note  also  that  the  effect  of  superheat  on  the 
solidification  time  diminishes  with  decreasing  Biot  number,  in 
agreement  with  our  previous  observation  that  0  has  a  negligible 
influence  on  the  temperature  distributions  at  the  start  of 
solidification  for  Bi  <  0.01. 

Figure  17  shows  the  effect  of  Biot  number  and  initial 
superheat  on  the  liquid-solid  interface  velocity  in  A1  droplets 
as  solidification  progresses  (i.e.  <p*  decreases).  As  in 

previous  figures,  the  values  reported  are  the  ratios  of 
non-Newtonian  to  Newtonian  velocities  calculated  for  the 
concentric  geometry.  Equation  (5.13).  Note  that  in  all  cases  the 
Newtonian  results  represent  an  upper  boundary  to  the  achievable 
growth  velocites  predicted  by  the  Enthalpy  model. 

Some  of  the  trends  observed  in  this  figure  can  be  better 
understood  from  the  boundary  condition 

R  =  (Q+  -  Q~) /ahm  =  (ksGs  -  kLGL)/AHM  (7.8) 

which  relates  the  interface  velocity  R  to  the  heat  fluxes  on  the 
liquid  (Q“)  and  solid  (Q+)  sides  respectively.  In  the  absence  of 
initial  superheat  (fig  =0),  both  models  dictate  the  same  velocity 
at  the  start  of  solidification  (<{>*  =  1)  regardless  of  Biot 


0.01  0.1  1.0 


BIOT  NUMBER  (Bi  =  hrQ/kL) 


Figure  17.  Normalized  interface  velocity  as  a  function  of  Biot 

number  and  initial  superheat.  Note  that  the  Newtonian 
model  gives  always  the  limit  to  the  achievable 
velocity  as  calculated  from  the  enthalpy  model. 
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number.  This  results  arises  when  no  gradient  GL  exists  prior  to 
solidification  (hence  Q'  =  0  in  the  above  equation),  and  the 
Q+  term  is  given  by  the  surface  heat  flow 

Q+  =  h  (Tm  -  Tg)  (7.9) 

which  is  identical  for  the  Newtonian  and  numerical  solutions. 

For  a  finite  superheat  ('  >  0)  there  is  a  positive  temperature 

gradient  GL  in  front  of  the  interface  which  slows  down  growth  by 
opposing  the  Gg  term  in  Equation  (7.8).  For  example,  an  initial 
superheat  of  190  K  (‘q  =  0.3Ste)  may  cut  the  initial  velocity  to 
less  than  half  the  value  calculated  for  =  o  (see  Figure  17) . 

It  is  also  clear  that  this  effect  decays  rapidly  as  fraction 
solid  increases ,  due  to  the  elimination  of  the  liquid  residual 
superheat  in  the  early  stages  of  solidification. 

Examples  of  interface  velocities  calculated  with  the 
Enthalpy  model  for  Al,  Fe  and  Ni  dropTets,  and  their 
corresponding  Newtonian  values  are  given  in  Table  10.  As 
anticipated,  the  two  solutions  diverge  with  increasing  Biot 
number  and  the  progress  of  solidification.  Furthermore,  the 
Enthalpy  model  predicts  a  reduction  in  velocity  as  the  particle 
size  increases,  whereas  the  Newtonian  model  indicated  no  such 
effect.  On  the  other  hand,  increasing  the  heat  transfer 
coefficient  at  constant  rQ  will  increase  R  at  any  location. 


93 


Table  10 

Interface  Velocities  for  Concentric  Solidification 
Enthalpy  Model:  No  Undercooling,  (cm  s-*)" 


•P  :(g) 

2r0  (i'm) 

A1 

Fe 

Ni 

0.9 

500 

4.0 

3.7 

3.0 

(0.271) 

50 

4.1 

4.7 

3.6 

Newtonian 

4.1 

4.8 

3.6 

0.5 

500 

10.0 

5.3 

6.5 

(0.875) 

50 

13.0 

11.3 

10.5 

Newtonian 

13.2 

15.4 

11.8 

0.1 

500 

44.5 

12.8 

21.0 

(0.999) 

50 

107.0 

46.6 

61.1 

Newtonian 

330.9 

385.8 

295.1 

t  h  =  5xl04  Wm  2K-1,  0q  =  0 . 3Ste . 

The  results  of  normalized  ratio  of  temperature  gradient  in 
the  liquid  side  of  the  interface,  GL,  to  interface  velocity  R  for 
the  three  metals  investigated,  are  presented  in  Figure  18.  G  at 

Li 

the  beginning  of  solidification  is  approximately  given  by  the 
boundary  condition  at  the  surface 

lim  Gl  =  <h/kL)  (Tm  -  Tg)  (7.10) 

/-I 

while  the  initial  interface  velocity  is  close  to  that  calculated 
from  Equation  (5.13),  hence,  the  initial  Gl/r  ratio  should 
asymptotically  approach  a  constant  value  with  decreasing  Biot 


number 


BIOT  NUMBER  (Bi  =  hrQ/KL) 


Figure  18.  Normalized  ratio  of  the  temperature  gradient  in  the 
liquid  at  the  liquid-solid  interface,  GL,  to  the 
interface  velocity  R,  during  solidification  of  Al, 
and  Ni  droplets. 
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lim  (Gl/R)(()*=1  =  AHm/(ul  Cl)  (7.11) 

Bi’-O 

For  example,  GL/R  ratios  of  2  x  107  and  4.3  x  107  K  s  m-2  were 
calculated  for  50  and  500  pm  A1  droplets  respectively,  at  the 
start  of  solidification.  The  limiting  value  given  by 
Equation  (7.11)  is  1.1  x  10^  K  s  m-2.  As  noted  in  Figure  18, 

Gl/r  decreases  markedly  with  increasing  fraction  solid, 
especially  at  low  Bi,  since  the  liquid  superheat  is  rapidly 
dissipated  (GL  -+  0)  and  the  interface  velocity  increases. 

Normalized  values  of  the  calculated  average  cooling  rates  in 
the  solid  and  liquid  during  solidification  of  Aluminum  droplets 
are  shown  in  Figure  19  and  Figure  20  respectively.  Note  that  the 
part  of  the  normalizing  factor  inside  the  brackets  is  the 
Newtonian  solidification  time  obtained  from  Equation  (5.9),  which 
is  directly  proportional  to  r^  and  inversely  proportional  to  h. 
Therefore,  increasing  the  Biot  number  by  increasing  the  droplet 
size  Tq,  will  increase  the  actual  cooling  rate  only  if  the  slope 
of  the  curve  at  that  point  is  larger  than  unity.  On  the  other 
hand,  increasing  the  heat  transfer  coefficient  h  will  always 
increase  the  cooling  rate.  Finally,  decreasing  the  initial 
superheat  will  increase  the  achievable  cooling  rates  in  the 
solid,  but  the  effect  is  only  felt  at  the  initial  stages  of 
solidification.  The  average  cooling  rate  in  the  liquid  ,  ~L, 
decreases  with  the  progress  of  solidification,  the  effect  being 
more  pronounced  at  small  Biot  numbers  (Figure  20) . 


NORMALIZED  AVERAGE  COOLING  RATE  IN  THE  SOLID 
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BIOT  NUMBER  (Bi  =  hrQ/kL) 


Figure  19.  Normalized  instantaneous  average  cooling  rate  in  the 

solid  during  solidification  of  aluminum  droplets,  as  a 
function  of  Biot  number  and  fractional  position  of  the 
liquid-solid  interface  4>*.  0o  and  Ste  denote  the 
dimensionless  initial  temperature  and  Stefan  number, 
respectively. 
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Figure  20.  Normalized  instantaneous  average  cooling  rate  in  the 
liquid  during  solidification  of  aluminum  droplets,  as 
a  function  of  Biot  number  and  fractional  position  of 
the  liquid-solid  interface  (J>*.  ()g  and  Ste  denote  the 

dimensionless  initial  temperature  and  Stefan  number, 
respectively. 
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Table  11 

Average  Cooling  Rates  During  Solidification 
of  Aluminum  Droplets  (K  s-1)+ 

.  Diameter  (ym) 


<p 

50 

500 

1.0 

7.3  x  105 

1.1  x  10 

0.9 

2.2 

1.2  x  10 

0.7 

- 

1.2 

0.7 

22.0  , 

5.3  x  10 

0.3 

2.4  x  10* 

6.1  x  10 

0 

1.3  x  106 

1.2  x  10 

t  h  =  5x10^  Wm-^K  6q  =  0.3Ste 

Examples  of  cooling  rates  during  solidification  of  A1 
droplets  without  undercooling  are  given  in  Table  11.  It  is 
important  to  mention  that  the  relatively  low  values  obtained  from 
this  model  are  due  in  part  to  the  prescription  of  a  discrete 
melting  temperature.  When  a  liquid-solid  temperature  range 
exists,  as  in  the  case  of  an  alloy,  and  no  substantial 
undercooling  is  considered,  the  solidification  times  will  be 
similar  to  those  of  the  pure  metal  (provided  the  heat  of  fusion 
is  reasonably  constant) .  Since  solidification  of  an  alloy  is  not 
isothermal,  the  cooling  rates  should  therefore  be  higher 
throughout  the  process. 


/ 
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CHAPTER  8 

SOLIDIFICATION  WITH  UNDERCOOLING:  ENTHALPY  MODEL 

This  final  chapter  of  the  non-Newtonian  heat  flow  analysis 
deals  with  the  modelling  of  the  solidification  process  in  an 
undercooled  droplet  from  a  single  nucleation  event  occur ing  at 
its  surface.  The  assumptions  of  constant  heat  transfer 
coefficient,  cooling  environment  temperature  and  pure  conduction 
inside  the  droplet  are  held  through  this  model,  but  since  the 
solid  grows  from  one  point  on  the  surface,  the  isotherms  are  no 
longer  concentric  and  the  problem  becomes  two-dimensional. 
Furthermore,  the  thermal  field  is  symmetrical  around  the  growth 
axis,  i.e.  the  line  that  joins  the  nucleation  point  on  the 
surface  with  its  antipode. 

Since  temperature  gradients  may  now  exist  inside  the 
droplet,  the  liquid-solid  interface  is  likely  to  develop 
different  thermal  conditions  at  different  locations.  Therefore, 
the  undercoolings  and  velocities  vary  along  the  interface  and  its 
shape  departs  from  the  spherical  configuration  prescribed  in  the 
Newtonian  analysis  for  axisymmetric  solidification.  Figure  3. 

8.1  Mathematical  Formulation 

We  divide  the  spherical  droplet  into  volume  elements  defined 
on  a  superimposed  bispherical  (bipolar  rotational)  coordinate 
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system  (u,v,w)  [89] .  By  making  the  rotational  axis  coincident 
with  the  growth  axis  and  following  a  pattern  similar  to  the  Wulff 
net  we  obtain  the  mesh  shown  in  Figure  21.  The  coordinate  surface 
u  =  tt/2  represents  the  droplet  surface,  u  =  it  is  the  growth  axis, 
v  =  -  oo  is  the  nucleation  point  and  (it,0)  is  the  droplet  center. 

A  typical  node  i,j  (1  <  i  £  n,  1  £  j  £  m)  with  coordinates 


u.  =  (rr/2)  fl+(i-l)/(n-l)] 


(8.1) 


Vj  =  0.5  ln{cot2[U-(2j-l)/2m)  (tt/2)]} 


(8.2) 


is  depicted  in  Figure  21.  The  volume  element  about  this  node  is 
bound  by  the  coordinate  surfaces 


ui  3  ui+l  =  (tt/2)  [1+  (2i-l) /2  (n-1)  ] 


i  /  n  (8.3) 


un  =  it 


(8.4) 


uj  =  tt/2 


(8.5) 


v1?  = 


vj+l 


=  0.5  ln{cot^[  (1-j/m)  (tt/2)  ]  j  /  m  (8.6) 


VN  =  co 

m 


(8.7) 


vf» 


—  oo 


(8.8) 
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Figure  21.  Mesh  defined  in  bispherical  coordinates  for  the 

Enthalpy  model  of  undercooled  solidification.  The 
mesh  has  rotational  symmetry  about  the  growth  axis, 
and  the  interface  is  approximated  by  a  series  of 
annular  steps  following  a  constant  v  surface  for  each 


— urn-. 


1 


103 


Nta*. 


oW,S  nE,S  ,,  S /0.  _  , 

Bi 1 1  =  B  ± L  =  tanh(v1/2)  =  -1 

also,  for  =  »’ 


bW,N 

i,m 


BE'N 

i  ,m 


tanh(v^/2) 


1 


(8.14) 


(8.15) 


We  now  recall  the  general  expression  for  the  nodal  enthalpy 
balance.  Equation  (7.1),  and  replace  the  gradient 

VT  =  (gl/^II)  OT/Du)  +  (e2/?'^2l)  (aT/3v)  (8.16) 

and  the  differentials  of  the  heat  transfer  areas,  namely 

dA  —  2  m3 3  (8.17) 

for  the  E  and  W  surfaces,  and 

d  A  —  /m mgg  (8.18) 

for  the  N  and  S  surfaces.  m^  and  are  the  metric 

coefficients  defined  by 

mll  =  m22  =  (cosh (v) -cos (u) ] ^  (8.19) 

2  2 
m33  =  [r0sin<u)l  /[cosh (v) -cos (u) ] 


(8.20) 
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and  e  is  given  by  -e^ ,  -§2  and  e 2  for  the  W,  E,  S  and  N 

surfaces  respectively.  Appropriate  rearrangement  and  integration 
of  Equation  (7.1)  yields  the  dimensionless  form  of  the  nodal 
enthalpy  balance  in  bispherical  coordinates 


<a*1(j/dFo)  -  <2»r;S/vi,J>  (Q?tj  +  of,,  ♦  of,,  +  Of,,, 


(8.21) 


where 


of,,  =  -Q?,j+1  -  (kNAL)  Oe/3v)N 


ln{ [cosh (  Vj) -cos (u|) ] /[cosh (Vj) -cos (u^) ] }  j  /  m  (8.22) 


qN  =  0 

wi,m 


(8.23) 


Qi,l  "  0 


(8.24) 


Q?  • 


-Qi+1, j  -  ^E/kL)  (99/3u) E 

[as in (B?' ^sin (uf ) ) -as in (B?'  ffsin  (u?)  )  ] 

1  I  J  -L  If  J  A 


i  /  n  (8.25) 


2E  • 

'n,j 


=  0 


(8.26) 


2^  ■  =  -Bi  (0 •  .+Ste) [tanh (v^) -tanh (v^) ] 
J-/J  Af  J  J  J 


(8.27) 


4>,  0,  Bi,  Ste  and  Fo  are  defined  as  before,  k  in 

Equations  (8.22)  to  (8.27)  is  assigned  the  value  of  kg  or  kL  if 

the  volume  element  is  solid  or  liquid  and  is  surrounded  by  nodes 


# 


I 
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of  the  same  phase.  If  either  the  volume  element  or  the  nearest 
neighbor  involved  in  the  particular  Q  term  under  evaluation  is 
undergoing  solidification,  or  if  they  are  not  of  the  same  phase, 
k  is  asigned  a  weighted  average  between  kg  and  k^. 

Writing  Equation  (8.21)  in  its  implicit  finite-difference 
form  we  obtain  an  expression  of  the  type 


^i, j , k+1  +  ci , j , k+1  °i , j , k+1  “  RHS 


(8.28) 


The  right  hand  side  RHS  of  this  equation  is  a  function  of  the 
nodal  enthalpy  in  the  previous  time  step  \p.  .  the  temperatures 

of  the  surrounding  nodes  and  the  conductivity  ratios  k/kL  in  the 
current  time  step.  The  coefficient  of  0.  .  .  .  in  the  left  hand 

1 t J t K+1 

side  is  also  unknown  for  the  time  Fok+1  since  it  contains  a 
series  of  ( k/k L)  factors.  Substituting  the  definition  of 
enthalpy.  Equation  (3.10),  into  (8.28)  we  obtain  the  expression 
for  the  nodal  temperature 


>i,J,k+l  -  lRHS-1+5i,j,k+l)/<ci,j,k+l+1+AC9i,j,k+l)  <8-29' 


For  a  mesh  with  N  =  n  x  m  nodes  there  are  N  equations  like  this, 
each  one  with  two  unknowns  (o  and  g) .  Since  g  is  known  for  the 
nodes  that  are  completely  liquid  (g  =  0)  or  solid  (g  =  1)  at  time 
FO|t+^ ,  we  only  need  to  define  a  set  of  relationships  between  g 
and  f<  for  those  nodes  undergoing  solidification  at  a  particular 
time.  As  in  Chapter  5,  it  is  necessary  to  prescribe  a  certain 


interfacial  geometry  for  each  node  in  order  to  translate  the 
linear  displacements  dictated  by  the  kinetic  model  into  changes 
in  nodal  fraction  solid. 

We  now  assume  that  the  interface  is  divided  into  a  series  of 
n  annular  steps,  each  one  propagating  along  a  line  of  nodes  of 
constant  u^f  as  indicated  in  Figure  21.  The  interface  takes  the 
shape  of  a  constant  v  coordinate  surface  for  each  level  i,  which 
is  identified  by  v*.  The  dimensionless  interface  velocity  is  then 
j# 

R^  =  (cosh (v^) -cos (u ^ ) ]  1  (3v^/ >Fo)  (8.30) 


and  for  a  time  interval  AFo  the  initial  and  final  positions  of 
the  interface  are  related  through  the  following  equations 


vi,k+l  =  ln[-b+(b2  +  01/21 

b  =  Bi,k+1  cos(Ui)/(l-B*fk+1) 

c=  (1+Bifk+i)/(l"Bi,  k+J.) 

*  r  * 

Bi,k+1  =  esc (u^) sin(RAFosin (u^) +asin [B^ , k 

*  ★  * 

B^  =  sinh(vi)/[cosh(vi)-cos(ui) ] 


i  /  n 


sin  (u^  ]  } 


(8.31) 

(8.32) 

(8.33) 

(8.34) 

(8.35) 


For  the  nodes  located  about  the  growth  axis,  i  =  n 
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vn, k+1  =  ln[(l+B;/k+i)/(l-B^k+l)]  (8.36) 

Bn,k+1  =  *'F°  +  Bn,k  (8.37) 

B*  =  tanh (v*/2)  (8.38) 

Note  that  for  the  initial  step  B*  .  =  -i  for  all  nodes. 

1  /  j 

The  nodal  fraction  solidified  can  now  be  easily  calculated 

*  N 

from  Equation  (8.9)  to  (8.13)  by  substituting  v^  for  Vj  and 
dividing  the  resulting  solid  volume  by  the  total  nodal  volume 

Vi  /  j  * 

Since  the  interface  velocity  R.  and  its  undercooling  Bi  are 
related  by  a  kinetic  law  of  the  type  discussed  in  Chapter  2,  the 

last  step  in  the  analysis  is  to  establish  the  relationship 

* 

between  and  the  temperature  of  the  solidifying  node 

* 

0 i ,  £  ( i )  .  If  one  simP1y  assigns  to  0^  the  value  of 

problems  arise  because  of  the  steep  negative  temperature  gradient 
developed  in  front  of  the  growing  solid.  When  the  interface 
moves  from  a  node  i,?(i)  to  i,>(i)+l.  Figure  22,  it  encounters  a 
much  cooler  environment  (hence  a  larger  driving  force)  because 
the  temperature  tb  p ^ is  usually  much  lower  than  0^  j(i)' 
especially  in  the  early  stages  of  recalescence .  This  condition 
would  translate  into  rather  large  peaks  in  the  velocity 
distribution  every  time  an  interfacial  step  crosses  a  node 
boundary. 


/ 
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A  further  complication  arises  because  the  interface  is  not 
coincident  with  the  same  v  coordinate  surface  for  all  nodes,  see 
Figure  22.  If  such  were  the  case,  most  of  the  latent  heat 
released  by  the  node  i ,  2,  ( i )  would  be  dumped  into  the  liquid  node 
in  front  of  it,  i,£(i)+l.  However,  if  the  interface  in  node 
i,£(i)  crosses  the  v!?  boundary  before  its  neighbors  i-l,£(i-l)  or 
i+l,£(i+l),  it  finds  itself  into  a  node  i,£(i)+l  that  has  two 
extra  heat  sinks,  namely  i-1 ,  £  { i— 1)  +1  and  i+1 ,  £  ( i+1) +1 .  Since 
the  driving  force  is  increased  and  the  heat  can  be  dissipated 
more  easily  for  the  step  i,  the  interface  is  likely  to  develop  a 
"perturbation"  that  grows  faster  into  the  undercooled  liquid  and 
slows  the  growth  of  the  neighbor  steps. 

The  stability  problem  described  above  is  mathematical  in 
nature,  caused  by  the  introduction  of  a  localized  interface  into 
the  Enthalpy  model.  The  behavior  of  the  resulting  "breakdown" 
may  be  compared  to  the  thermal  instability  predicted  for 
interfaces  moving  into  undercooled  liquids  (GL  <  0) ,  as  discussed 
in  Chapter  2. 

It  was  found  after  different  trials  that  the  best  way  to 
maintain  the  interface  stable  was  to  restrict  the  temperature 
gradient  on  the  solid  side,  Gg,  to  positive  values  (or  zero) , 
effectively  making 


ei  "  0i, £ (i) -1 


if  °i,£(i)  <  9 i ,  £  ( i )  - 1 


(8.39) 


9i  =  9 i , £  ( i ) 


if 


0i,£(i)  -  0 i , £ ( i) -1 


(8.40) 
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This  assumption  is  consistent  with  the  physical  situation.  A 
negative  Gg  is  not  likely  to  occur  in  a  recalescing  droplet  since 
it  requires  that  the  liquid  enhances  its  capacity  as  a  heat  sink 
during  solidification,  whereas  in  reality  it  is  exhausted.  The 
above  restriction  on  0*  manifests  as  an  averaging  of  the  growth 
velocities  throughout  the  node,  although  it  has  the  side  effect 
of  slightly  slowing  the  kinetics  of  the  process.  To  partially 
avoid  this  problem.  Equation  (8.39)  was  only  applied  in  the  first 
half  of  the  node,  that  is,  before  the  interface  crosses  the  nodal 
coordinate  v^.  Note  also  that  the  above  restriction  does  not 
affect  directly  the  temperatures  of  the  solid  nodes  behind  the 
interface. 

Equation  (8.29)  can  now  be  solved  with  the  aid  of 
Equations  (8.30)  to  (8.40).  Trial  runs  were  performed  for  84 
(12  x  7) ,  180  (18  x  10)  and  210  (20  x  11)  nodes  in  the  mesh.  The 
convergence  criterion  used  for  each  iteration  cycle  was 


n,m 


'i. j/ei, j 


-  1 


n  t 

l  lgi,£(i)/gi,Ui) 


-  1  |<  10 


-3 


(8.41) 


It  was  found  that  the  results  did  not  differ  significantly  among 
the  three  cases  but  the  computation  time  was  approximately 
quadrupled  for  a  twofold  increase  in  the  number  of  nodes.  The 
computation  time  was  also  increased  by  going  to  a  lower 
undercooling,  a  smaller  particle  size,  heat  transfer  or  kinetic 
coefficients.  It  was  concluded  that  for  the  range  of  values  of 


Ill 


interest  a  mesh  with  84  nodes  produced  satisfactory  results  in  a 
reasonable  time.  As  in  the  previous  cases,  calculations  were 
performed  on  a  CDC  Cyber  175  computer. 

8.2  Results  and  Discussion 

The  computer  model  was  used  to  study  the  solidification  of 
undercooled  liquid  aluminum  droplets,  focusing  on  the  departure 
from  the  Newtonian  trends  established  in  Chapter  5  and  the 
interfacial  behavior  during  solidification.  A  series  of 
calculations  with  both  linear  and  exponential  kinetics,  Equations 
(2.4)  and  (2.5)  respectively,  were  carried  out  for  different 
undercoolings.  Since  temperature  gradients  were  allowed  inside 
the  droplet,  the  enthalpy  temperature  curve  was  evaluated  in 
terms  of  the  weighted  averages  of  these  properties  over  the 
entire  particle. 

As  expected  from  the  Newtonian  analysis,  the  thermal  history 
was  found  to  be  strongly  dependent  on  the  initial  undercooling, 
but  instead  of  being  unique  for  a  constant  K/Bi  ratio,  it  was 
affected  by  the  particular  values  of  the  dimensionless  kinetic 
coefficient  and  the  Biot  number.  The  variables  involved  in  these 
dimensionless  paramenters,  namely  rQ,  km  and  h,  were  therefore 
considered  individually  in  terms  of  their  effect  on  the 
solidification  behavior. 


8.2.1  Effects  of  the  particle  size 

Although  the  K/Bi  ratio  is  independent  of  the  droplet 
dimensions,  increasing  rg  will  progressively  shift  the  ^-6  curve 
away  from  the  Newtonian  case,  as  shown  in  Figure  23.  The  Enthalpy 
model  consistently  predicts  larger  average  undercoolings  than  the 
corresponding  Newtonian  (same  K/Bi)  at  equivalent  fractions 
solidified.  This  is  due  to  the  "local"  recalescence  of  the 
interface  which  reduces  the  growth  rate  at  a  given  g  and 
consequently  enhances  the  role  of  the  heat  flow  during  the  "rapid 
solidification"  stage.  As  the  particle  radius  increases,  the 
external  heat  extraction  is  less  able  to  influence  the  withdrawal 
of  latent  heat  released  at  the  moving  solid-liquid  front.  The 
interface  then  reduces  its  undercooling  at  an  earlier  stage,  and 
the  amount  of  rapidly  solidified  material,  as  determined  from  the 
growth  rate,  is  diminished.  However,  the  "recalesced"  zone  is 
actually  a  smaller  fraction  of  the  overall  volume,  and  therefore 
the  average  droplet  undercooling  is  likely  to  be  higher.  It  is 
then  concluded  that  the  Ijj-e  path  will  be  truly  indicative  of  the 
extent  of  rapid  solidification  only  when  the  droplet  satisfies 
the  Newtonian  condition  of  negligible  temperature  differences 
during  recalescence. 

Figure  23  also  shows  that  as  rQ  decreases  the  curves 
approach  a  limit  slightly  different  from  the  Newtonian  case, 
probably  due  to  the  differences  in  geometry  and  the  "reduced" 
kinetics  discussed  in  the  previous  section.  For  example, 
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TEMPERATURE  (K) 
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DIMENSIONLESS  TEMPERATURE  0  =  Cl(T-Tm)AHm 


Figure  23.  Departure  from  the  Newtonian  thermal  history  as 

particle  size  increases.  Note  that  the  enthalpy  and 
temperature  are  now  averaged  over  the  entire  droplet. 
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characteristic  undercoolings  AT(ga)  of  35.7,  22.2  and  21.8  K  were 
calculated  from  this  figure  for  100,  10  and  1  um  A1  droplets 
respectively,  whereas  the  corresponding  Newtonian  value  is  18  K. 

Figure  24  shows  the  temperature  profiles  along  the  growth 
axis  (i  =  n),  as  a  function  of  fraction  sol  <  'or  a  fixed 
diameter  of  10  um.  The  corresponding  drop^. .  temperatures  and 
interface  positions  calculated  from  the  Newtonian  model  are  given 
for  comparison.  Note  that  the  internal  temperature  differences 
A (T)  go  through  a  maximum  during  recalescence.  For  example, 
maximum  A(T)  values  of  7,  52  and  148  K  occuring  at  fractions 
solid  of  0.08,  0.108  and  0.126  were  calculated  for  A1  droplets  of 
1,  10  and  100  pm  in  diameter,  respectively,  under  the  conditions 
of  this  figure. 

Figure  25  shows  the  effect  of  particle  size  on  the 
temperature  profiles  for  a  constant  fraction  solid  g  =  0.1.  It 
is  evident  from  this  figure  that  decreasing  the  particle  size 
drastically  reduces  A(T),  the  profile  being  almost  flat  for 
droplets  of  about  1  um  in  diameter  or  less.  Further  analysis  of 
Figure  25  reveals  that  as  the  droplet  size  increases,  the  liquid 
away  from  the  nucleation  point  becomes  less  sensitive  to  the 
recalescing  interface,  and  may  even  cool  slightly  below  the 
original  0^.  Since  the  time  scale  of  the  solidification  process 
is  roughly  proportional  to  rQ,  it  can  be  concluded  that  the 
larger  particles  can  maintain  liquid  regions  undercooled  for 
longer  times,  and  therefore  are  more  likely  to  develop  multiple 
nucleation  events  and  polycrystalline  structures. 
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Figure  24.  Thermal  profiles  developed  inside  the  droplet  as 

solidification  progresses,  for  a  constant  particle 
diameter  of  10  urn. 
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Figure  25.  Thermal  profiles  developed  inside  the  droplet  as  a 
function  of  particle  size  for  a  constant  fraction 
solid  g=0.1. 
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Although  the  temperature  differences  are  reduced  inside  the 
droplet  as  the  radius  decreases,  the  average  temperature  gradient 
at  the  interface, 


GSL  =  2  (kSGS  +  kLGL^^S  + 


(2.7) 


for  a  given  fraction  solid  does  not  change  significantly  in 
actual  magnitude,  as  shown  in  Table  12.  As  a  matter  of  fact,  the 
temperature  gradients  at  the  interface  should  increase  in  the 
smaller  particles  in  order  to  sustain  the  larger  velocities 
expected  in  them,  at  least  in  the  early  stages  of  recalescence. 


Table  12 

Average  Gradients  at  the  Interface  During 


-lv  + 


Solidification  of  Undercooled  Aluminum  Droplets  (K  pm  ) 


Fraction 
solid,  g 

1 

Diameter  (pm) 

10 

100 

0.1 

-15.3 

-12.4 

-4.7 

0.2 

-11.6 

-11.0 

-4.5 

0.3 

-  6.7 

-  6.7 

-4.3 

0.4 

-  2.5 

-  3.6 

-4.0 

0.5 

-  0.7 

-  1.2 

-1.9 

182  K,  h  = 

5xl04  Wm' 

'2K-1,  Km  =  2  cms 

-V1 

N 


The  average  gradient  GgL  is  always  negative  during 
recalescence  because  it  is  dominated  by  the  largely  negative 
Gl  term  in  Equation  (2.7).  It  then  follows  from  the  discussion 
on  morphological  stability  in  section  2.2.2  that  the  thermal 
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field  promotes  interface  breakdown  daring  solidification  of 
undercooled  droplets.  As  the  liquid  undercooling  is  exhausted 
Gl  and  Gsl  increase  towards  zero,  as  indicated  in  Table  12,  and 
the  latter  may  even  become  slightly  positive  near  the  end  of 
solidification,  therby  reducing  its  destabilizing  influence  on 
the  interface.  A  note  of  caution  is  necessary  with  regard  to  the 
gradient  values  calculated  from  this  program,  since  we  have 
introduced  a  localized  interface  in  a  volume  element  which  is 
assumed  to  have  uniform  properties.  Therefore,  the  values  will 
only  be  conceptually  correct  when  the  interface  position 
coincides  with  the  node  (v^  *  v?.jj)  and  its  temperature  is  equal 
to  that  of  the  node.  It  should  be  noted,  however,  that  the 
variations  in  the  gradients  as  the  interface  moves  inside  a  node 
are  within  an  order  of  magnitude. 

Figure  26  shows  the  effect  of  particle  size  on  interfacial 
temperatures  (hence  velocities)  at  the  growth  axis,  as  a  function 
of  fraction  solidified.  As  the  radius  increases,  the  interface 
recalesces  at  an  earlier  stage;  its  undercooling  and  velocity  are 
reduced  for  an  equivalent  g.  In  addition,  the  curves  tend  to 
develop  a  plateau  (e.g.  2r^  =  100  um)  of  relatively  constant 

velocity  as  they  approach  the  limiting  case  of  a  solid  growing 
into  an  infinite  liquid  bath.  On  the  other  hand,  as  rQ 
decreases,  the  interface  undercooling  should  approach  the 
Newtonian  behavior. 

The  extent  of  rapid  solidification  can  now  be  determined  by 
the  fraction  solidified  before  the  interface  reaches  a  critical 
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undercooling  ATC  necessary  to  produce  the  desired  modification  in 
the  microstr ucture •  The  metastable  effects  will  only  be  present 
throughout  the  particle  if  the  minimum  undercooling 

★ 

characteristic  of  the  post-recalescence  stage  (i.e.  iTmin  i-n 
Figure  26)  is  always  larger  than  AT^ .  This  minimum  undercooling 
is  fairly  independent  of  particle  size  and  close  to  the  Newtonian 
value  for  rQ  <  100  pm. 

Figure  27  compares  the  position,  shape  and  velocity  of  the 
liquid  solid  interface  as  solidification  progresses,  for  1,10  and 
100  urn  A1  droplets  and  the  corresponding  Newtonian  limit.  As 
expected,  the  growth  velocities  increase  as  the  radius  decreases, 
and  the  values  for  a  1  pm  particle  are  almost  identical  to  those 
calculated  from  the  Newtonian  model  (the  interface  and  droplet 
undercoolings  have  similar  values).  Note,  however,  that  the  Biot 
number  at  which  this  agreement  occurs  is  2.8  x  10-4,  far  below 
the  conventional  criterion  for  the  achievement  of  Newtonian 
conditions  (Bi  <  0.1).  Analysis  of  this  figure  also  reveals  that 
the  numerically  calculated  interface  profiles  are  similar  to  the 
spherical  configuration  assumed  in  the  Newtonian  model. 

The  interfacial  shape  is  basically  dictated  by  the  velocity 
distribution  along  the  solid  front.  If  the  undercoolings  are 
fairly  uniform  the  interface  is  expected  to  be  convex  towards  the 
liquid  because  the  distances  traveled  from  the  nucleation  point 
are  basically  constant  for  each  time  increment.  However,  if 
substantial  temperature  gradients  are  present  in  the  droplet 
prior  to  nucleation,  (e.g.  large  rQ  and  h  values),  and  the 
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Km  =  2  cm  s' 1  K" 1 

eN  =  -0.5 

K/Bi  =  400 

h  =  5x10“  W  m’2  K'1 

Figure  27.  Interface  profiles  and  velocities  predicted  from  the 
Enthalpy  model  for  undercooled  solidification.  The 
deviations  from  the  Newtonian  velocity  limit  increase 
with  particle  size,  but  the  profiles  for  the 
recalescence  stage  are  quite  close  in  both  models  for 
rQ  <  100  pm. 
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initial  undercooling  is  small,  the  solid  growing  along  the 
symmetry  axis  faces  a  superheated  liquid  whereas  that  growing 
along  the  surface  faces  a  supercooled  liquid  and  also  "feels" 
more  the  external  cooling.  It  is  then  expected  that  the 
interface  will  become  concave  and  eventually  approach  the  case  of 
no-undercooling  analyzed  in  Chapter  7.  Figure  28  shows  the 
transition  from  a  convex  to  a  concave  solid  front  as  ti  is 
reduced  and  rg  is  increased  above  100  pm.  In  addition, 
increasing  KM  will  further  magnify  the  concavity.  In  the  limit, 
then,  the  droplet  should  get  covered  very  rapidly  with  a  solid 
"shell"  after  nucleation,  and  growth  would  proceed  in  the 
concentric  pattern  discussed  before. 

In  closing  this  section,  it  should  be  noted  that  the  same 
effects  obtained  by  changing  the  particle  radius  are  produced  if 
the  heat  transfer  and  kinetic  coefficients  are  proportionally 
increased  or  decreased,  keeping  the  K/Bi  ratio  constant. 
Furthermore,  the  Biot  number  at  which  the  numerical  and  Newtonian 
velocities  converge  will  be  dependent  on  the  kinetic  coefficient 
and  initial  undercooling. 

8.2.2  Effects  of  the  heat  transfer  coefficient 

Increasing  the  heat  transfer  coefficient  at  constant  and 
rQ  enhances  the  role  of  the  external  cooling  during  recalescence , 
and  displaces  the  ijj-0  curve  following  the  trend  established  in 
the  Newtonian  analysis.  As  previously  indicated,  the  two  models 
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diverge  in  their  prediction  of  the  thermal  history,  but  the 
magnitude  of  this  deviation  becomes  relatively  insensitive  to  h 
as  Tq  decreases  below  100  pm.  It  was  also  noted  that  changes  in 
the  heat  transfer  coefficient  do  not  alter  significantly  the 
maximum  temperature  differences  A(T)  developed  inside  the 
droplet.  For  example,  the  maximum  A (T)  in  a  10  ym  A1  droplet 
with  an  initial  undercooling  of  182  K  (0N  =  -0.5)  and  a  kinetic 
coefficient  KM  =  2  cm  s~^  K-1  decreases  from  54  to  47  as  h  is 
reduced  from  2  x  10^  to  2  x  102  W  m-2  K-^  (K/Bi  =  102  and  10^ 
respectively) . 


Table  13 


Effect 

of  the  Heat 

Transfer  Coefficient 

on  the 

Interface 

Velocities  (cm  s 

-1)  + 

Fraction 

h  (kW  m-2  K-1) 

solid,  g 

200 

20 

2 

0.1 

266 

236 

233 

(323) 

(292) 

(289) 

0.3 

192 

135 

126 

(204) 

(149) 

(142) 

0.5 

111 

21 

3.1 

(102) 

(18) 

(3.1) 

Minimum 

62 

5.4 

0.5 

velocity 

(46) 

(4.3) 

(0.4) 

t  2rQ  =  10  ym,  ATN  =  182  K,  KM  =  2  cms_1K  1 

Values  in  parenthesis  are  the  Newtonian  results. 


Figure  29  illustrates  the  effect  of  h  on  interface 

undercooling  as  a  function  of  fraction  solid,  whereas  Table  13 

gives  examples  of  the  corresponding  interface  velocities.  As 

expected  from  the  Newtonian  analysis,  the  role  of  the  heat 

transfer  coefficient  is  larger  as  solidification  progresses, 

because  the  growth  rate  (and  hence  the  release  of  latent  heat)  is 

slowed  down  as  the  interface  recalesces.  In  general,  the 

fraction  of  material  solidified  above  a  certain  velocity 

increases  with  h,  but  the  effect  is  only  important  beyond  a 

relatively  large  value  (e.g.  2  x  10^  W  m-^  K-1  in  this  figure). 

Furthermore,  the  heat  transfer  coefficient  is  more  influential  on 

* 

the  minimum  undercooling  ATmin,  which  increases  in  a  direct 
proportion  to  h.  It  is  expected  that  the  importance  of  h  during 
recalescence  will  be  magnified  as  the  droplet  size  increases  and 
the  kinetic  coefficient  or  initial  undercooling  are  diminished. 

8.2.3  Effects  of  the  kinetic  coefficient 

Increasing  the  kinetic  coefficient  K„  at  constant  h  and 

M 

r q  reduces  the  role  of  the  external  cooling  on  the  thermal 
history,  in  agreement  with  the  trend  previously  indicated  by  the 
Newtonian  analysis.  The  differences  between  the  Ijj-e  curves 
predicted  by  the  two  models  are  small  for  the  range  of  values  of 
interest.  However,  it  was  found  that  changes  in  KM  significantly 
alter  the  temperature  profiles  and  maximum  a(T)  developed  during 
recalescence,  as  shown  in  Figure  30.  For  example,  increasing  the 
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FRACTION  SOLID  g 


Figure  29.  Dependence  of  the  interface  temperature  on  the  heat 
transfer  coefficient  during  solidification.  The 
curves  are  insensitive  to  h  below  2  kW  m~2  k-1  for 
these  particular  values  of  rQ,  0^  and  KM.  Note  that  h 
becomes  more  relevant  in  the  post  recalescence  stage. 
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DIMENSIONLESS  DISTANCE  ALONG  GROWTH  AXIS 


Figure  30.  Temperature  profiles  as  a  function  of  the  kinetic 
coefficient  and  particle  radius  for  a  constant 
fraction  solid  g=0.1. 
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kinetic  coefficient  from  2  to  50  cm  s--*-  K-1  increases  the  maximum 
temperature  difference  inside  a  1  ym  droplet  from  7  to  82  K,  and 
in  a  10  iim  droplet  from  52  to  154  K,  for  an  initial  undercooling 
of  182  K  (0N  =  -0.5)  . 

The  influence  of  KM  on  the  interfacial  temperature  during 
solidification  is  shown  in  Figure  31,  and  the  corresponding 
velocities  are  given  in  Table  14. 


Table  14 

Effect  of  the  Kinetic  Coefficient 
on  the  Interface  Velocities  (cm  s-1) ^ 


Fraction 
solid,  g 

2 

km  (cm  s'1 
8 

K-l) 

50 

0.1 

240 

640 

1456 

(294) 

(1163) 

(7221) 

0.3 

146 

405 

1074 

(159) 

(548) 

(3555) 

0.5 

41 

64 

127 

(37) 

(51) 

(79) 

Minimum 

14 

14 

12 

velocity 

(ID 

(ID 

(11) 

.4  __ 

+  2r Q  =  10  um,  ATN  =  182  K,  h  =  5xl04  Wm'2K_1 


Values  in  parenthesis  are  the  Newtonian  results. 

For  a  fixed  ATN  the  initial  velocity  is  purely  controlled  by  the 
growth  kinetics  (i.e.  KM  and  kinetic  model),  and  therefore 
independent  of  the  heat  flow.  On  the  other  hand,  the  minimum 
velocity  characteristic  of  the  post  recalescence  stage  is  nearly 
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Figure  31.  Effect  of  the  kinetic  parameter  on  the  interface 
undercooling  during  solidification.  Although 
decreases  AT*  at  a  particular  fraction  solid,  tne 
actual  velocity  increases. 
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independent  of  the  KM  value  and  close  to  the  Newtonian 
prediction,  meaning  that  the  interface  adjusts  its  undercooling 
to  satisfy  the  conditions  prescribed  by  the  surface  cooling. 
During  recalescence ,  the  interface  goes  through  a  transition 
between  the  above  situations.  It  alters  both  its  undercooling 
and  the  thermal  field  in  its  vicinity  (i.e  GL  and  Gs)  in  order  to 
accomodate  the  latent  heat  released  by  the  growth  process. 

Analysis  of  Table  14  reveals  that  the  interface  velocities 
calculated  from  the  Enthalpy  model  are  markedly  less  sensitive  to 
the  value  of  the  kinetic  coefficient  than  those  obtained  from  the 
Newtonian  simulation.  This  result  arises  from  the  limitations  to 
heat  flow  inside  the  droplet  (finite  thermal  conductance)  which 
lead  to  the  "local"  recalescence  of  the  interface.  As  KM 
increases  the  system  tends  to  develop  a  regime  of  steady  state 
growth,  resembling  the  case  of  solidification  into  an  infinite 
liquid  pool.  It  is  then  concluded  that  increasing  KM  decreases 
the  Biot  number  at  which  the  Newtonian  model  provides  a  good 
estimate  of  the  interface  velocities.  Furthermore,  it  is 
expected  that  the  effects  of  will  be  enhanced  by  reducing  the 
particle  size  or  increasing  the  initial  undercooling  and  heat 
transfer  coefficient. 

8.2.4  Effects  of  the  nucleation  temperature 

Increasing  the  undercooling  level  prior  to  nucleation  is 
perhaps  the  most  efficient  way  of  improving  the  extent  of  rapid 
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solidification,  as  concluded  from  the  Newtonian  analysis.  The 
. curve  calculated  from  the  numerical  model  progressively  shifts 
towards  the  corresponding  Newtonian  curves  as  /.T  decreases. 
Furthermore,  increasing  ATN  magnifies  the  temperature  gradients 
inside  the  droplet  at  a  given  fraction  solid  during  recalescence , 
as  shown  in  Figure  32.  Maximum  temperature  differences  of  25  and 
94  K  were  calculated  for  10  tun  A1  droplets  with  initial 
undercoolings  of  191  and  364  K  respectively,  given  the  and  h 
values  of  this  figure. 


Table  15 

Effect  of  the  Nucleation  Temperature 
on  the  Interface  Velocities  (cm  s~l) + 


Fraction 

Average 

Nucleation 

Temperature 

solid,  g 

-0.25 

-0.50 

-1.0 

0.01 

154 

306 

634 

(189) 

(371) 

(742) 

0.25 

45 

169 

466 

(45) 

(198) 

(553) 

0.50 

13 

41 

341 

(11) 

(40) 

(364) 

1.00 

95 

95 

110 

(36) 

(36) 

(54) 

0  =  10  um. 

h  =  5xl04 

Wm~2K_1 ,  K, 

J 

u  =  2  cms 

Values  in  parenthesis  are  the  Newtonian  results. 


Figure  33  shows  the  interfacial  undercoolings  during 
solidification  for  several  0N  values  and  two  kinetic  models. 
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DIMENSIONLESS  DISTANCE  ALONG  GROWTH  AXIS 


Figure  32.  Temperature  profiles  as  a  function  of  initial 
undercooling  for  a  fraction  solid  g=0.1. 
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Examples  of  the  corresponding  interface  velocities  are  given  in 
Table  15.  It  is  evident  from  this  figure  that  the  amount  of 
material  solidified  above  a  critical  undercooling  ATC  markedly 
increases  as  nucleation  is  suppressed  to  larger  undercoolings. 

The  local  recalescence  is  also  faster  as  0N  decreases  due  to  the 
consequent  increase  in  the  driving  force  for  solidification. 
However,  the  assumption  of  an  exponential  kinetic  law  reduces 
this  effect  because  it  consistently  predicts  lower  velocities 
than  the  linear  law  for  a  given  temperature.  As  in  the  Newtonian 
case,  the  kinetic  models  tend  to  produce  similar  results  as  the 
undercooling  is  reduced.  Finally,  it  is  expected  that  reducing 

km  and  rg,  as  well  as  increasing  h  will  enhance  the  role  of  the 
initial  undercooling  on  the  recalescence  behavior  of  the 
droplets. 
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CHAPTER  9 


MICROSTRUCTURAL  ANALYSIS 


Aluminum  alloy  powders  were  produced  in  external  facilities 
from  ingots  made  in  our  laboratories,  and  their  microstructures 
examined  by  Transmission  (TEM)  and  Scanning  Electron  Microscopy 
(SEM) .  Determination  of  lattice  parameters  was  also  performed  on 
some  samples  by  X-ray  dif f ractometry . 

Most  of  the  studies  were  focused  on  the  Al-Si  system,  using 
alloys  of  3  and  6%  Si,  as  well  as  a  high  purity  (99.999%) 
aluminum  standard.  Some  supplementary  observations  were  made  on 
Al-4 . 5%Cu .  The  alloys  were  prepared  from  A1  ingot  (99.97%),  Cu 
plate  (99.95%)  and  Si  lumps  (99.98%)  by  melting  batches  of  20  kg 
in  a  clay-graphite  crucible,  degassing  with  chlorine  and  casting 
in  sand  or  graphite  molds,  depending  on  the  projected  atomization 
technique . 

9.1  Atomization 

Powder  samples  were  produced  by  two  different  methods, 
namely,  electrohydrodynamic  atomization  in  vacuum  (Phrasor 
Scientific  EHD  process)  and  high  speed  centrifugal  atomization  in 
helium  (Pratt  and  Whitney  RSR  process) .  These  techniques  have 
been  extensively  discussed  elsewhere  [77,90,91]  and  are  briefly 
described  in  what  follows. 
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9.1.1  Electrohydrodynamic  atomization  (EHD) 

The  EHD  method  for  droplet  generation  involves  the 
application  of  very  intense  electric  fields  on  a  conductive 
liquid  at  the  tip  of  a  capillary  nozzle  (emitter) .  The  molten 
metal,  contained  in  a  pressurized  ceramic  reservoir,  is  given  a 
positive  potential  (usually  +10kV)  while  delivered  to  the  nozzle, 
typically  of  50  um  in  diameter.  The  emmiter  is  exposed  to  a 
negative  (-1  kV)  annular  electrode  (extractor) ,  concentrating  the 
electric  field  on  the  very  small  liquid  meniscus.  The 
interaction  of  the  electrostatic  stresses  with  the  surface 
tension  forces  produces  a  liquid  spike  that  breaks  into  charged 
droplets.  These  are  then  accelerated  by  the  voltage  difference 
and  fly  through  the  orifice  in  the  extractor  to  be  collected  on  a 
target  after  solidification.  Velocities  above  200  m  s-1  have 
been  reported  for  particles  below  1  |im  in  diameter,  with 
calculated  times  of  flight  lower  than  2.5  ms  [91]. 

The  atomizer  is  enclosed  in  a  vacuum  chamber  (lO--*  torr)  , 
thereby  restricting  the  heat  loss  from  the  droplets  to  radiation 

—  9—1 

(h  <  100  W  m  K  ) .  Maximum  cooling  rates  prior  to 
solidification  are  estimated  to  range  from  9.6  x  104  K  s-1  at  the 
melting  temperature,  to  1.3  x  104  K  s-1  at  the  onset  of 
hypercooling,  for  a  1  um  droplet  of  aluminum  with  an  emissivity 
e  ~  1.  Actual  values  may  be  lower  by  more  than  one  order  of 
magnitude  if  we  consider  the  e  values  reported  in  the  literature 
and  mentioned  in  Chapter  3.  In  addition,  the  cooling  rates  for  a 


100  nm  particle  will  be  two  orders  of  magnitude  lower,  as 
indicated  by  Equation  (6.4). 

EHD  powders  range  in  size  from  less  than  100  nm  to  about 
100  nm.  They  were  collected  by  deposition  on  a  125  pm  thick 
plastic  tape  (of  the  type  used  to  prepare  electron  microscopy 
replicas)  placed  on  the  target.  By  appropriate  manipulation  of 
the  process  variables,  the  substrate  may  be  covered  with  a 
monolayer  consisting  mostly  of  submicron  powders.  Particles 
larger  than  ~10  nm  are  usually  limited  in  number,  presumably 
because  they  bounce  off  or  are  later  detached  from  the  tape. 

This  loose  material  is  gathered  from  the  chamber  at  the  end  of 
the  experiment. 

Figure  34  shows  a  typical  sample  of  powders  as  received  on 
the  plastic  substrate.  They  are  nearly  spherical  and  have  smooth 
surfaces,  although  some  appear  to  have  solidified  on  the  tape. 
Mainly  intended  as  a  laboratory  tool,  this  process  is  able  to 
produce  only  a  few  grams  of  material  per  hour.  Nevertheless,  it 
has  proven  to  yield  a  substantial  number  of  samples  under 
reproducible  conditions,  making  it  specially  suitable  for 
fundamental  solidification  studies. 

9.1.2  Centrifugal  atomization  ( RSR) 

The  RSR  (for  Rapid  Solidification  Rate)  process  for  the 
production  of  metal  powders  is  based  on  a  turbine-driven  atomizer 
capable  of  reaching  velocities  in  excess  of  24000  rpm.  The 
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atomizer  consists  of  a  water  cooled  copper  disc,  7.5  cm  in 
diameter,  rotating  about  a  vertical  axis  and  enclosed  into  a 
helium  filled  chamber.  A  controlled  stream  of  liquid  metal  is 
poured  near  the  center  of  the  disk,  spreading  into  a  layer  due  to 
the  centrifugal  acceleration.  At  the  edge,  the  liquid  breaks 
down  into  droplets  with  tangential  velocities  of  the  order  of 
100  m  s”l.  These  are  then  quenched  by  a  large  cross-flow  of 
helium  injected  into  the  chamber,  and  the  solid  powders  are 
collected  from  the  gas  exhaust  line. 

Convection  is  then  the  dominant  mechanism  of  heat  extraction 
in  these  powders,  and  heat  transfer  coefficients  above 
10^  w  m-2  k-1  are  claimed  to  be  achievable  by  this  method. 

Cooling  rates  of  6.8  x  108  K  s-1  at  the  melting  temperature  and 
2.9  x  108  K  s"1  at  the  onset  of  hypercooling  were  calculated  for 

tr  _2  —1 

a  1  pm  droplet  using  the  h  =  4.8  x  10  W  m  K  given  in 
Table  2.  Note  in  this  table  that  increasing  the  particle  size  to 
100  pm  reduces  h  roughly  by  two  orders  of  magnitude,  hence 
decreasing  the  achievable  cooling  rates  by  four  orders  of 
magnitude,  as  indicated  by  Equation  (6.4). 

RSR  powders,  shown  in  Figure  35,  are  fairly  spherical  in 
shape  but  with  rough  surfaces  (as  if  they  had  undergone 
collisions  during  or  after  solidification) .  They  apparently  lack 
the  "orange  peel"  topography  typical  of  materials  produced  by 
atomization  [92]  and  range  in  size  from  slightly  below  1  pm  to 
above  100  pm.  Standard  size  classification  of  the  Al-Si  powders 
reveals  that  about  80%  of  the  material  received  is  smaller  than 


Figure  35.  SEM  view  of  Al-3%Si  powders  produced  by  centrifugal 
atomization;  (a)  particle  size  below  mesh  325 
(<44  i  m) ,  (b)  particle  size  between  100  and  120  mesh 


(125-149  tun)  . 


141 


125  pm,  and  15%  passes  mesh  325  (<44  , m) .  However,  SEM  analysis 
revealed  that  the  coarse  fraction  of  the  powders  (+100  mesh)  tend 
to  be  formed  by  conglomerates  rather  than  individual  particles. 
Yields  of  1-2  kg  are  common  for  a  batch  aluminum  run,  but  the 
process  is  aimed  to  produce  0.5  kg  s-1  of  material  in  a 
continuous  fashion  at  commercial  scale  [77] . 

9.2  Specimen  Preparation 

Samples  of  pure  aluminum  and  the  three  alloys  were  obtained 
from  EHD  atomization,  but  only  the  alloys  were  processed  through 
the  RSR  technique.  EHD  powders  were  received  on  sections  of  the 
tape  located  on  the  target  during  atomization.  Small  quantities 
of  loose  particulate  material  containing  both  powders  and  splats 
were  also  provided  for  most  runs.  RSR  powders  were  furnished  in 
bulk  form. 

Aluminum  powders  below  1  pm  in  diameter  are  transparent  to 

200  keV  electrons,  and  therefore  lend  themselves  to  conventional 

TEM  analysis  without  further  thinning.  EHD  specimens  were 

prepared  by  cutting  the  tape  into  3  mm  discs  and  coating  on  the 

deposition  side  with  a  20-50  nm  film  of  evaporated  carbon.  The 

substrate  was  then  dissolved  in  acetone  and  the  powders 

attached  to  the  carbon  film.  This  was  carefully  rinsed 

isopropyl  alcohol  and  placed  on  a  3  mm,  400  mesh  copper 

The  powder  density  was  commonly  about  100  particles  per 
3  2 

opening  (~10  um  )  . 
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Preparation  of  TEM  specimens  of  submicron  RSR  powders  was 
achieved  by  separating  the  finest  fraction  (-325  mesh)  from  each 
alloy,  suspending  a  small  amount  of  this  material  in  methanol, 
and  allowing  it  to  settle  with  time.  Samples  were  taken  with  a 
thin  pipet  at  regular  intervals  and  a  drop  of  each  one  was  placed 
onto  a  3  mm  disc  of  blank  replicating  tape  (previously  treated  to 
obtain  a  smooth  surface) .  The  operation  was  repeated  a  few  times 
for  each  specimen.  Once  the  methanol  had  evaporated,  the  discs 
were  coated  with  carbon  and  the  films  recovered  in  the  same 
manner  described  for  the  EHD  powders.  It  should  be  noted  that 
the  density  of  electron-transparent  material  in  these  latter 
specimens  was  only  of  a  few  particles  per  disc. 

Loose  powders  from  both  processes  were  mounted  on  resin, 
polished  with  diamond  following  standard  metallographic  practices 
and  examined  by  SEM  (25  keV) .  Since  the  atomic  number  contrast 
in  the  Al-Si  alloys  is  rather  poor,  the  as-polished  sections  were 
etched  in  1%  NaOH  for  90  seconds  to  enhance  the  surface 
topography. 

9.3  Results  and  Discussion 

The  powders  were  divided  for  their  study  in  submicron  and 
conventional  (>10  ym)  material.  Although  the  former  are  not  of 
commercial  importance,  their  electron- transparency  permits  the 
three-dimensional  examination  of  a  complete  solidification 
structure.  Furthermore,  our  analysis  revealed  that  substantial 


undercolings  were  achieved  in  most  submicron  droplets  prior  to 
nucleation,  and  the  resulting  microstructures  were  in  agreement 
with  the  trends  predicted  by  the  heat  flow  model. 

9.3.1  EHD  submicron  powders 

The  observations  resulting  from  the  analysis  of  EHD  powders 
will  be  presented  depending  on  whether  they  pertain  to  their 
segregation  or  their  nucleation  behavior.  For  the  purposes  of 
the  ensuing  discussion,  we  will  call  segregates  those  second 
phases  originating  from  the  melt,  whereas  precipitates  will  refer 
to  phases  formed  after  solidification. 

Most  of  the  Al-Si  EHD  powders  conform  to  a  variant  of  the 
basic  segregation  pattern  shown  in  Figure  36.  The  microstructure 
exhibits  a  region  (A)  containing  a  uniform  distribution  of  very 
fine  precipitates  or  clusters  of  Si  which,  according  to  the 
diffraction  evidence,  have  an  orientation  relationship  to  the 
primary  A1  phase.  Dark  field  imaging  of  precipitates  slightly 
coarser  than  those  in  the  present  figure  suggest  that  they  are 
either  needles  or  platelets  on  the  {ill}  A1  planes.  This  finding 
is  consistent  with  the  decomposition  kinetics  of  Al-Si  solid 
solutions;  the  plates  evolve  from  spherical  clusters  and  needles 
that  form  in  the  early  stages  of  precipitation  [31-36] .  It  was 
then  concluded  that  this  zone  was  originally  a  supersaturated 


solid  solution. 


Typical  segregation  pattern  for  an  Al-6%Si  EHD 
submicron  particle  (2rQ>  0.5  urn)  showing  3  distinct 
regions.  The  initially  rapid  solidification  produces 
a  supersaturated  zone  (A),  which  breaks  later  into 
cells  (C)  after  the  droplet  recalesces.  (B)  is  a 
transition  zone  where  solute  builds  up  in  front  of  th 
interface.  The  growth  axis  is  close  to  a  <111> 
direction . 


Figure 
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Region  (A)  is  followed  by  a  precipitate-free  zone  (B)  ,  whose 
configuration  closely  resembles  the  axisymmetr ic  geometry 
discussed  in  the  heat  flow  analysis  and  depicted  in  Figure  27. 

The  boundary  between  regions  (A)  and  (B)  corresponds  to  the  trace 
of  the  solidification  front,  which  is  symmetric  about  the  growth 
axis.  The  latter  was  consistently  found  to  be  close  to  a  <111> 
direction. 

Finally,  the  particle  in  Figure  36  also  has  a  well  developed 
cellular  region  (C) ,  with  segregate  spacings  in  the  order  of 
0.1  ym.  Similar  structures  in  Al-4.5%Cu,  shown  at  a  different 
angle  in  Figure  37,  reveal  more  clearly  the  cellular  pattern. 

Figure  38  illustrates  the  effect  of  particle  size  on  the 
extent  of  the  different  zones  in  the  segregation  pattern.  As  the 
diameter  decreases,  the  cellular  structure  becomes  less  defined 

(b)  ,  and  eventually  disappears,  leaving  only  a  few  segregates 

(c)  .  At  the  lower  end  of  the  size  distribution,  the  particles 
appear  almost  featureless  (d) ,  showing  only  the  supersaturated 
region  (A) .  The  intensity  of  the  silicon  spots  is  also  reduced 
with  diminishing  diameter,  although  the  crystallographic 
relationship  remains. 

The  structures  in  Figure  38  are  by  no  means  restricted  to  a 
particular  size  range.  For  example,  segregate-f ree  powders  as 
large  as  1  ym  were  found  in  some  instances,  and  segregated  Si  was 
detected  in  some  of  the  smallest  particles.  The  trend 
illustrated  in  this  figure  was  established  from  the  TEM  analysis 
of  over  400  particles  between  0.1  and  1.0  ym  for  each  Al-Si 


Figure  38.  Transition  from  a  partially  segregated  structure  to  a 
homogeneous  one  as  particle  size  decreases.  The 
diameters  of  the  powders  shown  are  (a)  570  nm, 

(b)  470  nm,  (c)  220  nm  and  (d)  170  nm. 


PERCENTAGE  WITHIN  SIZE  RANGE 


PARTICLE  DIAMETER  (yrn) 


Figure  39.  Segregation  behavior  of  aluminum-silicon  submicron 
powders  produced  by  EHD.  Results  are  from  over  400 
particles  in  each  alloy  composition. 
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composition.  The  results,  given  in  Figure  39,  confirm  that 
increasing  the  particle  size  or  the  solute  content  enhance  the 
incidence  of  segregation.  Since  not  all  of  the  segregated 
powders  exhibit  a  cellular  pattern  (see,  for  example.  Figures 
38(c)  and  40),  the  occurrence  of  the  latter  has  been  also  plotted 
in  Figure  39. 

The  previous  microstructures  reveal  that  substantial 
undercoolings  must  have  been  reached  in  the  liquid  droplets  prior 
to  solidification,  which  is  believed  to  proceed  as  follows.  In 
general,  nucleation  occurs  at  the  surface,  as  evidenced  by  the 
growth  direction.  The  subsequent  rapid  solidification  into  a 
highly  undercooled  liquid  produces  a  homogeneous  zone  (A)  of 
extended  solubility.  As  the  droplet  recalesces,  the  interface 
velocity  and  undercooling  are  reduced,  thereby  decreasing  the 
driving  force  for  solute  "trapping".  When  the  interface 
temperature  surpasses  some  critical  temperature  for  the 
supersaturation  process,  Tc,  solute  is  rejected  ahead  of  the 
solidification  front.  The  interface  is  eventually  rendered 
unstable  and  breaks  into  cells. 

The  extent  of  homogeneous  solidification  can  be  rationalized 
in  terms  of  the  heat  flow  concepts  previously  developed. 

Assuming  that  the  thermophysical  properties  are  reasonably 
constant,  the  behavior  of  an  alloy  growing  in  a  plane  front  mode 
without  solute  partitioning  is  expected  to  resemble  that  of  the 
pure  metal. 


t 


Equation  (2.11)  for  constant  specific  heat  indicates  that 
the  required  nucleation  temperature  for  adiabatic  homogeneous 
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solidification  is 


TN  1  TH  -  TC  -  AVCL  (9-D 

where  T^  is  now  the  hypercooling  temperature  for  an  alloy. 

As  indicated  in  Chapter  2,  the  critical  temperature  for 
homogeneous  solidification  has  a  value  between  the  solidus  Tg  and 
the  Tq  curve.  For  purposes  of  illustration,  we  have  arbitrarily 
assigned  to  Tq  the  average  value  of  TL  and  Ts,  and  estimated  on 
this  basis  an  upper  bound  to  the  hypercooling  temperature  TH*  see 
Table  16.  The  lower  bound  to  TR  corresponds  to  Tc  =  Tg.  The 
values  of  0  were  calculated  by  substituting  T^  for  tm  in 
Equation  (3.12)  and  hence  are  proportional  to  the  undercoolings. 


Table  16 

Effect  of  the  Alloy  Composition  on  the 
Critical  and  Hypercooling  Temperatures^ 


3%  Si  6%  Si 


T 

e 

T 

e 

Liquidus 

tl 

913 

0 

893 

0 

Critical 

TC<To> 

848 

-0.18 

762 

-0.36 

Tc<Ts> 

782 

-0.36 

631 

-0.72 

Hyper coo ling 

TH(To> 

484 

-1.18 

398 

-1.36 

W 

418 

-1.36 

267 

-1.72 

t  Estimated 

from  the 

equilibrium  phase  diagram,  Te 

are 

from  the  metastable 

solidus . 

To  -  ITL 

+  Tg)/2 
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The  values  in  this  table  indicate  that  dimensionless 
temperatures  in  the  order  of  -1.18  (ATn/tl  =  0.47)  and  -1.36 
(ATn/tl  =  0.55)  may  be  necessary  to  eliminate  segregation  in 
aluminum  alloys  with  3  and  6%  Si  respectively.  Note  that  the 
ratios  of  ATN/TL  are  much  higher  than  the  maximum  of  about  0.35 
that  has  been  obtained  in  the  controlled  undercooling  experiments 
with  emulsified  metals  [49-51] . 


Table  17 

Extent  of  Rapid  Solidification  in  Al-Si  Alloys  ^ 

as  Fraction  Solidified  Below  a  Critical  Temperature  g(6„) 


Critical 

Temperature 

eC 

(Wm  •*■) 

0N 

CO 

r— 1 
• 

O 

1 

-0.36 

-0.72 

20 

-1.0 

0.80 

0.62 

0.26 

(EHD) 

-0.75 

0.56 

0.38 

0.03 

-0.50 

0.31 

0.13 

- 

-0.25 

0.07 

- 

- 

2xl05 

-1.0 

0.92 

0.67 

0.28 

(RSR) 

-0.75 

0.67 

0.43 

0.05 

-0.50 

0.43 

0.20 

- 

-0.25 

0.21 

- 

- 

t  Newtonian  cooling,  KM  =  2  cms  ^ 


Table  17  shows  calculated  Newtonian  values  of  the  fraction 
solidified  below  different  critical  temperatures,  g(ec)  (i.e., 
the  extent  of  super saturation) ,  for  different  nucleation 
temperatures  and  heat  transfer  coefficients.  For  an 
h  =  20  W  m"2  K-^,  typical  of  EHD  atomization,  heat  extraction 


from  the  droplet  does  not  have  a  significant  role  during 
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recalescence,  as  previously  shown.  It  is  concluded  from  this 
table  that  the  rapid  solidification  effects  are  diminished  when 
the  critical  temperature  decreases  as  a  result  of  increasing  the 
solute  content.  In  addition,  the  achievable  undercooling  prior 
to  nucleation  is  likely  to  increase  as  the  droplet  size 
decreases,  and  the  initial  zone  of  super saturation  is 
consequently  extended  until  it  covers  the  entire  powder. 

It  should  be  noted  that  the  present  Newtonian  data  are  upper 
limits  to  g (ec)  and  independent  of  particle  size.  Increasing 
Tq  or  Km  will  promote  "local"  recalescence  of  the  interface  and 
therefore  reduce  g(6c)  for  a  given  h  and  0N.  Finally,  one  must 
also  conclude  that  undercoolings  near  or  into  the  hypercooled 
regime  were  obtained  for  those  particles  where  segregation  was 
absent. 

Figure  39  also  shows  the  existence  of  a  large  proportion  of 
powders  in  the  lower  sizes  that  experienced  segregation  but  no 
cellular  breakdown.  This  phenomenon  may  simply  reflect  that  the 
onset  of  segregation  (AT  <  AT^,)  occurred  very  late  in  the 
solidification  process.  Moreover,  interface  breakdown  may  also 
be  prevented  by  the  stabilization  of  short  wavelength 
perturbations  due  to  surface  tension  effects  [64] ;  the  smallest 
perturbation  that  can  grow  is  of  dimensions  comparable  to  the 
particle  size. 

A  minor  proportion  of  the  powders  (<5%)  exhibit  quite  a 
different  type  of  segregation  microstructure,  where  large  faceted 
Si  second  phases  appear  distributed  throughout  the  surface,  as  in 
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Figure  40. 
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morphology  of  Si  segregates  observed  in  some 
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Figure  40.  Diffraction  analysis  did  not  reveal  any  particular 
orientation  relationship  with  the  matrix,  but  rather  strong 
reflections  of  individual  silicon  crystals.  The  results  indicate 
that  these  phases  were  probably  segregation  products,  but  the 
growth  direction  and  geometry  were  not  evident,  and  the 
supersaturated  zone  (A)  typical  of  most  submicron  particles  was 
apparently  absent.  Furthermore,  this  anomalous  segregate 
morphology  was  not  observed  in  the  Al-Cu  system,  and  its 
occurrence  was  more  common  in  the  larger  particles  and  higher  Si 
alloys. 

It  was  concluded  that  powders  showing  the  faceted  silicon 
were  unlikely  to  have  reached  substantial  undercooling  prior  to 
nucleation,  hence  being  extensively  segregated.  Although  the 
absence  of  cells  could  be  ascribed  to  interface  stabilization  by 
surface  tension  effects,  it  is  not  clear  why  most  of  the 
particles  of  equivalent  sizes  were  cellular.  In  addition,  the 
distribution  of  segregates  was  not  consistent  with  an 
axisymmetric  growth  pattern,  in  which  the  Si  should  be 
accumulated  on  one  side  of  the  powder. 

High  purity  aluminum  submicron  powders,  such  as  those  in 
Figure  41,  do  not  exhibit  any  evidence  of  cellular  structures 
that  could  have  been  originated  from  thermal  instability  due  to 
the  negative  gradient  in  front  of  the  interface.  This 
observation  lends  support  to  the  fact  that  the  breakdown  observed 
in  the  alloys  was  mainly  of  constitutional  nature.  It  could  also 
be  argued  that  thermal  cells  in  a  pure  metal  would  be 
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Figure  41.  Typical  structures  in  high  purity  aluminum  (99.999%) 
EHD  submicron  powders.  (a)  shows  a  monocrystalline 
particle  and  (b)  a  set  of  twin  crystals. 
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indistinguishable  after  solidification,  unless  they  are  decorated 
by  impurities  or  strained  areas  caused  by  crystalline 
misor ientation.  Nevertheless,  pure  aluminum  particles  show  some 
platelike  features  that  can  be  interpreted  either  as  clusters  of 
impurities  or  vacancy  condensation  loops.  The  existence  of 
denuded  zones  near  the  surface  and  crystal  boundaries  suggests 
that  vacancies  are  involved  in  the  formation  of  these  features. 

The  microstructural  analysis  also  revealed  that  most  EHD 
powders  are  monocrystalline,  indicating  that  nucleation  was 
usually  limited  to  a  single  event  per  droplet.  A  substantial 
proportion  of  them,  however,  consist  of  two  or  more  crystals 
which  could  be  either  randomly  oriented  with  incoherent 
boundaries,  as  in  Figure  42,  or  twins,  as  in  Figure  43.  It  was 
concluded  that  each  rendom  grain  originates  from  an  independent 
nucleation  event,  and  its  boundaries  result  from  impingement  with 
other  growing  crystals.  On  the  other  hand,  twin  crystals 
normally  exhibit  {ill}  coherent  boundaries  and  are  believed 
unlikely  to  be  independently  nucleated. 

The  nucleation  behavior  of  Al-Si  EHD  submicron  powders  is 
quantitatively  summarized  in  Figure  44.  In  general,  the 
proportion  of  polycrystalline  material,  the  number  of  nucleation 
events  per  particle  and  the  incidence  of  twinning  all  increase  as 
the  powder  diameter  decreases  below  1  pm.  It  was  therefore 
concluded  that  larger  undercoolings  promote  the  occurrence  of 
twins  and  multiple  nucleation.  This  is  consistent  with  the 
qualitative  observation  that  polycrystalline  powders  are  commonly 


Figure  42, 


Typical  polycrystalline  powder  in  Al-3%Si  where  two 
nucleation  events  occurred  simultaneously  in  the 
droplet  and  formed  an  impingement  boundary. 
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Figure  44.  Nucleation  behavior  of  aluminum-silicon  powders 

produced  by  EHD  atomization.  Results  are  from  over 
400  particles  in  each  alloy  composition. 
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less  segregated  than  their  monocrystal  counterparts  of  equivalent 
size. 

The  incidence  of  multiple  nucleation  can  be  accounted  for  in 
terms  of  Equation  (2.2),  which  may  be  also  written  as  follows 

T 

N  =  V  /  N  [ J (T) / f J  dT  (9.2) 

T 

aM 

This  relationship  indicates  that  the  number  of  nuclei  depend  on 
the  volume  of  the  assembly  (droplet) ,  the  nucleation  frequency 
J (T)  expressed  now  as  a  function  of  temperature,  and  the  cooling 
rate  c.  Reducing  the  particle  size  will  decrease  the  volume 
(V  «  rQ)  and  increase  the  cooling  rate  (e  a  l/r0) ,  thereby 
reducing  N  roughly  with  a  4th  power  of  the  radius.  On  the  other 
hand,  J(T)  increases  with  higher  undercoolings  AT^,  and  therefore 
opposes  the  volume  effect  on  N.  Since  J(T)  behaves  approximately 
like  a  step  function  of  ATN,  small  undercoolings  (such  as  those 
expected  in  large  droplets)  might  not  affect  J  significantly  and 
the  number  of  grains  should  decrease  with  particle  size. 

However,  when  substantial  undercoolings  are  achieved  (as  in  the 
submicron  powders)  J(T)  may  increase  rapidly  with  ATN,  which 
varies  inversely  with  Cq.  There  is  conceivably  a  size  range  in 
which  the  increase  in  nucleation  frequency  overcomes  the  volume 
effect,  enhancing  the  probability  for  multiple  nucleation. 

The  occurrence  of  twins  is  a  phenomenon  more  difficult  to 
rationalize  than  multiple  nucleation.  Twins  are  considered 
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rather  rare  in  A1  [93] ,  although  they  have  been  reported  to  occur 
frequently  in  evaporated  films  [94]  presumably  due  to  stacking 
errors  in  the  faceted  growth  from  the  vapor.  Furthermore,  most 
A1  twins  reported  in  the  literature  show  boundaries  with  coherent 
and  incoherent  sections. 

Twins  observed  in  EHD  powders  are  definitely  formed  during 
solidification,  as  evidenced  by  the  segregation  pattern  in 
Figure  45.  Note  that  an  incipient  cell  within  the  larger  grain 
has  penetrated  the  smaller  one,  disrupting  the  otherwise  planar 
boundary.  In  those  twins  where  a  growth  direction  could  be 
identified,  the  nucleation  site  apparently  lays  on  the  junction 
of  the  twin  boundary  and  the  droplet  surface.  It  was  then 
concluded  that  twins  originate  from  a  stacking  "accident"  in  the 
nucleation  process,  and  the  boundary  so  created  propagates  later 
as  the  solid  front  advances. 

Although  it  is  anticipated  that  the  preservation  of  the  twin 
boundary  will  be  energetically  more  favorable  than  its 
degeneration  into  a  noncoherent  type,  its  apparent  perfection  is 
remarkable  in  view  of  the  high  solidification  rates  involved  in 
these  powders.  Furthermore,  twin  boundaries  are  sometimes 
preferentially  decorated  with  precipitates,  as  in  Figure  43,  and 
may  show  denuded  zones  in  their  vicinity,  as  in  Figure  41.  If  the 
boundary  were  perfectly  coherent  at  the  atomic  level,  these 
effects  could  be  explained  by  assuming  that  the  segregation 
conditions  where  the  boundary  emerges  to  the  liquid  are  different 
from  those  at  the  rest  of  the  interface.  It  could  also  be  argued 
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that  small  irregularities  occur  at  the  boundary  during  growth, 
and  their  further  rearrangement  in  the  solid  state  (to  reduce 
surface  energy)  can  accomodate  vacancies  or  solute  atoms. 

A  small  fraction  of  the  powders  exhibiting  this  phenomenon 
( < 10% )  show  more  than  one  twin  boundary  (see  Figure  46)  with  as 
many  as  4  crystallogr aphically  related  grains  being  detected  in 
some  cases.  In  addition,  twins  with  boundaries  showing 
incoherent  segments  were  also  found.  These  variations  of  the 
structure  seem  unaccountable  at  the  moment  if  we  ascribe  the 
origin  of  the  twins  only  to  the  nucleation  process. 

9.3.2  RSR  submicron  powders 

RSR  powders  exhibit  quite  different  structures  from  their 
EHD  submicron  counterparts.  They  are  not  as  spherical,  and  many 
appear  extensively  polycrystalline  with  evidence  of  a  texture 
reflected  in  the  diffraction  pattern,  e.g.  Figure  47.  However, 
only  in  a  few  instances  was  it  possible  to  distinguish  the 
individual  crystals,  and  even  under  dark  field  the  grain 
boundaries  were  not  sharply  delineated.  The  presence  of  silicon 
was  revealed  by  diffraction,  but  could  not  be  discerned  in  the 
images.  Some  of  this  Si  was  crystallogr aphically  related  to  the 
matrix  and  probably  formed  after  solidification,  see  Figure  48. 

Extensive  dislocation  substructures  were  detected  in  most 
powders,  as  opposed  to  the  EHD  material  which  is  commonly 
dislocation-free.  The  relatively  high  dislocation  density. 


Figure  46.  TEM  images  of  multiple  twins  in  an  Al-3%Si  EHD  powder. 
Two  of  the  boundaries  are  {ill}  planes  and  seem 
perfectly  coherent.  A  smaller  grain  with  partially 
faceted  boundaries  is  also  evident  but  out  of 
contrast. 


probably  stemming  from  collisions  of  the  particle  after 
solidification,  largely  contributes  to  the  rather  confusing 
aspect  of  these  powders. 
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About  the  only  conclusive  observation  at  the  present  time 

seems  to  be  the  absence  of  cellular  or  other  coarse  segregates  in 

the  examined  specimens.  Table  17  indicates  that  even  for  a  heat 

transfer  coefficient  of  2  x  10^  W  m-^  K-1-,  typical  of  RSR 

atomization,  the  external  cooling  is  expected  to  extend  the  rapid 

solidification  regime  only  slightly  over  that  of  the  EHD  powders. 

However  this  increase  is  not  enough  to  shift  the  "nose"  of  the 

recalescence  curve  beyond  any  of  the  assumed  critical 

temperatures.  The  apparent  decrease  in  segregation  is  then 

likely  to  reflect  the  achievement  of  larger  undercoolings  Atn  in 

the  RSR  process  at  comparable  particle  sizes.  This  conclusion  is 

in  line  with  the  much  larger  cooling  rates  prior  to 

8  4  -1 

solidification  estimated  for  RSR  (e.g.  ~10  vs.  ~10  K  s  for  a 
1  um  droplet) ,  and  also  with  the  higher  incidence  of  multiple 
nucleation  noted  in  these  powders.  Nevertheless,  further  work  on 
the  solidification  behavior  of  RSR  powders  is  necessary  to 
confirm  that  cooling  rates  capable  of  suppressing  the  onset  of 
nucleation  are  achievable  during  atomization. 

9.3.3  Conventional  powders 

SEM  analysis  of  RSR  powders  of  conventional  size  (>10  um) 
revealed  the  presence  of  at  least  two  distinct  structures  ,  as 
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shown  in  Figure  49.  Regular  dendritic  morphologies  with  DAS  of 
~2  pm  or  less,  as  well  as  duplex  structures,  were  observed  in  the 
Al-4.5%Cu  alloy.  The  latter  consist  of  large  u-Al  regions 
surrounded  by  a  finer  dendritic  matrix  or  eutectic.  Duplex 
configurations  were  not  detected  in  the  Al-Si  alloys,  which  were 
atomized  with  higher  initial  superheat  than  the  Al-4.5%Cu  sample. 
It  has  been  claimed  in  a  parallel  investigation  on  IN-100  alloy 
[95]  that  the  duplex  structure  results  from  partial 
solidification  of  the  metal  on  the  atomization  disc.  This 
hypothesis  was  supported  by  the  fact  that  increasing  the  initial 
melt  temperature  reduced  the  proportion  of  powders  showing  this 
particular  configuration. 

Figure  49  also  suggests  that  in  some  instances  nucleation 
might  have  occurred  inside  the  droplet,  rather  than  at  the 
surface,  but  the  evolution  of  the  structure  was  not  readily 
apparent.  Furthermore,  there  was  no  clear  indication  of  whether 
these  powders  achieved  any  substantial  undercooling  prior  to 
nucleation,  and  homogeneous  particles  were  not  observed  in  this 
size  range. 

Analysis  of  the  microstructur al  characteristics  as  a  , 

function  of  particle  size  revealed  that  the  dendrite  arm  spacings 
increase  from  1  pm  for  2rg  <  44  pm  to  about  2  pm  for 
105  <  2r q  <  125  pm,  as  shown  in  Figure  50.  This  is  a  consequence 
of  the  increase  in  solidification  time  (or  reduction  in  cooling 
rate)  for  the  larger  particles.  In  addition,  decreasing  the 
solute  content  reduces  the  amount  of  segregation  but  does  not 
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significantly  change  the  scale  of  the  microstructure.  Finally, 
some  conventional  EHD  powders  show  finer  segregate  spacings  than 
RSR  material  of  equivalent  size.  This  finding  is  not  consistent 
with  the  lower  cooling  rates  expected  from  radiation,  and 
probably  means  that  the  particles  analized  solidified  in  contact 
with  a  substrate. 


Table  18 


Effect  of  Particle  Size  on  the  Lattice  Parameter 


t 


of  the  Primary  Phase  in  Al-Si  Alloys 


Diameter  (pm)  Al-3%Si 


Al-6%Si 


<44 

4.0445 

45 

-  74 

4.0453 

75 

-  105 

4.0460 

105 

-  125 

4.0462 

4.0439 

4.0447 

4.0450 

4.0454 


t  Average  of  [220],  [331]  and  [222]  peaks. 


The  effect  of  particle  size  on  the  potential  extension  in 
solid  solubility  was  investigated  by  X-ray  dif fractometry .  The 
results,  shown  in  Table  18,  indicate  that  reducing  the  particle 
size  (hence  increasing  the  achievable  interface  velocities, 
regardless  of  geometry)  reduces  the  lattice  parameter,  aQ,  and 
increases  the  solute  content  in  the  primary  phase.  In  addition, 
increasing  the  nominal  alloy  composition  also  enhances  the 
solubility  of  Si.  According  to  previous  research  by  Itagaki  et 
al.  [33] ,  an  increase  of  1%  in  Si  content  shrinks  the  lattice 
parameter  approximately  by  .001  A.  Therefore,  it  was  estimated 


that  an  extra  1.3%Si  was  dissolved  in  A1  by  going  from  the 
coarser  (>100  pm)  to  the  finer  (<44  urn)  powders.  It  should  be 
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noted,  however,  that  although  the  trend  in  aQ  vs.  rQ  is 
consistent  with  that  expected  from  heat  flow  considerations,  the 
changes  detected  are  in  the  order  of  magnitude  of  the 
experimental  error.  Nevertheless,  the  approach  is  conceptually 
sound,  and  with  careful  implementation  may  provide  a  reasonable 
way  of  relating  microstructural  modifications  to  heat  flow 
predictions. 
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CHAPTER  10 

SUMMARY  AND  CONCLUSIONS 

The  infinite-series  analytical  solution  to  the  problem  of 
radially  uniform,  unidimensional  heat  flow  in  an  initially 
superheated  metal  droplet  was  coupled  to  numerical  formulations 
based  on  the  Enthalpy  model  for  studying  their  solidification 
behavior.  The  Newtonian  model  was  reviewed  and  extended  to  use 
it  as  a  basis  of  comparison  and  to  develop  trends  for  the 
limiting  case  where  the  resistance  to  heat  flow  inside  the 
droplet  is  negligible. 

The  interface  was  assumed  morphologically  stable  under  all 
conditions,  and  the  analysis  covered  two  basic  growth  modes. 
First,  isothermal  solidification  at  the  melting  temperature  (no 
kinetic  barrier  or  undercooling)  with  concentric  isotherms.  In 
the  second  mode,  undercooling  is  necessary  for  nucleation  and 
growth.  Solidification  starts  at  a  point  on  the  droplet  surface 
and  follows  an  axisymmetric  geometry.  Although  most  of  the  work 
was  focused  on  powders  growing  from  a  single  nucleation  event, 
the  implications  of  multiple  nucleation  were  also  discussed  for 
the  case  of  microcrystalline  structures. 
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10.1  Liquid  Cooling  Stage 

Analysis  of  the  pre-solidification  stage  using  the  Newtonian 
and  series  solutions  resulted  in  the  following  conclusions: 

1.  Both  models  predict  that  the  instantaneous  average 
cooling  rate  in  the  liquid  is  directly  proportional  to 
the  heat  transfer  coefficient,  h,  inversely  proportional 
to  the  droplet  radius  Tq,  and  independent  of  initial 
superheat  6q. 

2.  The  cooling  time  estimated  from  the  analytical  solution 
is  always  lower  than  the  Newtonian  prediction  for  a 
fixed  surface  temperature  (0g  =  0N)  ,  and  higher  for  a 
required  enthalpy  content  in  the  droplet  prior  to 
nucleation  <PN.  tN  increases  with  r0  and  0Q,  but 
decreases  with  increasing  h  and  0N. 

3.  Reaching  a  desired  value  becomes  more  difficult  as 
the  droplet  size  increases,  not  only  because  nucleation 
must  be  suppressed  to  lower  temperatures,  but  also  due 
to  the  consequent  decrease  in  cooling  rate. 

4.  As  expected,  the  two  models  converge  with  decreasing 
Biot  number,  but  substantial  temperature  gradients  in 
the  droplet  were  estimated  for  Bi  as  low  as  0.01. 

10.2  Solidification  without  Undercooling 

The  trends  obtained  from  the  Enthalpy  and  Newtonian  models 
for  solidification  without  significant  undercooling  are 


summarized  below: 
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Reducing  diminishes  the  net  solidification  time, 
tf,  in  a  more  pronounced  way  than  the  direct 
proportionality  predicted  by  the  Newtonian  model. 
Furthermore,  it  increases  the  interface  velocity  R  at 
equivalent  fractions  solidified  towards  the  limiting 
Newtonian  values  (which  are  independent  of  rQ) . 

Reducing  Tq  also  decreases  the  G^/R  ratios  and  increases 
the  achievable  cooling  rates  in  the  liquid  (eL)  at  the 
beginning  of  solidification,  and  in  the  solid  (eg)  at 
the  end  of  it. 

Increasing  h  reduces  t^  and  increases  R  in  a  less 
pronounced  manner  than  that  predicted  by  the  Newtonian 
model.  It  always  increases  Gl/r  as  well  as  the  cooling 
rates  at  any  stage  during  solidification. 

Increasing  initial  superheat  increases  tf  GL/R  and 
eL,  but  decreases  R  and  eg.  The  effects  of  superheat  are 
mostly  felt  in  the  initial  stages  of  solidification  and 
decay  rapidly  with  decreasing  Biot  number. 

The  relatively  low  cooling  rates  calculated  from  the 
enthalpy  model  are  due  in  part  to  the  prescription  of  a 
discrete  melting  temperature.  Alloys  solidifying  under 
similar  conditions  are  expected  to  experience  larger 
cooling  rates  since  the  process  takes  place  through  a 
liquid-solid  range. 
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10.3  Solidification  with  Undercooling 

The  dimensionless  enthalpy-temperature  diagram  provides  a 
convenient  way  to  describe  and  compare  the  thermal  history  of 
undercooled  droplets  during  solidification.  The  approach  is 
particularly  useful  in  the  Newtonian  analysis,  which  leads  to  the 
following  results: 

1.  The  thermal  history  of  the  droplet  is  independent  of 
particle  size  r^,  for  a  given  growth  geometry,  and  can 
be  readily  characterized  by  two  parameters:  the  initial 
undercooling  0^,  and  the  ratio  of  a  dimensionless 
kinetic  coefficient  to  the  Biot  number,  (K/Bi) .  The 
latter  reflects  the  degree  of  predominance  of  the 
recalescence  rate  (release  of  latent  heat)  over  the 
external  cooling. 

2.  The  thermal  history  normally  consists  of  two  regimes. 

The  first  one  is  the  recalescence  or  "rapid 
solidification"  stage,  where  the  undercooled  droplet  is 
absorbing  most  of  the  released  heat  of  fusion,  and  the 
heat  loss  to  the  environment  plays  only  a  minor  role. 

The  second  regime  is  one  of  slower  growth,  roughly 
isothermal,  and  is  controlled  by  the  external  rate  of 
heat  extraction. 

3.  For  typical  values  of  K/Bi  encountered  in  atomization 
(K/Bi>100) ,  the  external  cooling  has  little  or  no  effect 
in  extending  the  "rapid  solidification"  regime  to  a 
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fraction  solid  substantially  larger  than  that  dictated 
by  the  limiting  case  of  adiabatic  solidification, 
ga  =  "Ojj.  Increasing  the  number  of  nucleation  events 
further  reduces  the  role  of  the  external  cooling  during 
recalescence . 

4.  Although  the  net  solidification  time,  t^,  decreases  with 
increasing  undercooling  prior  to  nucleation,  the  total 
time  spent  by  the  droplet  in  cooling  and  solidification 
increases,  as  much  as  twice  from  the  case  of 
no-undercooling,  thereby  requiring  longer  flight  times 
in  a  given  atomizer.  As  before,  the  times  increase  with 

Cg  and  decrease  with  larger  h  values. 

The  Enthalpy  model  for  undercooled  solidification  revealed 
that  departures  from  the  Newtonian  model  may  occur  at  Biot 
numbers  as  low  as  10“^.  This  is  a  direct  consequence  of  the 
"local"  recalescence  of  the  interface  which  is  now  able  to  alter 
both  its  velocity  and  the  thermal  field  in  its  vicinity  in  order 
to  accomodate  tha  latent  heat  released.  In  general,  this  model 
predicts  that: 

1.  For  a  given  undercooling,  the  local  recalescence  of  the 
interface  is  accelerated  (and  hence  the  extent  of  rapid 
solidification  reduced)  as  rQ  and  the  kinetic  parameter 
Km  are  increased. 

2.  Interface  velocities  and  undercoolings  at  any  fraction 
solid  are  increased  by  reducing  rQ  or  increasing  0N  and 
h,  although  this  latter  variable  becomes  important  only 
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10.4 


towards  the  end  of  recalescence.  Furthermore, 
increasiang  will  reduce  the  interface  undercooling 
but  increase  its  velocity. 

3.  Both  models  indicate  that  if  the  critical  temperature, 
Tc,  for  a  certain  metastable  effect  is  reasonably  close 
to  Tm,  rapid  solidification  may  be  obtained  with 
moderate  initial  undercoolings  provided  the  K/Bi  ratio 
is  small  enough  to  make  AT*^n  <  AT^.  However,  if 

Tc  «  tm  the  extent  of  rapid  solidification  is 
limited  by  0^.  in  any  case,  the  Newtonian  model  always 
provides  a  limit  to  the  achievable  fraction  solidified 
below  Tc. 

4.  Maximum  temperature  differences  inside  the  droplet  (i.e. 
departure  from  Newtonian  conditions)  increase  with  Biot 
number,  KM  and  Atn. 

5.  The  interface  shapes  predicted  by  the  enthalpy  model  are 
close  to  the  convex  spherical  front  configuration 
assumed  for  the  Newtonian  analysis.  At  large  Bi  and 

Km,  and  small  0N  the  interface  becomes  concave  towards 
the  liquid,  approaching  the  concentric  geometry  of  the 
no-undercooling  case. 

Microstructural  Observations 


Examination  of  rapidly  solidified  aluminum  alloy  powders 
produced  by  two  different  atomization  techniques  (EHD  and  RSR) 
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revealed  some  novel  microstructures,  especially  in  the  submicron 
size  range.  The  results  were  coupled  with  heat  flow  concepts  in 
developing  an  understanding  of  the  thermal  history  during 
solidification.  The  relevant  findings  are  listed  below: 

1.  Decreasing  the  powder  size,  therefore  increasing  the 
achievable  undercoolings  and  interface  velocities, 
reduces  the  extent  of  segregation.  This  translates  into 
an  increased  solute  content  of  the  primary  phase, 
detected  in  the  conventional  size  powders,  and 
eventually  in  homogeneous  solidification  of  submicron 
droplets  with  Si  contents  as  high  as  6%. 

2.  Increasing  h  from  radiation  (EHD)  to  convection  (RSR) 
appears  to  reduce  segregation,  presumably  by  suppresing 
nucleation  towards  larger  undercoolings  due  to  an 
increase  of  4  orders  of  magnitude  in  the  calculated 
cooling  rate. 

3.  Submicron  powders  showing  segregation  but  not  cellular 
breakdown  may  be  evidence  of  interface  stability  due  to 
surface  tension  effects. 

4.  Increasing  undercooling  by  reducing  particle  size 
promotes  multiple  nucleation  and  the  incidence  of  twin 
crystals  in  the  submicron  range.  The  latter  are 
believed  to  originate  mainly  from  stacking  "accidents" 
in  the  nucleation  process.  The  twin  boundaries  appear 
to  be  coherent,  although  a  few  instances  of  incoherent 
segments  were  observed.  In  general,  twins  and 
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polycrystalline  powders  are  less  segregated  than  their 
monocrystal  counterparts. 
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APPENDIX  A 


METAL  PROPERTIES  USED  IN  THE  CALCULATIONS* 


Al 

Fe 

Ni 

Melting  Temperature 
tm  <k> 

933 

1810 

1726 

Hypercooling  Temperature 
th  (K) 

569 

1464 

1315 

Heat  of  Fusion 

AHM  (10 5  j  kg-1) 

3.95 

2.72 

3.02 

Average  Density 
p  (kg  m-3) 

2.46 

7.19 

8.07 

Specific  Heat 

(10 ^  j  kg-1K"1)CT 

CS 

1.09 

1.19 

0.787+ 

0.764+ 

0.734 

0.619 

Thermal  Conductivity 
(W  m-1  IT1)  kT 

kS 

90.8 

210.8 

40.3 

34.6 

39.3  + 
87.0  + 

Thermal  Diffusivity 
(10-5m2  s-1)  aL 

aS 

3.4 

7.18 

0.712 

0.630 

0.664 

1.74 

Stefan  Number 

Ste  (Tq  =  300) 

1.74 

4.37 

3.47 

*  Thermophysical  properties  at  T^ 
t  Estimated  by  extrapolation. 
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APPENDIX  B 

NOMENCLATURE 


Dimensioned  variables 

A  heat  transfer  area  between  two  neighbor  volume  elements 

aQ  lattice  parameter  of  the  aluminum 

C  specific  heat 

DAS  dendrite  arm  spacing 

self-diffusion  coefficient  in  the  liquid 

dMX  distance  between  two  nearest  neighbors  in  a  close  packed 
microcrystalline  structure 

Ed  activation  energy  of  diffusion 

e  emissivity 

F  free  energy 

G  temperature  gradient  at  the  liquid  solid  interface 
g  fraction  solid 

H  enthalpy  content  over  that  of  the  solid  at  the  melting 
temperature  HgM  ^ 

Hm  heat  of  fusion 

h  heat  transfer  coefficient 

J  nucleation  frequency 

Km  kinetic  parameter  relating  the  velocity  and  undercooling  in 
the  linear  law,  evaluated  at  the  melting  temperature 

k  thermal  conductivity 

m  metric  coefficients 


f 
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N 

Q 

R 

r 

T 

Tc 

T0 

AT 

A  (T) 

t 

V 

V 


z 


a 

6 


Y 


v 


► 1 

a 

1 1 


number  of  nuclei  or  grains 

heat  flow  term  in  the  enthalpy  balance  equations 

interface  velocity 

radial  coordinate,  0  <_  r  <_ 

temperature 

critical  temperature  for  a  metastable  effect  to  be  achieved 
by  rapid  solidification 

initial  temperature 

undercooling 

maximum  temperature  difference  inside  the  droplet  during 
solidification 

time 

volume 

velocity  of  the  droplet  relative  to  the  cooling  environment 
(gas) 

distance  from  the  center  of  the  droplet  along  the  growth 
axis 

thermal  diffusivity 
length  of  an  interface  step 
cooling  rate 

average  jumping  distance  for  an  atom  to  join  the  solid 

kinematic  viscosity 

density 


Stefan-Boltzmann  constant 
jumping  frequency  of  an  atom 
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Dimensionless  variables 

Bi  Biot  number,  hrQ/|<L 

2 

Fo  Fourier  number,  at/r^j 

K  Kinetic  parameter,  (rg/aL)  (AHM/CL)  Km 

M  Mach  number,  V00/Vsound 

Nu  Nusselt  number,  2hrg/kG 

Pe  Peclet  number,  Rr./2aT 

t'  L 

Pr  Prandtl  number,  v/a 

R  interface  velocity,  (rg/aL)  R 
Ste  Stefan  number,  Cl(tm  -  TG) /AHM 

u  coordinate  in  the  bispherical  system  perpendicular  to  the 
growth  direction 

v  coordinate  in  the  bispherical  system  along  the  growth 
direction 

Z*  position  of  the  interface  along  the  growth  axis,  z/rQ 
6  Cahn's  parameter  in  the  kinetic  equations 

0  temperature,  Cl(t  - 

<})  radial  coordinate,  rAg 

,p  enthalpy,  (H  -  HgM)/AHM 

Subscr ipts 

c  convection 

E  equilibrium 

f  end  of  solidification 

G  gas  environment 


t 
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H  hypeccooling 

i  nodal  index  in  the  radial  coordinate  for  the  case  of 

no-undercooling  and  in  the  u  direction  for  the  bispherical 
mesh. 

j  nodal  index  in  the  v  direction  for  the  bispherical  mesh 
k  time  step  index 

l  denotes  the  node  containing  the  liquid-solid  interface 
L  liquid 

M  at  the  melting  temperature 

m  number  of  nodes  in  the  v  direction 
N  at  the  moment  of  nucleation 

n  number  of  nodes  in  the  <t>  direction  for  the  concentric  model 

and  in  the  u  direction  for  the  bispherical  mesh 

r  radiation 

S  solid 

s  at  the  surface 

t  at  the  dendrite  tip 


Superscripts 
*  at  the  interface 

+  in  the  positive  r  direction  for  the  concentric  geometry 
in  the  negative  r  direction 

N  in  the  positive  u  direction  for  the  bispherical  mesh 

S  in  the  negative  u  direction 

E  in  the  positive  v  direction 

in  the  negative  v  direction 


W 


